A simple mathematical model utilizing topological invariants for automatic detection of tumor areas in digital tissue images
Diagnostic Pathology volume 8, Article number: S27 (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 ().
We treat two dimensional pathological images, the Betti numbers are consisting of two numbers. One is b0 (the 0-dimensional Betti number), which is the number of such isolated solid component as each cell or cell nucleus. The other is b1 (the 1-dimensional Betti number), which is the number of windows in the fenestrated area. That area is created by incomplete fusion of neighboring isolated solid component. The schematic illustration is shown in Figure 1.
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.
In Figure 2, we show a digital image of tissue of colon and its binarized image. In the tumor area, we can see many windows.
Results and discussion
The Betti numbers are obtained per unit binarized test image areas of normal and tumor tissues of colon. The values, b1 and b1/b0, of representative each five unit areas of normal and tumor colonic tissues are listed in Table 1.
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.
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.
Arkowitz M: “Introduction to Homotopy Theory (Universitext)”. 2011, Springer Press
CHomP (Computational Homology Project). [http://chomp.rutgers.edu/]
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.
The authors declare that they have no competing interests.
KN carried out most of experiments, participated in the design of the study and drafted the manuscript. TY and MN participated in the design of the study and helped the manuscript. All authors have read and approved the final manuscript.
About this article
Cite this article
Nakane, K., Tsuchihashi, Y. & Matsuura, N. A simple mathematical model utilizing topological invariants for automatic detection of tumor areas in digital tissue images. Diagn Pathol 8 (Suppl 1), S27 (2013). https://doi.org/10.1186/1746-1596-8-S1-S27
- Pathological Diagnosis
- Binarized Image
- Simplical Complex
- Colonic Tumor
- Solid Component