Desmoplasia measured by computer assisted image analysis: an independent prognostic marker in colorectal carcinoma
- 1Department of Pathology, Dokuz Eylül University School of Medicine, Inciraltı, Izmir 35340, Turkey
- 2Department of General Surgery, Dokuz Eylül University School of Medicine
- 3Department of Biomedical Engineering, Dokuz Eylül University School of Medicine
- Correspondence to: Dr S Sarioglu Dokuz Eylül University School of Medicine, Department of Pathology, Inciraltı, Izmir 35340, Turkey;
- Accepted 27 July 2004
Aims: The assessment of desmoplasia by traditional semiquantitative methods does not provide reliable prognostic data. The aim of this study was to quantify desmoplasia by computerised image analysis in primary colorectal carcinomas and to investigate its ability to predict overall survival.
Methods: In total, 112 colorectal adenocarcinomas, with a median follow up of 66 months, were studied. The representative tumour sections were stained by the van Gieson method, which stains collagen rich stroma red. For quantitative histochemical measurement, digital images were analysed by a computerised image analysis program to calculate the percentage of red stained tissue area. The percentage of desmoplasia (PD) was related to conventional clinicopathological prognostic factors and overall survival.
Results: The mean (SD) PD was 4.85 (3.37). PD was found to be significantly associated with lymph vessel and venous invasion. By Kaplan–Meier analysis, PD was associated with survival—patients with PD > 4 had a shorter survival than those with PD ⩽ 4. In multivariate analysis, tumour stage, distant metastasis, and PD emerged as independent prognostic factors.
Conclusion: Desmoplasia measured by image analysis seems to be a significant prognostic indicator in patients with colorectal carcinoma and the improved method described in this study would be useful for routine prognostication.
During epithelial carcinogenesis, it is not only the epithelial cells that change, but the stroma also. The morphological signs of cancer associated stromal alterations are desmoplasia, angiogenesis, and inflammatory cell infiltration.1,2 Cancer cell invasion usually induces new extracellular matrix formation and the activation of fibroblast-like cells, known as the desmoplastic reaction.1–3 This desmoplastic reaction is first observed with the invasive phenotype, and is not seen in carcinoma in situ or benign tumours.2 Some invasive tumours trigger a severe desmoplastic response characterised by abundant stroma rich in collagen.3
“We have devised an objective and simple method to measure desmoplasia in colorectal carcinomas by computer assisted image
The biological relevance of the desmoplastic host reaction is not fully understood. Initially, it was suggested that this response could play a protective role and limit the process of tumour invasion.3 However, recent evidence suggests that molecular crosstalk between neoplastic cells and stromal cells, and cancer induced changes in the stroma, modify the differentiation, proliferative capacity, and invasive capacity of tumour cells.1,4–6 In a xenograft model, injection of human prostate cancer cells combined with stromal components into athymic nude mice resulted in a high tumour incidence as a result of the stimulation of early angiogenesis and tumour growth, whereas prostate cancer cells alone were essentially non-tumorigenic.5 In another experimental study, the transplantation of colon carcinoma cells into reactive subcutaneous granulation tissue induced invasive colon carcinomas, but transplantation of colon carcinoma cells into undisturbed subcutaneous connective tissue resulted in the formation of encapsulated, non-invasive, well differentiated tumours.6
A variety of clinical studies of human lung and breast carcinomas and squamous cell carcinomas of the skin confirmed these experimental data, and suggested an association between desmoplasia and tumour invasion, development of metastases, or recurrence.7–9 In colon cancer, studies of desmoplasia gave more conflicting results. A small number of studies suggested that a peritumorous stromal response was an adverse prognostic factor.10,11 However, other studies of desmoplasia have failed to demonstrate prognostic significance for this factor.10 These discrepancies seem to be mostly the result of technical errors, such as subjective grading systems (including semiquantitative methods), interobserver and intraobserver variation, or sampling error. The case volume is also an important factor. The amount of desmoplasia in colorectal carcinoma was traditionally graded as: little, moderate, and extensive.12 Such a semiquantitative grading system is easy to use but seems to be unreliable. The College of American Pathologists recommended that desmoplasia should be evaluated by a uniform method as an individual variable in large studies on prognostic factors using multivariate analysis.12
To date, no quantitative measures have been used to assess desmoplasia in colorectal carcinomas, perhaps because of its unknown biological relevance. Hence, we have devised an objective and simple method to measure desmoplasia in colorectal carcinomas by computer assisted image analysis. The purpose of our study was to evaluate the prognostic relevance of desmoplasia in primary colorectal carcinomas using this novel method.
MATERIALS AND METHODS
Tissue sections were cut from the archival, paraffin wax embedded resection specimens of 112 patients with primary colorectal adenocarcinoma treated between 1991 and 1997. The clinical records of all patients were reviewed for age, sex, primary tumour site, tumour size, and the occurrence of distant metastasis. Haematoxylin and eosin stained sections of the specimens were reviewed by one investigator (BS) and conventional histopathological prognostic parameters in colorectal carcinoma, such as tumour grade, depth of invasion, lymph node involvement, and vascular and perineural invasion, were noted. Conventional semiquantitative grading for desmoplasia12 was also carried out. The tumours were staged according to the modified Astler-Coller staging system,13 and graded according to the World Health Organisation criteria.14 The clinical stage of disease was determined according to the AJCC staging system.15 Table 1 details the distribution of clinicopathological variables in the patients. All patients were followed up until death or for a median of 66 months (range, 48–120). During follow up, 38 (33.9%) patients died—37 as a result of distant metastasis and one as a result of postoperative complications. The median survival time was 42.5 months (range, 1–108). Twelve (10.7%) patients were out of follow up.
The most representative paraffin wax block from each tumour was chosen and cut into 5 μm thick sections. The sections were dewaxed, rehydrated, and then stained with van Gieson’s solution (saturated aqueous picric acid and 1% aqueous acid fuchsin) to demonstrate collagen rich stroma.16
For each case, the section stained by van Gieson’s method was observed by light microscopy and scored semiquantitatively for desmoplasia as: little, moderate, and extensive.
Digital images obtained from the van Gieson stained tumour sections using a light microscope (Labophot-2; Nikon, Tokyo, Japan) at ×4 magnification and a digital colour video camera (WK-C220E; Hitachi, Tokyo, Japan) with image analysis software (Mediscope; Dokuz Eylul University Clinic Engineering Department, Izmir, Turkey) were collected and stored on a computer (Vestel, Manisa, Turkey). In each case, 10 digital images (total area, 26.1 mm2) were selected and analysed by automatic processing.
The software was intended for colour segmentation including a determined neighbourhood of colours. The colour selected by the pathologist is accepted as the index for automatic choosing of the stained area.
For each captured area, the percentage of red stained area was determined semiautomatically. First, the red stained area was selected by the pathologist using visual observation by multiple clicking on the image. Subsequently, the system selected the areas with the same degree and configuration of staining, converted into pixel density units, and presented the area of selection as a percentage within each field. To allow the pathologist to repeat the procedure until the optimum selection was achieved, each selected view could be examined by the pathologist. The mean of 10 fields was determined and called the percentage of desmoplastic stroma (fig 1A, B), and the percentage of desmoplasia (PD) was determined by dividing the selected area by the total test area.17,18
One representative van Gieson stained tumour section of 13 resection specimens was measured separately by two observers, and seven specimens were measured twice at intervals by one observer. Next, we determined the interobserver variation, intraobserver variation, and coefficient of variation (CV (SD/mean)).19–21 To test the reproducibility, we also used scatter plots and calculated Spearman’s rank correlation coefficient (r). In addition, we also performed a reliability analysis by calculating the coefficient of intraclass correlation using ANOVA.19
Statistical analysis was performed using computer software (SPSS 12.0; SPSS, Chicago, Illinois, USA). A probability level of 0.05 or less was chosen to represent significance. Differences between groups were analysed with the Pearson χ2 test, one way ANOVA, and independent samples t tests. The relation between variables was investigated using Spearman’s rank correlation and Kendall’s tau-B correlation tests. Univariate survival analyses were performed with the Kaplan–Meier method and multivariate analyses with the Cox proportional hazards model. For Kaplan–Meier, comparisons among groups were carried out using the log rank test.
Reproducibility analysis of the measurement of desmoplasia
Table 2 shows the mean PD variance and CV for interobserver and intraobserver reproducibility studies. The CV of PD was 0.73 for intraobserver variation and 1.1 for interobserver variation. Reproducibility is equal to 100 − CV.21 As shown in table 2, the coefficient of intraclass correlation of 0.98 of both intra-assay and interassay studies showed that the measurement of desmoplasia by image analysis was a highly reliable method, and similar results of reliability analysis for computer assisted measurement of fibrosis have been obtained before.19 Spearman’s rank correlation coefficients also showed a substantial intraobserver (r = 0.97) and interobserver (r = 0.98) agreement (figs 2 and 3).
Percentage of desmoplasia
The mean PD was 4.85 (SD, 3.37; range, 1–18). The cases were also grouped according to the median value (4.0): PD ⩽ 4 (group 1) and > 4 (group 2).
Table 1 showed the relation between PD and the clinicopathological parameters. The mean (SD) PD of patients with stage I disease was 4.09 (2.34) and for stage IV patients this increased to 6.0 (4.52); however, the difference was not significant (p > 0.05; table 1; fig 4). PD was found to be significantly associated with lymph vessel and venous invasion (p = 0.006 and p = 0.02, respectively; figs 5 and 6). No significant associations were found between PD and tumour grade, perineural invasion, lymph node involvement, or distant metastasis (p > 0.05).
Using the semiquantitative method, 11 (10%) cases had little desmoplasia, 88 (78%) moderate, and 13 (11%) extensive. No significant relations were found between semiquantitative desmoplasia values and other clinicopathological parameters (p > 0.05). Table 3 provides a comparison of the desmoplasia values obtained by quantitative image analysis and the semiquantitative method. PD values were not significantly different in cases with different desmoplasia grades (p > 0.05). The median PD was similar among the different grades (fig 7).
In univariate survival analysis, tumour stage (p ⩽ 0.001), distant metastasis (p ⩽ 0.001), venous invasion (p = 0.006), perineural invasion (p = 0.001), semiquantitative desmoplasia values (p = 0.03), and PD (p = 0.008) were found to be significantly related to survival (table 4; fig 8). Multivariate analysis showed tumour stage (p ⩽ 0.001), distant metastasis (p ⩽ 0.001), and PD (p = 0.01) to be independent determinants of survival.
The biological relevance of the desmoplastic response to cancer is still unclear. Initially, it was suggested to act as a barrier and inhibit tumour invasion.3 Inhibition of the desmoplastic response in mice has been reported to increase the invasion and metastasis of the BL6 melanoma cell line.22 Hewitt and colleagues3 found little evidence of a desmoplastic response at the invasive edge of colorectal carcinomas—there was limited collagen I immunostaining and sparsely distributed or absent myofibroblasts—and suggested that desmoplasia fails to prevent cancer invasion because it is lacking at the invasive edge. It has also been proposed that myofibroblasts around in situ ductal breast carcinoma prevent cancer invasion by stabilisation of the peritumorous matrix through lysyl oxidase secretion.23
In contrast, there is growing evidence that the stromal reaction may play an important role in triggering cancer invasion.24,25 As reviewed by Ohtani,2 cancer invasion usually induces new extracellular matrix formation by the activation of stromal cells, which is a desmoplastic reaction, and this process includes both matrix production and degradation. It is already known that cancer stroma resembles granulation tissue because of the presence of myofibroblasts and active matrix degrading processes.1,2,6 Degradation of extracellular matrix plays a crucial role in tumour invasion and metastasis, and matrix degrading proteases, such as matrix metalloproteinases, are synthesised mainly by the stromal cells of cancer tissue.2,26
An interesting experimental study by Dingemans and colleagues6 showed that the transplantation of colon carcinoma cells into reactive subcutaneous granulation tissue induced invasive colon carcinomas, whereas transplantation into undisturbed subcutaneous connective tissue formed encapsulated, non-invasive, and well differentiated tumours. As reviewed by De Wever and Mareel1 and Bosman et al,27 the molecular crosstalk between cancer cells and the stroma may contribute to cancer invasion and may modulate differentiation in cancer cells. Martin et al proposed that stromal myofibroblasts alter the migratory properties of colon cancer cells and play a role in the invasion and metastasis of colorectal cancer.28 Similarly, Dimanche-Boitrel et al showed that stromal cells, more specifically myofibroblasts, contribute to the invasiveness of a rat colon cancer cell line in vivo and in vitro.29 In vitro experiments in rat colorectal tumour models by Lieubeau and colleagues30 showed that tumour activated myofibroblasts prevent the penetration of T cells and macrophages within tumour nodules, and they suggested that successful immunotherapeutic strategies should include complementary treatments against tumour associated myofibroblasts.
A few clinical studies point out an interesting association between intratumorous hypoxia and the presence of fibrosis in human malignant tumours, particularly in breast cancer.31,32 Intratumorous hypoxia significantly correlates with tumour invasion, metastasis, and poor prognosis in several malignant tumours, and it is also the driving force behind angiogenesis.31–35 Colpaert et al reported the presence of fibrotic foci in breast cancer as a marker of intratumorous hypoxia and hypoxia driven angiogenesis, and showed that it was associated with a worse survival.31,32
Despite the growing experimental evidence of the biological importance of tumour stroma, both in vitro and in vivo, a small number of clinical studies have shown the prognostic value of desmoplasia in colorectal carcinoma.10 These conflicting data are probably attributable to the small number of cases, the interobserver and intraobserver variation in semiquantitative scoring systems, and/or sampling variation.
In our present study, quantitative histochemical analysis of collagenous stroma (desmoplasia) was performed using van Gieson’s staining method, which stains collagen fibres red. Other collagen binding stains, such as picro-Sirius red or Masson’s trichrome stain, could be used with this computer assisted quantification method. We chose the van Gieson method because it is commonly used in most routine pathology laboratories, it is very easy to apply, and it does not need polarising microscopy.
“Only the measurement of desmoplasia by the quantitative method retained a clear independent association with overall survival”
Different morphometric methods have been described for the measurement of fibrosis, particularly in kidney and liver disease.19–21,36 In our current study, we measured the amount of desmoplasia in 112 primary colorectal carcinomas using computer assisted image analysis. The method was easy to perform, although it is more laborious and time consuming than semiquantitative techniques, and requires more complex equipment, including a digital capture system, computer, and software. However, it is the preferred method because it is reliable and allows desmoplasia to be measured more accurately. Using image analysis, the amount of red colour, representing the desmoplastic tumour stroma, was calculated and expressed as the PD. As shown in the cumulative survival curve, patients with a PD > 4 had a shorter survival than those with a PD ⩽ 4. Our results showed that the current semiquantitative scoring for desmoplasia is not a reliable method, because quantitative desmoplasia values were not significantly different among the different semiquantitative desmoplasia scores. Although univariate survival analysis showed both semiquantitative desmoplasia grades and PD to be relevant to prognosis, only the measurement of desmoplasia by the quantitative method retained a clear independent association with overall survival.
With a better understanding of the contribution of stromal tissue to cancer progression, new therapeutic agents might be targeted at inhibiting the crosstalk between cancer cells and the stroma.1,24 A potential target in anti-stroma treatment is transforming growth factor β, a well described cancer cell derived cytokine that activates the stroma and directly transdifferentiates fibroblasts into myofibroblasts.1 Other examples include the induction of apoptosis in myofibroblasts, antiangiogenic factors, and matrix metalloproteinase inhibitors.1,24,30 In addition, the extracellular matrix may contribute to the drug resistance of a solid tumour by preventing the penetration of therapeutic agents or reducing chemotherapy induced apoptosis in cancer cells.37,38 Netti et al detected an extended collagen network in the more penetration resistant tumours, and identified collagen as a potential target of treatment to improve drug penetration.37
Take home messages
The percentage of desmoplasia (PD) measured by image analysis was found to be significantly associated with lymph vessel and venous invasion and was associated with survival—patients with PD > 4 had a shorter survival than those with PD ⩽ 4
Thus, PD measured by image analysis seems to be a significant prognostic indicator in patients with colorectal carcinoma and is a more accurate predictor than semiquantitative methods
The improved method described in this study would be useful for routine prognostication
In our present study, the measurement of desmoplasia by computer assisted image analysis indicated that the amount of desmoplasia is an independent prognostic marker in primary colorectal carcinoma. With advances in our understanding of the contribution of the stroma to cancer progression, pathologists will probably be asked to assess the extent of desmoplasia in surgical pathology reports, for both prognostic and therapeutic purposes. Reliable, reproducible, and uncomplicated/easy to use methods for the measurement of desmoplasia, as proposed in our present study, should be established for further clinical trials.