A simple mathematical model utilizing topological invariants for automatic detection of tumor areas in digital tissue images
© Nakane et al; licensee BioMed Central Ltd. 2013
Published: 30 September 2013
Pathological diagnosis starts usually from detection of morphological deviations of tissues and cells in question from their normal counterparts. Measurement and analysis of how they are deviated are the processes of actual pathological diagnosis and many of them can be achieved by computer using certain mathematical models. Computer assisted pathological diagnosis is now an issue of importance in the current situations of shortage of diagnostic pathologists in Japan. In the present study we propose a simple mathematical model to differentiate tumor tissues from their normal counterparts utilizing changes in the Betti numbers in tumorigenesis. We tested our method in several normal and tumor tissue images and preliminary results are reported.
Materials and methods
The Betti numbers are the numbers coined after Enrico Betti, an Italian mathematician of topology. The Betti numbers are one of the invariants in a homology. The invariant implies the quantity that is unchangeable by continuous transformation ().
Morphological changes during the course of tumorigenesis
Normal cells are usually controlled so as to grow and be reproduced only when a tissue composed of the cells requires new cells. That is, a new cell is generated when the old cell dies or is damaged. However, if the foregoing process of cell growth and cell division becomes out of order, the cells grow and divide, excessively. The excessive cells then produce a mass of tissues called “a tumor” or “a neoplasm”. The disordered cell growth causes a variety of morphological characteristics of tumor tissues.
Mathematical representation of morphological change in tumorigenesis
Uncontrolled cell growth usually makes cell nuclei become larger. This causes a change in contact condition between components in a tissue. This change appears along with neoplastic changes and tumorigenesis. The following description deals with how to numerically express the changes.
As shown in Figure 1, when the number of contact points between individual solid components (cells or cell nuclei) increases, i.e., 7 solid components (b0) fuse each other through contact, 2 solid components with fenestration are created. They have 4 windows (b1). Namely, b1 increase and the ratio b1/b0 significantly change.
Numerical procedure of our method
Hematoxylin and eosin stained mucosal biopsy sections of colonic tumors are used as test samples. First they are taken as whole slide images and then binarized setting a threshold value depending on each HE stained condition. We regard the binarized image as a sum of a simplical complex. By using a free calculation software (cf. ), we can calculate the Betti numbers. We compare the values, b1 and b1/b0, of colonic tumor tissues with those of normal tissues.
Results and discussion
b1 and b1/b0 of the normal and tumor tissue.
In normal tissue, b1 is less than 1400 (565.6 average), b1/b0 is less than 1.4 (0.251 average). If tissue is malignant, b1 is large (3266.4 average), b1/b0 is increased (3.79 average).
Based on the results obtained above, the indices for a normal tissue and a tumor tissue can be assumed.
Marks are put on each unit area as the value of index
None (~ 1400)
▽ (1400 ~ 2000)
× (2000 ~ )
△ (1.4 ~ 1.7)
○ (1.7 ~ )
The Betti numbers can express the degree of the connection of the figure. When the number of contact points between individual components increases, irrespective of their shape, b1 and a ratio b1/b0 significantly change. So we hypothesize that these numbers can be used as indexes to represent the cellular “accumulation” which is one of the characteristics of tumor tissue.
The numerical results indicate that the difference in our indices can differentiate tumor tissues from their normal counterparts.
We would like to take this opportunity to thank Dr. Nagumo (research professor of graduate school of medicine, Osaka University) for her valuable advice from a pathological standpoint. Further, we would like to thank Prof. Suematsu J. Nobuhiko (graduate school of advanced mathematical science Meiji University) for his valuable advice on how to display calculated results on an image.
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.