Volume 8 Supplement 1

Proceedings of the 11th European Congress on Telepathology and 5th International Congress on Virtual Microscopy

Open Access

Impact of tumor heterogeneity on disease-free survival in a series of 368 patients treated for a breast cancer

  • Myriam Oger1Email author,
  • Mohamed Allaoui2,
  • Nicolas Elie3,
  • Jacques Marnay2,
  • Paulette Herlin1,
  • Benoît Plancoulaine1,
  • Jacques Chasle1,
  • Véronique Becette1 and
  • Catherine Bor-Angelier1, 2
Diagnostic Pathology20138(Suppl 1):S43

DOI: 10.1186/1746-1596-8-S1-S43

Published: 30 September 2013


Tumor heterogeneity [14] is an old concept but its impact on the cancerogenesis process is poorly understood. Breast cancer is a noteworthy model for its frequency, and for the diversity of its phenotypes and of its evolution. This study examines the influence of the heterogeneity of tumor proliferation on disease-free survival of patients with a breast carcinoma.

Material and methods

Histological slides

The study involved a series of 368 patients from the François Baclesse Cancer Centre (Caen) treated for a breast carcinoma between 1991 and 1995, whitout neoadjuvant therapy and with a follow-up of more than 15 years. The table 1 contains the description of the series.
Table 1

Univariate Analysis of Disease Free Survival – 368 Eligible Patients. In grey: Follow up (2011)


No. of patients

P value



Mean = 58.7 yr






237 (65%)



124 (34%)




Right breast

176 (47.8%)


Left breast

187 (50.8%)


Synchronous bilateral

5 (01.4%)


Tumor size


Mean = 25 mm






212 (57.6%)



141 (38.3%)



15 (04.1%)


Excision quality



202 (57.2%)



47 (13.3%)



104 (29.5%)


Histological type



297 (80.71%)



37 (10.05%)


SBR grade


Grade 1

50 (14%)


Grade 2

180 (49%)



Grade 3

137 (37%)


Mitotic index (/1.7mm²)


Mean = 10 mitosis



Tumor vascular emboli



293 (79.6%)



75 (20.4%)


Lymph node metastasis



153 (44%)




195 (56%)


Hormone receptor status (at least 1)



265 (73%)




98 (27%)


Visceral or lymph node (other than axillary) metastasis



136 (37%)



232 (63%)


Local recurrence



66 (18%)



302 (82%)


Oncological event (Metastasis and/or local recurrence)



157 (42.5%)



211 (57.5%)





90 (24.46%)



278 (75.54%)

Histological sections, representative of each tumor, have been stained with the anti-phosphohistone-H3 antibody (PHH3: Ser10, MILLIPORE®, dilution 1/600) [5, 6]. With this specific immuno-stain, cells presenting mitosis figures are more easily identifiable (Figure 1).
Figure 1

dyes used to stain the histological sections. a: histological slide stained with HES; b: histological slide of the same case immunohistochemically stained with PHH3 [68] (anti-phosphohistone-H3 antibody); c: thanks to this specific immunohistochemical, identification of cells with mitosis figures, from prophase to telophase, is improved.


Histological slides have been scanned with a high resolution slide scanner to obtain virtual slides with a final resolution of 0.5 µm (ScanScope® CS from Aperio Technologies (20x NA 0.7 objective). The true color images obtained (color RGB 24 bits) have been saved in the tiled pyramidal TIFF file format.

Region of interest (ROI)

Before the automatic image analysis, the user can discard “normal” tissue surrounding the tumor by drawing a region of interest on the high resolution virtual slide with the Aperio ImageScope® software.

Image processing

The image processing was performed in two steps on a personal computer with a 1.6 GHz Pentium IV processor and a 1 GB of random access memory (RAM). The first step being a sub-sampling of virtual slide done with a specific algorithm 'Daubechies' second moment orthogonal wavelet decimation developed in C++ language which creates a low resolution image of the virtual slide (divided by 8: from 0.5µm to 4µm/pixels). In a second step, the low resolution image is automatically processed thanks to chaining operators of image analysis toolbox software (Aphelion, ADCIS).

In addition to estimating the frequency of mitotic figures, the program detects “hot spots” and measures 9 features representing the tumor heterogeneity, including the Haralick texture features and Fisher’s index. The zones of influence of each stained nuclei have been determined using Voronoï’s pavement principle. When nuclei are close, the size of pavements is small, highlighting the “hot spots”.

Feature selection

A principal component analysis has been done in order to select the most relevant features.

Statistic analysis

These features have been statistically analyzed, combined with classic clinic-pathological prognostic factors (age, tumor size, grading, mitotic index, vascular emboli and metastatic lymph nodes).


Principal component analysis

Thanks to the principal component analysis (PCA) 4 features representing tumor heterogeneity have been chosen then combined into three new features: CP1, CP2 and CP3, corresponding to the three principal directions of the PCA.

The four selected features are:
  • 2 Haralick’s texture indexes (correlation and energy);

  • Fisher’s index;

  • variance of the size of Voronoï pavements.

The variance of the size of Voronoï pavement (named Voronoï) and the Fisher’s index are regional features whereas the Haralick’s texture indexes are local features. Indeed, Voronoï and Fisher features are “cutting” the tissue into pieces and analyzing each of them compared to the others, whereas Haralick is dealing with relations between neighbor pixels, each pixel representing a cell at this resolution.

Prognostic study

In the analysis of prognostic factors, disease free survival was used as the end point.

Univariate statistical analysis (DFS)

Univariate analysis of disease free survival was performed with the features of age, tumor location, initial tumor size, pathologic lymph node status (N), histological type, SBR grade, mitotic index, vascular emboli, metastatic lymph nodes and hormone receptor status. The results are shown in Table 1 for usual features, in Table 2 for heterogeneity features.
Table 2

Results of the univariate analysis


P value



Normalized variance of density


Energy (Haralick)








The CP2 feature correlated highly with disease free survival, whereas the variance of the Voronoï pavements was borderline significant.

Multivariate statistical analysis (Cox)

The above features that correlated with disease free survival in univariate analysis were combined with clinic-pathologic factors and included in the multivariate analysis. Cox’s regression analysis highlighted 3 independent prognostic factors: tumor heterogeneity feature CP2 (RR = 1.46; p = 0.03), mitotic index (RR = 1.71; p = 0.004) and lymph node metastasis (RR = 2.20, p < 0.0001) correlated highly with disease free survival.

The construction of this model has individualized 3 groups of patients: 0 factor, 1 or 2 factors and 3 poor prognostic factors (mitotic index > 10, lymph node metastasis in the axillary dissection, upper tercile of CP2; p < 0.0001).

Disease free survival according to this model is shown in Figure 2.
Figure 2

Model build using the results of the multivariate analysis. Multivariate analysis highlighted 3 independent prognostic factors: tumor heterogeneity (CP2), mitotic index, lymph node metastasis. The construction of this model has individualized 3 groups of patients: 0 factor, 1 or 2 factors and 3 poor prognostic factors (mitotic index > 10, lymph node metastasis in the axillary dissection, upper tercile of CP2; p < 0.0001).

Discussion and conclusion

To characterize tumor heterogeneity in the presented series of breast cancer, 9 features were computed. 4 nonredundant of them have been selected by principal component analysis (PCA).

PCA was also used to create 3 new composite features: CP1, CP2 and CP3, corresponding to the 3 principal directions of the PCA.

The univariate analysis made for each feature from image analysis has first highlighted that only the combination CP2 and Voronoï’s feature had a prognostic value. It has to be noted that a high value of heterogeneity index is associated with a poor prognosis.

In multivariate analysis, CP2 was found to be an independent prognostic feature just like the mitotic index and the lymph node status. The lymph node status is a well-known clinical factor; the two other features are intrinsic factors of tumor growth, at cellular level for mitotic index and at the tissue level for heterogeneity.

Surprisingly, age, tumor size, Scarff and Bloom Grade and hormone receptor status are of secondary importance compared to these 3 features.

This result encourages to confront the heterogeneity feature CP2 to clinic information, such as recent or late oncologic event or the nature locoregional or distant visceral of the recurrence, and to the absence of lymph node metastasis.


Authors’ Affiliations

Imagin' Team of EA 4656 in François Baclesse Cancer Centre
Pathology laboratory, François Baclesse Cancer Centre
Plateau d’Histo-Imagerie Quantitative, SF ICORE, CMABIO, University of Caen Basse-Normandie


  1. Dodd LG, Kerns B-J, Dodge RK, Layfield LJ: Intratumoral Heterogeneity in Primary Breast Carcinoma: Study of Concurrent Parameters. Journal of Surgical Oncology. 1997, 64: 280-288. 10.1002/(SICI)1096-9098(199704)64:4<280::AID-JSO6>3.0.CO;2-5.View ArticlePubMedGoogle Scholar
  2. Sharifi-Salamatian V, de Roquancourt A, Rigaut JP: Breast carcinoma, intratumour heterogeneity and histological grading, using geostatistics. Analytical Cellular Pathology. 2000, IOS Press, 20: 83-91. 0921-8912Google Scholar
  3. Bertucci F, Birnbaum D: Reasons for breast cancer heterogeneity. Journal of Biology. 2008, 7: 6-10.1186/jbiol67.PubMed CentralView ArticlePubMedGoogle Scholar
  4. Fisher B, Redmond CK: Evolution of Knowledge Related to Breast Cancer Heterogeneity: A 25-Year Retrospective. Journal of Clinical Oncology. 2008, 26 (13): 2068-2071. 10.1200/JCO.2007.14.1804.View ArticlePubMedGoogle Scholar
  5. Hans F, Dimitrov S: Histone H3 phosphorylation and cell division. Oncogene. 2001, 20: 3021-3027. 10.1038/sj.onc.1204326.View ArticlePubMedGoogle Scholar
  6. Bossard C: Phosphohistone H3 labelling for histoprognostic grading of breast adenocarcinomas and computer-assisted determination of mitotic index. J Clin Pathol. 2006, 59 (7): 706-10. 10.1136/jcp.2005.030452.PubMed CentralView ArticlePubMedGoogle Scholar
  7. Gosme M: Modélisation du développement spatio-temporel des maladies d’origine tellurique, thèse de doctorat. 2007, Rennes: Université de Rennes I, 198-Google Scholar
  8. Fisher RA: The use of multiple measurements in taxonomic problems. Ann. Eugenics. 1936, 7: 179-188. 10.1111/j.1469-1809.1936.tb02137.x.View ArticleGoogle Scholar
  9. Marcelpoil R, Usson Y: Methods for the study of cellular sociology: Voronoi diagrams and parametrization of the spatial relationships. J. Theor. Biol. 1992, 154: 359-369. 10.1016/S0022-5193(05)80176-6.View ArticleGoogle Scholar
  10. Haralick RM, Shanmugam K, Dinstein I: Textural features for image classification. IEEE Transaction on Systems, Man and Cybernetics. 1973, 3 (6): 610-621.View ArticleGoogle Scholar


© Oger et al; licensee BioMed Central Ltd. 2013

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.