|
|
||||||||
Cancer Diagnostics |
1 Department of Medical Oncology, Erasmus Medical Center, Rotterdam, The Netherlands.
2 Department of Veterinary Pathobiology, Faculty of Life Sciences, University of Copenhagen, Copenhagen, Frederiksberg C, Denmark.
aAddress correspondence to this author at: Erasmus MC, Josephine Nefkens Institute, Rm. Be 400, PO Box 2040, 3000 CA Rotterdam, The Netherlands. Fax 31-10-408-8377; e-mail a.sieuwerts{at}erasmusmc.nl.
| Abstract |
|---|
|
|
|---|
Methods: We measured mRNA transcripts of full-length TIMP1 (TIMP1-v1) and the novel splice variant lacking exon 2 (TIMP1-v2) by use of real-time RT-PCR in frozen primary tumor samples. Transcript concentrations are correlated with histomorphological and biological factors, TIMP-1 protein, and distant metastasis-free survival (MFS) and OS time.
Results: TIMP1-v1 and TIMP1-v2 alone were not informative with respect to predicting prognosis. However, the PCR assay designed to measure the combination of v1 + v2 showed that high concentrations of this combination were associated with good prognosis. In Cox multivariate regression analysis, which also included the traditional prognostic factors, increasing concentrations were independently associated with prolonged MFS (P = 0.004) and OS (P = 0.048). Including TIMP-1 protein and TIMP1-v1+v2 mRNA together in the multivariate model revealed that protein and mRNA were both independently associated with prognosis, with hazard ratios pointing in opposite directions.
Conclusion: High concentrations of TIMP1-v1+2 mRNA are associated with good prognosis in patients with primary breast cancer. Since high concentrations of TIMP-1 protein are associated with poor prognosis, the presence of possible posttranscriptional mechanisms requires further investigation.
| Introduction |
|---|
|
|
|---|
In analogy with free and complexed TIMP-1 protein, the prognostic value of TIMP-1 may be improved by detection of specific splice variants of TIMP1 mRNA. Furthermore, biological understanding of TIMP-1 protein and its gene might help in understanding the controversial findings about TIMP-1 with respect to tumor development and prognosis. To address this, we analyzed mRNA concentrations of the common full-length variant of TIMP1 (v1) and a newly discovered splice variant (v2) lacking exon 2 in a large cohort of 1301 primary breast tumors. We related TIMP1-v1 and v2 expression with histomorphological and clinical factors, mRNA expression of the proliferation marker Ki-67, and total TIMP-1 protein concentrations. Finally, we investigated whether mRNA expression of the TIMP1 splice variants adds to the prognostic value of total TIMP-1 protein.
| Patients and Methods |
|---|
|
|
|---|
The remaining 1301 patients were treated either with breast-conserving surgery (44%) or with modified mastectomy (56%); 931 patients (72%) received adjuvant radiotherapy. During this period, 195 of the 620 lymph nodepositive (LNP) patients did not receive adjuvant systemic therapy; 425 of the LNP patients were treated with adjuvant systemic therapy, of these patients 197 received hormonal therapy, 210 chemotherapy, and 18 received combination therapy. Routine postsurgical follow-up and definition of time to metastasis were as described (20). Median follow-up was 92 months (range 3248 months). Six hundred sixty-nine (51%) patients developed a distant metastasis and count as events in the analysis for metastasis-free survival (MFS). Seventy-two patients (6%) died without evidence of disease and were censored at last follow-up in the analysis of MFS. Five hundred twenty-six patients (40%) died after a previous relapse. Thus, 598 patients (46%) were counted as events in the analysis of overall survival (OS). Tumor staging was according to the Union Internationale Contre le Cancer tumor node metastasis classification. Other relevant patient and tumor characteristics are listed in Table 1
(see also Table 1 in the Data Supplement that accompanies the online version of this article at http://www.clinchem.org/content/vol53/issue7).
|
tissue processing and estimation of the amount of invasive tumor cells
We processed tissue and estimated amount of invasive tumor cells as described (21)(22). Only specimens with at least 30% of the nuclei of epithelial tumor cell origin and distributed uniformly over at least 70% of the hematoxylin-eosinstained tissue section area were included. Furthermore, we dichotomized our tumor cohort at the median of 70% tumor cell nuclei in stroma-rich tumors (primary tumors containing
30% stromal components) and stroma-poor tumors (primary tumors containing at least 70% tumor cells).
rna isolation, CDNA synthesis, and quantification of specific MRNA species
RNA isolation, cDNA synthesis, quantification of specific mRNA species, and quality control checks were done as described in detail (21). We performed real-time RT-PCR in an ABI Prism 7700 Sequence Detection System (Applied Biosystems) and a Mx3000PTM Real-Time PCR System (Stratagene) using both Assay-on-Demand from Applied Biosystems and the intron-spanning forward and reverse primer combinations at the conditions shown in Supplemental Data Table 2. RT-PCR products generated with these primers by cell lines and cancer tissue samples were gel-purified and sequenced (AGOWA, Berlin, Germany), and sequences were confirmed by BLAST search as described (23). Primer sequences for ESR1, PGR, and the housekeeping genes have all been described, as have the PCR reactions and validations performed to ensure PCR specificity (21). To measure concentrations of the proliferation marker Ki-67, we used the Hs00606991 m1 Assay-on-Demand from Applied Biosystems. Concentrations of the target genes, expressed relative to our housekeeping set [low-abundance housekeeping gene hydroxymethylbilane synthase (HMBS, formerly porphobilinogen deaminase, PBGD), medium-abundance hypoxanthine phosphoribosyltransferase (HPRT), and high-abundance ß2-microglobulin (B2M)], were quantified as follows: mRNA target = 2(mean Ct housekeeping mean Ct target), as described (21).
esr1 and pgr MRNA receptor status
We established mRNA cutpoints to define tumors as steroid hormone positive at 0.2 for estrogen receptor (ER) and 0.1 for progesterone receptor (PGR). We compared these mRNA cutoffs with the established protein cutpoints of 10 fmol/mg protein in the 1203 samples with known protein concentrations as established by ELISA. For ER, sensitivity and specificity were 93% and 72%, respectively, and positive and negative predictive accuracy were 90% and 81%. For PGR, sensitivity and specificity were 84% and 83%, respectively, and positive and negative predictive accuracy rates were 90% and 75%.
elisa
To compare TIMP1 mRNA concentrations with TIMP-1 protein, we used protein concentrations that were previously measured in cytosol preparations of the same tumors (6).
statistics
Computations were done with the use of the STATA statistical package, release 9 (STATA). We assessed differences in concentrations with the MannWhitney U-test or KruskalWallis test. In these tests, patient and tumor characteristics were used as grouping variables. We tested the strengths of the associations between continuous variables with the Spearman rank correlation (rs). Variables were either log-transformed or Box-Coxtransformed to reduce the skewness. We investigated the prognostic values of the clinical and biological variables, with MFS and OS as the endpoints in the univariate, multivariate, and interaction analyses, with the use of the Cox proportional hazards model. The hazard ratio (HR) and its 95% CI were derived from these results. The proportionality assumption was investigated with a test based on the Schoenfeld residuals. We used Kaplan-Meier survival plots and log-rank tests for trend to assess the differences in time of the predicted high-risk, intermediate-risk, and low-risk groups of patients. All tests were 2-sided, and P <0.05 was considered statistically significant.
| Results |
|---|
|
|
|---|
|
All tumors readily expressed TIMP1-v1 mRNA; for only 16 of 1301 tumors were we unable to detect v2 transcripts within 45 amplification rounds. Because TIMP-1 protein overexpression has been inversely associated with cell proliferation (8), we matched our TIMP1 PCR data with those of the proliferation marker Ki-67 measured in the same preparations. The strength of the associations between TIMP1-v1, TIMP1-v2, TIMP1-v1+2, and Ki-67 mRNA are summarized in Table 3
. Whereas TIMP1-v1+2 showed a statistically significant association with TIMP1-v1, no correlation was observed with TIMP1-v2 (Spearman rs = 0.51 and 0.04, respectively). TIMP1-v1+2, compared with TIMP1-v1 and TIMP1-v2 separately, showed the strongest (inverse) correlation with Ki-67 (rs = 0.31, 0.09, and 0.01, respectively).
|
associations with histomorphological and clinical factors
Associations of mRNA expression of TIMP1-v1, TIMP1-v2, and TIMP1-v1+2 with patient and tumor characteristics are shown in Table 1
. Most notable are the opposing results of the TIMP1-v2 assay with the TIMP1-v1 and TIMP1-v1+2 assays. Whereas concentrations assessed with the TIMP1-v2 assay were negatively associated with age, and were lower in postmenopausal patients, results were the opposite for the other 2 TIMP1 assays. Furthermore, only for the TIMP1-v2 assay were concentrations significantly lower in the group of LNN patients compared with the group of LNP patients, and higher in the stroma-rich tumors compared with the group of stroma-poor tumors. On the other hand, concentrations measured with the TIMP1-v1 and TIMP1-v1+2 assays were significantly higher in ER-positive tumors compared with ER-negative tumors. No such difference was observed for concentrations measured with the TIMP1-v2 assay. Also with respect to grade and tumor size, only the TIMP1-v1 and TIMP1-v1+2 assays showed a relationhigher concentrations in the prognostically more favorable tumors. Finally, for all TIMP1 variants, concentrations were higher in infiltrating lobular carcinoma compared with infiltrating ductal carcinoma.
univariate and multivariate analysis for mfs and os
To assess a possible relationship of TIMP1 mRNA with prognosis, we first performed Cox univariate analyses for MFS and OS as a function of continuous TIMP1 mRNA concentrations. Because the proportional hazards assumption was violated for the total follow-up time, we restricted our exploration of the relationships of TIMP1 with MFS to the first 5 years of follow-up, with 571 failures in the 1301 patients. In these analyses, only concentrations measured with the TIMP1-v1+2 assay were significantly associated with MFS (HR 0.80, P <0.001) and OS (HR 0.85, P <0.001). Next, the TIMP1 variants were separately introduced to the base multivariate model that included the factors of age, menopausal status, nodal status, tumor size, grade, ER, and PGR (Table 4
). Of the 3 TIMP1 assays, only the assay measuring both variant 1 and variant 2 (TIMP1-v1+2) contributed significantly to the multivariate model for MFS (HR 0.89, P = 0.004) (Table 4
) and OS (HR 0.92, P = 0.048) (data not shown in a table). Adding Ki-67 and adjuvant systemic therapy or radiotherapy to the multivariate model that included TIMP1-v1+2 did not alter the coefficients of TIMP1-v1+2. To visualize the prognostic value of TIMP1-v1+2 in Kaplan-Meier curves, we divided mRNA concentration curves into 3 equal parts (low, intermediate, and high). These curves are shown in Fig. 1A
for MFS and Fig. 1B
for OS.
|
|
nodal and er status
We next performed exploratory Cox univariate analyses for MFS and OS as a function of mRNA expression in the clinically relevant subgroups of LNN, LNP, ER+, and ER patients. Only those analyses that gave significant results after concentrations were entered as transformed continuous variables are shown in Table 5
. TIMP1-v1 mRNA concentrations alone were not significantly associated with nodal status and steroid hormone receptor status. But the following 2 divergent observations between the various assays are notable. First, whereas concentrations measured with the TIMP1-v2 assay were associated with good prognosis only in the subgroup of LNP patients, the association of increasing concentrations measured with the TIMP1-v1+2 assay and good prognosis were independent of nodal status. Second, whereas increasing concentrations of TIMP1-v1+2 were associated with good prognosis exclusively in the group of ER+ tumors, concentrations measured with the TIMP1-v2 assay were associated with good prognosis only in the group of ER tumors.
|
correlations between timp-1 Mrna and protein
To compare TIMP1 mRNA with total TIMP-1 protein, we made use of protein concentrations that were previously measured with ELISA in cytosol preparations of the same tumors (6). In the 839 tumors with both measures, the highest correlation between total protein and mRNA was observed for the real time RT-PCR assay able to measure TIMP1-v1 (rs = 0.34, P <0.001). The strength of the association was lower for the PCR assay able to measure both variant 1 and 2 (TIMP1-v1+2) (rs = 0.28, P <0.001) and inverse for the TIMP1-v2 assay (rs = 0.11, P = 0.01).
To ensure that our cohort of 1301 patients did not differ from the cohort of 2984 patients with protein data, we repeated all analyses for the overlapping cohort of 839 patients. We divided the protein concentrations in 3 equal parts to classify the tumors at the protein level as TIMP-1 low, intermediate, or high. High and intermediate vs low concentrations of TIMP-1 protein were significantly associated with shorter MFS in univariate analysis (HR intermediate vs low 1.56, HR high vs low 1.37, P = 0.002) and multivariate analysis (HR intermediate vs low 1.42, HR high vs low 1.42, P = 0.010), which is in agreement with the original study, in which high tumor tissue concentrations of TIMP-1 protein were identified as an independent marker of poor prognosis in 2984 primary breast cancers (6).
Next, we similarly compared the prognostic value of mRNA concentrations measured with our TIMP1-v1, v1+2, and v2 PCR assays in these 839 patients. Only increasing concentrations of TIMP1-v1+2 mRNA were significantly associated with a prolonged MFS in univariate analysis (HR intermediate vs low 0.81, HR high vs low 0.59, P <0.001) and multivariate analysis (HR intermediate vs low 0.89, HR high vs low 0.66, P = 0.020).
Finally, we explored a potential interaction between TIMP-1 protein and TIMP1 mRNA with respect to MFS. No such interaction was observed (P = 0.56). Including log-transformed continuous concentrations of both TIMP-1 protein and TIMP1-v1+2 mRNA to the base multivariate model for MFS revealed that TIMP-1 protein and TIMP1-v1+2 mRNA were both independently associated with prognosis, with HRs pointing in opposite directions (HR 1.36, 95% CI 1.161.61, P <0.001, n = 839 for protein and HR 0.78, 95% CI 0.690.87, P <0.001, n = 839 for mRNA).
| Discussion |
|---|
|
|
|---|
Our main purpose was to investigate the potential prognostic value of TIMP1 mRNA and a newly discovered splice variant to gain more knowledge on the biology of TIMP-1 in breast cancer. To address this, our retrospective study that included RNA preparations from tumor tissue obtained from 1301 patients suffering from primary breast cancer is, to the best of our knowledge, the largest study performed on the mRNA concentrations of TIMP1 to date. Because of different assay conditions, absolute values of real-time RT-PCR assays can be compared only within an assay, and values from different assays, such as our 2 assays measuring individual transcripts of TIMP-v1 and TIMP-v2, cannot simply be added. With our multiplex TIMP1-v1+2 assay that measures both transcripts in the same reaction with the same primer pairs, we corrected as much as possible for such differences in assay conditions. However, by measuring both transcripts in the same reaction, we cannot exclude that the shorter (v2) variant was favored relative to the larger (v1) variant.
All tumors expressed full-length TIMP1-v1 mRNA; we were unable to detect v2 transcripts for only 16 of 1301 tumors. Sequence analysis of TIMP1-v2 has already revealed that this variant lacking exon 2, if translated, is probably a soluble, intracellular protein lacking part of the region that directs the main inhibitory MMP activity (unpublished data). It is therefore unlikely that a putative protein of TIMP1-v2 forms complexes with MMPs, and thus it probably exhibits a biological function different from full length TIMP-1indeed, our data suggest that. We found that TIMP1-v1 mRNA concentrations increase with age, are higher in ER/PGR-positive tumors, and are higher in smaller-sized and moderately differentiated to well differentiated tumors. In contrast to some reports (3)(4)(9), but in agreement with another report (10), TIMP1-v1 mRNA concentrations were not different in our cohort of 620 LNP patients compared with the group of 681 LNN patients. These discrepancies might be due to the relatively small sample sizes in the earlier studies [n = 30 LNN and 24 LNP (3); n = 49 LNN and 66 LNP (4)].
Separate evaluation of TIMP1-v2 in our patient cohort revealed a strong inverse correlation with age and no correlation with ER, PGR, grade, and tumor size. In addition, only for TIMP1-v2, concentrations were higher in stroma-rich compared with stroma-poor primary breast tumors. The lack of a correlation between ER and TIMP1-v2 concentrations suggests that TIMP1-v2, unlike TIMP1-v1, is regulated by an ER-independent mechanism. Moreover, the lower TIMP1-v2 and higher TIMP1-v1 mRNA concentrations in the older age group support our hypothesis that v2 is regulated by a different mechanism. Another observation we made is the relatively strong negative correlation between the proliferation marker Ki-67 and TIMP1-v1+2. No such correlation was observed for the TIMP1 assays able to measure the variants separately. This finding suggests that only the combined action of full-length TIMP-1 and its del-2 variant are effectively able to downregulate proliferation or monitor-reduced proliferation.
We recently raised the hypothesis that high concentrations of total TIMP-1 protein are not necessarily associated with poor prognosis but that the association depends on the ratio of uncomplexed/total TIMP-1 (19). In analogy with this, our present study shows that TIMP1-v1 mRNA and TIMP1-v2 mRNA alone were not associated with prognosis. However, our real-time RT-PCR assay developed to measure both transcripts at the same time revealed that high mRNA concentrations of the combination of both variants were associated with low tumor aggressiveness. Whether changing the balance between full-length TIMP-1 and its variant lacking exon 2 has potential as a possible therapeutic approach to reduce tumor aggressiveness remains to be investigated. To establish this, and since it is only the actual protein that is biologically active, variant-specific immunohistochemistry and a quantitative assay (ELISA) able to measure the putative del-2 protein in relation to full-length TIMP-1 protein are required.
In conclusion, this retrospective study on a large cohort of primary breast cancers provides evidence that the combined expression of full-length TIMP1-v1 mRNA and its v2 variant lacking exon 2 are associated with low tumor aggressiveness. This splice variant-dependent association might help our understanding of the role of TIMP-1 with respect to breast cancer.
| Acknowledgments |
|---|
Financial disclosures: None declared.
Acknowledgments: This work is the result of an EORTC-PathoBiology Group collaboration. We thank Miranda Arnold, Anneke Goedheer, Roberto Rodriguez-Garcia, Anita Trapman, Vanja de Weerd, and Henk Portengen for their technical support.
| Footnotes |
|---|
2 Human genes: TIMP1, TIMP metallopeptidase inhibitor 1; TIMP1-v1, TIMP1 full-length variant; TIMP1-v2, TIMP1 variant lacking exon 2; ESR1, estrogen receptor 1 (formerly ER-
); PGR, progesterone receptor; HMBS, hydroxymethylbilane synthase, formerly porphobilinogen deaminase, PBGD; HPRT, hypoxanthine phosphoribosyltransferase; B2M, ß-2-microglobulin. ![]()
| References |
|---|
|
|
|---|
The following articles in journals at HighWire Press have cited this article:
![]() |
T. A. O'Mara, J. A. Clements, and A. B. Spurdle The Use of Predictive or Prognostic Genetic Biomarkers in Endometrial and Other Hormone-Related Cancers: Justification for Extensive Candidate Gene Single Nucleotide Polymorphism Studies of the Matrix Metalloproteinase Family and their Inhibitors Cancer Epidemiol. Biomarkers Prev., September 1, 2009; 18(9): 2352 - 2365. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. van Agthoven, A. M. Sieuwerts, M. E. Meijer-van Gelder, M. P. Look, M. Smid, J. Veldscholte, S. Sleijfer, J. A. Foekens, and L. C.J. Dorssers Relevance of Breast Cancer Antiestrogen Resistance Genes in Human Breast Cancer Progression and Tamoxifen Resistance J. Clin. Oncol., February 1, 2009; 27(4): 542 - 549. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |