ORIGINAL RESEARCH article
The Quantitative Criteria Based on the Fractal Dimensions, Entropy, and Lacunarity for the Spatial Distribution of Cancer Cell Nuclei Enable Identification of Low or High Aggressive Prostate Carcinomas
- 1Department of Urology, Alb Fils Kliniken, Goeppingen, Germany
- 2The Bȩdlewo Institute for Complexity Research, Poznań, Poland
Background: Tumor grading, PSA concentration, and stage determine a risk of prostate cancer patients with accuracy of about 70%. An approach based on the fractal geometrical model was proposed to eliminate subjectivity from the evaluation of tumor aggressiveness and to improve the prediction. This study was undertaken to validate classes of equivalence for the spatial distribution of cancer cell nuclei in a larger, independent set of prostate carcinomas.
Methods: The global fractal capacity D0, information D1 and correlation D2 dimension, the local fractal dimension (LFD) and the local connected fractal dimension (LCFD), Shannon entropy H and lacunarity λ were measured using computer algorithms in digitalized images of both the reference set (n = 60) and the test set (n = 208) of prostate carcinomas.
Results: Prostate carcinomas were re-stratified into seven classes of equivalence. The cut-off D0-values 1.5450, 1.5820, 1.6270, 1.6490, 1.6980, 1.7640 defined the classes from C1 to C7, respectively. The other measures but the D1 failed to define the same classes of equivalence. The pairs (D0, LFD), (D0, H), (D0, λ), (D1, LFD), (D1, H), (D1, λ) characterized the spatial distribution of cancer cell nuclei in each class. The co-application of those measures enabled the subordination of prostate carcinomas to one out of three clusters associated with different tumor aggressiveness. For D0 < 1.5820, LFD < 1.3, LCFD > 1.5, H < 0.7, and λ > 0.8, the class C1 or C2 contains low complexity low aggressive carcinomas exclusively. For D0 > 1.6980, LFD > 1.7644, LCFD > 1.7051, H > 0.9, and λ < 0.7, the class C6 or C7 contains high complexity high aggressive carcinomas.
Conclusions: The cut-off D0-values defining the classes of equivalence were validated in this study. The cluster analysis suggested that the number of the subjective Gleason grades and the number of the objective classes of equivalence could be decreased from seven to three without a loss of clinically relevant information. Two novel quantitative criteria based on the complexity and the diversity measures enabled the identification of low or high aggressive prostate carcinomas and should be verified in the future multicenter, randomized studies.
The spatial distribution of cancer cell nuclei in digitalized histological images of prostate carcinomas reflects complex intercellular interactions within malignant tumor underlying its biological aggressiveness. Since this distribution can be characterized by a number of complexity and diversity measures, those measures define the quantitative criteria for low complexity carcinomas with low aggressiveness or for high complexity carcinomas with high aggressiveness.
Histological evaluation of tumor aggressiveness, that is a potential of cancer cells for local growth and metastasis formation, is a key element in both diagnosis and prognosis of any malignancy. In the case of prostate carcinomas, pathologists developed about 40 subjective grading systems scaling aggressiveness of those tumors (reviewed in Humphrey, 2003). The Gleason system (Gleason, 1977) is the most commonly used one (Humphrey, 2004). It was validated with a database of a few thousand prostate cancer patients. Since its revision in 2005 and in 2010, the Gleason system comprises three patterns of tumor growth called the Gleason grades (Epstein et al., 2005; Epstein, 2010). The Gleason score is calculated as a sum of two or, sometimes, three Gleason grades representing most predominant patterns of growth in a given carcinoma. This information co-determines the risk stratification of cancer patients. According to that stratification, cancer patients will undergo follow-up or will be treated. That latter option comprises a choice of the optimal therapy, extent of surgical resection with or without nerve sparing or dosing of radiation therapy (Epstein et al., 2005; Epstein, 2010; Lotan and Epstein, 2010).
A prostate biopsy is a procedure to obtain small samples of prostate tissue for the microscopic examination using a thin needle. Owing to both the limitations in obtaining a representative tissue specimen and a low specificity of this procedure, there is no good correlation between the grading score of the prostate carcinomas in biopsy tissues and in prostatectomy specimen (Fine et al., 2012). Using a prostate biopsy, one cannot predict the existence of insignificant (indolent) cancer, that is cancer without a significant influence on survival (Epstein et al., 1994), the final histological grade, size of tumor, its extracapsular extension, the existence of positive margins or lymph node involvement (D'Amico et al., 2000; Boccon-Gibod et al., 2005; Fine et al., 2012). In addition, the modification of the criteria in the Gleason system (Epstein et al., 2005; Epstein, 2010) caused both a discrepancy in tumor grading and difficulties in comparison of treatment results for patients treated before and after the change (Lotan and Epstein, 2010). Although tumor aggressiveness is a key parameter in the prediction models, the subjectivity of histological evaluation is a source of the prediction weakness. Human eye is not able to evaluate correctly many details in tumor images, such as the ratio of cells bound in glands to the number of infiltrating cancer cells, or some configurations, such as the co-existence of glands of different size and shape. It is a challenge for pathologists to match the images to the definition of a grade, especially in the borderline cases (reviewed in Montironi et al., 2013). The subjective tumor grading has a significant inter- and intraobserver variability in the range of 38–80% (Nguyen et al., 2004; van der Kwast et al., 2010; McKenney et al., 2011; Netto et al., 2011; Egevad et al., 2013; Berney et al., 2014) and the coefficient κ for interobserver agreement 0.15–0.7 (McLean et al., 1997; Allsbrook et al., 2001; van der Kwast et al., 2010; McKenney et al., 2011; Scott Lucia et al., 2013). Even though some DNA-, RNA-, or protein biomarkers were co-applied, it did not improve the accuracy of grading. The value of the correlation coefficient was low in the range of 0.1–0.6 (McDunn et al., 2013; Pin et al., 2013). A novel diagnostic tool, prostate magnetic resonance imaging reinforced by the multiwavelet filters was supposed to enable the automated diagnosis of prostate carcinoma. The results remain nonspecific. Biopsies must be performed to evaluate tumor aggressiveness nonetheless (reviewed in Turkbey and Choyke, 2012).
A decision whether a prostate cancer patient should be treated or not depends on the fulfillment of two criteria [Klein, 2014; Mottet et al., 2015; National Comprehensive Cancer Network (NCCN), 2015]. The first one is the patient's overall life expectancy determined by age, health status, and comorbidities. In particular, patients older than 75 years, with life expectancy less than 10 years, in a reduced health status or with serious comorbidities will not be treated with the intention to cure. This category of patients will be offered a program of watchful waiting with a follow-up schedule. A palliative treatment can be initiated if local or metastatic progression with clinical symptoms occurs. The second criterion plays a pivotal role in the risk assessment of cancer progression. It considers three parameters (Partin et al., 1997), such as biological aggressiveness of tumor described by the subjective tumor grading according to Gleason (Gleason, 1977; Humphrey, 2003, 2004; Epstein et al., 2005; Epstein, 2010), PSA concentration in serum (Stamey et al., 1989; D'Amico et al., 2004, 2005; Thompson et al., 2006; ElShafei et al., 2013; Wu et al., 2014), and tumor stage (Prout et al., 1980; Petros and Catalona, 1992). The statistical accuracy of the nomograms based on those parameters approximates 70% only (Partin et al., 1997; Miller et al., 2006; Thompson et al., 2006; Klotz et al., 2010).
According to the clinical guidelines [Klein, 2014; Mottet et al., 2015; National Comprehensive Cancer Network (NCCN), 2015], low risk cancer patients with PSA < 10 ng/ml, Gleason score of 6 or less, and pT1-2a stage do not need any treatment, nor an adjuvant therapy. Despite of those criteria, the overtreatment occurs at the indolent stadium of the disease in about 55% of patients (Miller et al., 2006; Danzig and McKiernan, 2014; Daskivich et al., 2014). Urologists recommend frequently treatment even for 80-year-old patients with limited life expectancy (Dall'Era et al., 2008; Albertsen, 2010; Danzig and McKiernan, 2014; Daskivich et al., 2014). Rather, those patients should be followed-up to identify local tumor progression and aggressive metastatic disease in the long term (Klotz, 2006; Denis, 2007; Dall'Era et al., 2008; Albertsen, 2010); phenomena especially frequent about 15 years from the diagnosis (Albertsen et al., 2005; Cooperberg et al., 2011; Cooperberg, 2015; Klotz et al., 2014). The same strategy might also be offered to some patients with intermediate risk of progression (PSA 10–20 ng/ml, Gleason score 7, pT2b-2c) (Cooperberg et al., 2011). Unfortunately, only up to 49% of eligible patients participate in the follow-up programs (Womble et al., 2015). The lack of the histological, molecular or clinical criteria for the unequivocal identification of the low risk cancer patients increases a fear for tumor progression among those involved in active surveillance as well as among those who already underwent some kind of aggressive therapy, without even knowing if they needed it (Womble et al., 2015).
Patients with high-risk (PSA >20 ng/ml, Gleason score 8–10, stage T3a) or very high risk carcinomas (stage T3b–T4 N0 or any T N1) should be treated aggressively with all available modalities, such as surgery, radiotherapy, brachytherapy, or hormonal therapy, even though a choice of a modality seems to increase survival in about 10% of patients (Heidenreich et al., 2014). For example, SEER, a population-based cohort study of 66,717 patients 66 years of age or older with localized prostate cancer demonstrated that primary hormonal therapy did not improve long-term overall or disease-specific survival (Lu-Yao et al., 2014).
The relevant question, whether a given prostate carcinoma is an indolent, low grade, low risk malignancy rather than a high grade, high risk one can be answered only if prostate gland is removed surgically. Although this answer appears after the initiation of the aggressive treatment, it is still important to get the unequivocal evaluation of tumor aggressiveness based on some objective, quantitative criteria. Those criteria might help to improve the accuracy of prognosis, to plan more individualized therapy, or to facilitate a search for some features of cancer cells that would improve predictions made on the basis of biopsies. A variety of approaches, such as graph models (Altunbay et al., 2010), gland segmentation with morphometric analysis and application of probabilistic pairwise Markov models (Naik et al., 2008; Loeffler et al., 2012; Nguyen et al., 2012), color channel histograms, texture analysis, fractal algorithm with fractal code or multifeature analysis with a number of classifiers, such as Bayesian one (Doyle et al., 2007), multiwavelet method (Jafari-Khouzani and Soltanian-Zadeh, 2003; Candes et al., 2006; Huang and Lee, 2009; Hong-Jun et al., 2011; Veltri et al., 2012; Lopez and Again, 2013), k-NN one (Tabesh et al., 2005), support vector machine (Tabesh et al., 2007), or linear discriminant analysis method (Tabesh et al., 2005, 2007) were proposed for the automated diagnosis of prostate carcinoma or for the evaluation of tumor aggressiveness. Those statistically elegant and mathematically advanced approaches resulted in the maximal accuracy of 78–81% for grading in the case of multiwavelet method (Tabesh et al., 2005, 2007). The accuracy of 95–97% was reported only for the low number of carcinomas (Lopez and Again, 2013). This discrepancy between the Gleason grading and the quantitative classifications generated by the sophisticated computer algorithms had at least one reason. Two different categories of data, the subjective and the objective one, should not be compared each other if the subjective tumor grading according to Gleason is chosen as the absolute frame of reference.
Since, there is no golden standard in pathology of adenocarcinomas that might be used as the objective frame of reference for the quantitative studies, a novel approach based on the fractal geometric model of prostate carcinomas and evaluation of complexity in the spatial distribution of cancer cell nuclei was proposed (Waliszewski, 2013; Waliszewski et al., 2014, 2015). Cell nuclei are fragments of the surface in two-dimensional histological images. They compose irregular, yet self-affine configurations (Waliszewski, 2013). The fractal geometrical model of prostate carcinomas introduced the idea of the distinct classes of equivalence called complexity classes (Waliszewski et al., 2015). Those classes were defined by the values of the global capacity fractal dimension D0. Using that model, all basic patterns of growth seen in prostate cancer were subordinated to the appropriate class. In particular, low grade prostate carcinomas with well-preserved glandular structure were re-stratified to the class C1, C2, or C3. High grade carcinomas with the predominating infiltrating patterns of growth were stratified to the class C6 or C7.
This study was undertaken to validate the classes of equivalence (Waliszewski et al., 2015) using utterly novel set of carcinomas than that used in the study (Waliszewski, 2013; Waliszewski et al., 2014, 2015). A question arises whether the additional complexity or diversity measures, such as the local fractal dimensions, Shannon entropy or lacunarity allow a similar stratification of prostate carcinomas into the same classes of equivalence? Is it possible to define some quantitative criteria based on all those parameters that define prostate carcinomas with low or high complexity of the spatial distribution of cancer cell nuclei associated with low or high tumor aggressiveness?
The Renyi Family of the Global Fractal Dimensions
The Renyi family of the global spatial fractal dimensions comprises the capacity (D0), information (D1), and correlation (D2) dimension. The appropriate introduction to the theory of dimension, fractals, or the formal mathematical definitions of the fractal dimensions can be found elsewhere (Engelking, 1978; Czyz, 1994). Briefly, the capacity dimension D0 is defined as a relationship between the logarithm of a number of boxes N(ε) covering the geometric object and the logarithm of a box size ε (Equation 1; Vicsek, 1992; Kinsner, 2007):
The information dimension D1 measures how the average information needed to identify an occupied box, scales, as the scale of boxes gets smaller. The algorithm for the information dimension will search for the linear relationship between the logarithm of a box size (ε) and the logarithm of the probability p that a given box contains the element of the object (Equation 2; Baker and Gollub, 1996).
The correlation dimension D2 measures the number of points M used to generate a representation of the fractal and the number of pairs of points closer than ε to each other. The correlation dimension is a probability measure that two pixels within the object are close to each other less than ε. The appropriate algorithm calculates not only if a given box is occupied by the pixel of the analyzed geometric object, but also how many pixels in the box are, and how close to each other they are (Equation 3):
in which M is a number of pixels of the analyzed geometric object and gε denotes a number of pairs of pixels that are closer to each other than ε (Grassberger and Proccacia, 1983; Baker and Gollub, 1996).
Local Connected Fractal Dimension
The local fractal dimension (LFD) is the fractal capacity dimension calculated according to Equation 1 for every pixel in the image. This dimension is a local and indirect measure of intercellular connectivity, i.e., interconnectedness which denotes the existence of complex, dynamic relationships in a population of cells leading to the spatial and temporal emergence of global features in the system that would never appear in a single cell existing out of the system (Waliszewski, 1997, 2005, 2009; Waliszewski et al., 1998, 2001; Waliszewski and Konarski, 2001).
The local connected fractal dimension (LCFD) characterizes local irregularities of geometry of heterogeneous geometrical objects, such as clusters of cancer cells present either in glands or in infiltrates. Instead of a single value of the global fractal capacity dimension D0 calculated for the entire image, one gets a set of values calculated for each pixel that belongs to the analyzed object. This is done according to the following Image J algorithm:
Step 1: choose a pixel P that belongs to the analyzed object and possesses eight neighbor pixels.
Step 2: define the local connected set of pixels by finding all the pixels connected to the pixel P within the increasing s-pixel-side window centered at P.
Step 3: count how many pixels N(s) of the analyzed object are within the window.
Step 4: use the least square method to compute the slope of the log-log curve composed by the coordinates [log(N(s), log(s); Landini et al., 1995].
Entropy and the Global Information Fractal Dimension D1
The Shannon entropy H is a statistical measure of both the information content and diversity in a given cell configuration in the image. It can be used to characterize gray images or binary images of the spatial distribution of cancer cell nuclei. The Shannon entropy describes the following Equation (4):
In the above equation, pi denotes the probability that the difference between two adjacent pixels of the image is equal to i, and is calculated from the histogram counts. Low entropy images, such as those containing mostly black background, e.g., in benign prostatic hyperplasia or in well-differentiated adenocarcinomas, have large numbers of pixels with the same zero value. The black image that contains no white pixels will have entropy equal zero.
The entropy HD of a set with D dimensions is defined by the following Equation (5):
ín which H stands for the Shannon entropy (see Equation 4), ε denotes size. In consequence, entropy H is a linear function of log ε with slope D and intersection HD. It is important to notice that the slope D in the limit for ε = 0 becomes the global information fractal dimension D1 (see Equation 2; Sethna, 2011).
Cancer cells form various morphological patterns. The areas of the image complementary to the lumen of glands correspond to the gaps on the binary image. Investigating those gaps provides important information about homogeneity of cell distribution in malignant tumor. Lacunarity λ, that is, a degree of gappiness, inhomogeneity, or translational and rotational invariance in the image characterizes quantitatively those gaps as well as differences in their spatial distribution between the patterns of growth. This parameter can also quantify differences between self-similar structures having the same global fractal dimension (Plotnick et al., 1996).
In general, lacunarity λ is equivalent to the square of the coefficient of variation, CV, and expresses a relationship between standard deviation σ and mean μ for the pixels in a given area of an image, and a number of pairs of pixels g that are closer to each other than ε:
More specifically, lacunarity λ is a relationship between variation for the function P(m, r), that defines probability of the localization of the pixels m in the square of size r to the square of the mean value of that function (Equation 7):
The computer algorithms that calculate the mean value of lacunarity for images use the same principle as in the case of the global fractal dimensions, i.e., they analyze the digital image from different scaled levels of resolution to examine how certain geometric features change with the size of the element used to inspect the image. The first approach, called box counting, is identical with the principle underlying the capacity dimension (see above) in the sense that the box for each ε is placed as though it were part of a grid overlaid on the image so that the box does not overlap itself. In the second one, called the sliding box algorithm, the box is slid over the image so that it overlaps itself and the sliding box lacunarity is calculated. The values are usually expressed as logarithms. Lacunarity analyses using the types of values discussed above have shown that data sets extracted from geometric fractals, or from patterns that change little when rotated have low lacunarity. The larger and more irregular are the gaps in the image, as those seen in glands in benign prostatic hyperplasia, the greater is the value of lacunarity. In some instances, fractal dimensions and values of lacunarity can be correlated, but it does not hold for all types of patterns and measures of lacunarity.
Patients, Preparation of Images, and Isolation of Cell Nuclei
This study was performed according to the ethical standards outlined in the WMA Declaration of Helsinki “Ethical Principles for Medical Research Involving Human Subjects” (http://www.wma.net). The data bank was constructed and analyzed according to the ethical requirements of the Ethics Committee for the Medical Faculty of the Justus-Liebig University in Giessen, Germany (ethical vote 49/05). Since this project was based solely on computer-aided image analysis of the digitalized anonymous tissue slides of prostate carcinomas, no written consent from the prostate cancer patients was necessary. The surgical procedures were performed from 2007 to 2010. Tissues were stained with hematoxylin and eosin as well as with the appropriate antibodies against the markers, i.e., PSA (DAKO, Germany) and AMACR (DAKO, Germany) for prostate cancer cells, cytokeratin 34βE12 (Leica Biosystems, Canada) for basal cells in order to confirm the diagnosis (Kristiansen and Epstein, 2014), and classified by the Gleason criteria (Gleason, 1977; Humphrey, 2003, 2004; Epstein et al., 2005; Epstein, 2010).
A reference set of prostate carcinomas comprised cases representing the homogeneous patterns with the Gleason score 3+3 = 6 defined as a regular gland architecture with no basal cells and no areas of cellular infiltration (n = 20), the Gleason score 4+4 = 8 defined as dominating small glands spread freely in the prostate tissue or organized in the cords (n = 20), and the Gleason score 5+5 = 10 defined as a random infiltration of the prostate tissue by cancer cells without any glands (n = 20), and benign prostatic hyperplasia (n = 10). We analyzed and re-stratified a test set of prostate carcinomas comprising 208 carcinomas with the following Gleason score: 3+3 = 6 (n = 70), 3+4 = 7a (n = 18), 4+3 = 7b (n = 14), 4+4 = 8 (n = 23), 4+5 = 9 (n = 28), 5+4 = 9 (n = 20), 5+5 (n = 35). All prostate carcinomas were subordinated to the appropriate subsets called structural Gleason classes using the Gleason criteria. We also analyzed benign prostatic hyperplasia (n = 20). None of the patients had lymphatic or distant metastases.
Histological slides of prostate carcinomas stained with hematoxylin and eosin were digitalized using both a microscope Axioskop 5.0 with the halogen lamp 12V 50W 2800K 950 lm and the objective Plan Neofluar 20x with the numerical aperture 0.5, Zeiss, Germany and a camera Nikon Coolscope, Japan. The optimal conditions for the magnification (20x) and lighting intensity (41.19 × 106 lx) were chosen in such a manner that the values of the capacity fractal dimension for the test image were in the plateau-area of the test curve (data not shown).
The color images in the bmp format with 1240 × 1000 pixels and resolution 150 × 150 dpi were first resized to the jpg format with 648 × 432 pixels and resolution 150 × 150 dpi. Cancer cell nuclei were isolated electronically using a package Definiens Tissue Map ver. 7.0 (Definiens, Munich, Germany), and stored as RPG images in the jpg format with 648 × 432 pixels, resolution 150 × 150 dpi. Both conversion to the eight-bit images and their thresholding to the binary images, i.e., images having only 0 (black) and 255 (white) pixel values was performed using the Renyi Entropy filter of the open-source software Image J ver 1.48v (NIH, Bethesda, USA, http://imagej.nih.gov/ij/). The images were analyzed by the computer algorithms measuring the global fractal dimensions of the Renyi family, i.e. capacity dimension D0 and information dimension D1 (Benoit 1.3, True Soft, USA, http://www.trusoft-international.com/benoit.html), and correlation dimension D2 (Fractalyse 2.4, CNRS Universite de Franche-Comte, France, http://www.fractalyse.org). The values of the global fractal dimensions generated by the computer algorithms were verified against the geometrical model (Waliszewski, 2013; Waliszewski et al., 2014, 2015). Both the local fractal dimension and the local connected fractal dimension were measured in the same digitalized images using the open-source software Image J ver 1.48v and open-source plugin FracLac 2013 Janb420 by A. Karperien (http://imagej.nih.gov/ij/plugins).
The statistical analysis including the ROC analysis was performed by Sigma Plot ver. 10.0 (Systat Software Inc., San Jose, USA, http://www.sigmaplot.com). The cluster analysis was performed by Mathematica ver. 10.0 (Wolfram Research Inc., USA) and Statistica ver. 8.0 (Statsoft, Oklahoma City, USA). The following aspects are important while analyzing a ROC curve. The ROC curve is a two dimensional graph. It represents a relationship between sensitivity of a statistical test (y−axis) in a function of 1–specificity of that test (x−axis). Sensitivity is defined as a rate of true positive events. The second quantity represents a rate of non-events that were falsely identified as positive events for the different cut-off values of a binary classifier, such as D0, D1, or D2, LFD, etc. Second, the closer the ROC curve follows the y-axis, and then the top border of the ROC diagram, the more accurate the statistical test. The area under the ROC curve (AUC) is a probability of the test accuracy, i.e., a probability that the classifier discriminates between the elements of two sets correctly. In particular, the AUC of 1.0 denotes a perfect relationship between the data that belong to a given category and those that do not. Third, the closer the ROC curve comes to the diagonal of the ROC diagram, the AUC is closer 0.5, what denotes that the test is not able to discriminate between two sets of data. It should be noted that a task of the re-stratification of prostate carcinomas on the basis of the quantitative characteristics of the spatial distribution of cancer cell nuclei by the fractal dimensions D0, D1, and D2 into the complexity classes is equivalent to the task of the identification the cut-off values of the best classifier, i.e., the classifier generating the AUC of 1.0 (Hill and Lewicki, 2006).
Cluster analysis is a data mining tool that uses a number of computer algorithms and statistical methods to sort elements of a set into groups (clusters) in such a manner that the elements belonging to the same cluster have a maximal degree of association and a minimal degree otherwise. In this work, a v-fold cross-validation algorithm was applied to determine the number of clusters in the data. Then, the k-means algorithm was applied to produce exactly k different clusters predicted by the first algorithm with the greatest possible distinction and maximal distances between those clusters. One examines here the means for each cluster on each dimension (parameter) to assess how distinct the k-means clusters are. Ideally, the means for all dimensions should be different, and the magnitude of the F-values from the analysis of variance performed on each dimension indicates how well the dimensions discriminate between clusters (Hill and Lewicki, 2006).
The Global Capacity Fractal Dimension (D0) and the Information Fractal Dimension (D1) Discriminate Well between the Structural Gleason Classes (Grades)
Both the global capacity and information fractal dimension, D0 and D1, have the sufficient statistical power to discriminate between the Gleason classes, that is the subsets of prostate carcinomas with the Gleason score 3+3, 4+4, or 5+5. Figures 1, 2 demonstrate that the dimension D0 was slightly better than the dimension D1 (Figure 1) and much better than the global correlation fractal dimension D2. In that latter case, the area under the curve (AUC) was almost 1.00 for the dimension D0 vs. 0.8 for the dimension D2. In addition, coefficient of correlation between the D0 and D1 was very high R = 0.963. The local fractal dimension, LFD, its variance, the local connected fractal dimension, LCFD, its variance, entropy H, and lacunarity λ were not able to discriminate between those classes at the statistically significant level (data not shown).
Figure 1. The ROC analysis of prostate carcinomas with Gleason score 3+3, 4+4, and 5+5 of the test set demonstrating a discriminating power of the fractal dimensions D0, D1, and D2. AUC stands for the area under the curve, and has values in a similar range as in the case of the reference set.
Figure 2. The ROC analysis of carcinomas with Gleason score 3+3, 4+4, and 5+5 of the test set demonstrating a discriminating power of the local fractal dimensions LFD. AUC stands for the area under the curve, and has values in a similar range as in the case of the reference set.
It is worth to notice that none of the complexity or the diversity measures can discriminate between the spatial distribution of cell nuclei in benign prostatic hyperplasia and in prostate carcinomas of the complexity class C1 containing carcinomas with the Gleason score 3+3 (compare Tables 1, 3, and Figure 2). The mean values overlap (Table 3), and differences between those two groups of tissues are not statistically significant (t-Student test, p > 0.1).
Table 1. Results of ROC analysis for the reference set (A) or the test set (B) of prostate carcinomas and their stratification into the structural Gleason classes.
The Gleason Classes and Complexity Classes are not Identical
Table 1 contains results of ROC analysis for two sets of prostate carcinomas. The first one is a reference set. It contains carefully selected carcinomas with relatively homogeneous structure. Those carcinomas were stratified into three distinct structural classes with the Gleason score 3+3, 4+4, or 5+5. The spatial distribution of cancer cell nuclei in those carcinomas was also characterized quantitatively by the global fractal dimensions D0, D1, and D2. The ROC analysis enabled calculation of the cut-off D0-values that defined the numerical limits of those structural Gleason classes (see Table 1). As it turned out, even very careful histological evaluation of tumor structure brought together cases with a variety of infinitesimal structural alterations that could not be defined unambiguously. In consequence, each structural Gleason class of the reference set contained a low number of carcinomas with the quantitative characteristic typical of the other classes. Those cases were re-stratified using the cut-off D0-values for the structural Gleason classes as a frame of reference. After that operation, the new cut-off D0-values were calculated. They define the complexity classes C1, C4, and C7 (compare Table 1A with Table 2A).
Table 2. Results of the ROC analysis for the re-stratified reference set (A) or the re-stratified test set (B) of prostate carcinomas.
The Cut-off D0-values of the Structural Gleason Classes 3+4, 4+3, 4+5, or 5+4 are Located between the Values for the Structural Gleason Classes 3+3, 4+4, or 5+5
The test set contains prostate carcinomas representing all Gleason scores. Those carcinomas were chosen randomly from the archive (Table 1B). The cut-off D0-values for the structural classes Gleason 3+3, 4+4, or 5+5 in that set are close to those of the reference set if analyzed without the structural classes 3+4, 4+3, 4+5, and 5+4. The same is true of the AUC values (Figures 1, 2) as well as sensitivity and specificity (data not shown). It appears that human eye deals much better with the evaluation of the distinct homogeneous structures rather than with the heterogeneous ones, where the number of details is greater, more difficult to grasp, and to set in order.
If the structural classes 3+4, 4+3, 4+5, and 5+4 are taken into consideration, the cut-off D0-values locate them between the values for the structural classes Gleason 3+3, 4+4, or 5+5. The D0 intervals defining those heterogeneous classes are smaller now (compare with Tables 1A,B). However, both sensitivity and specificity are not equal 1.00. The AUC values are also far from 1.0. This is because all structural classes of the test set including those with seemingly homogeneous carcinomas contain some number of cases with complexity measures typical of the other structural Gleason classes. As expected, those “impure” cases were re-stratified according to the D0-cut-off values, as described in the next paragraph, to the adjacent complexity classes exclusively.
Re-Stratification of Prostate Carcinomas into the Seven Classes of Equivalence
Re-stratification puts aside a problem of accuracy of the subjective evaluation as well as troubles connected with a matching of tumor structure under examination to the grade definition. This was achieved by the application of the cut-off D0-values of the re-stratified reference set (see Table 2A). Carcinomas with the quantitative characteristics that did not match well the characteristics of the other carcinomas of the same class, nor of the adjacent class were identified and clustered in the novel classes of equivalence defined by the novel D0-values. As a result, carcinomas of the test set were re-stratified to the seven classes of complexity, i.e., C1, C2, C3, C4, C5, C6, and C7 (Table 2B). Table 3 provides a percentage of cases of each structural Gleason class that was a subject to re-stratification. This was from 30 to 70% of cases per structural class.
Table 3. The values of the mean, median and standard deviations for the complexity or diversity measures in the re-stratified test set of prostate carcinomas.
It should be emphasized that all complexity classes are unequivocally defined by the cut-off D0-values. Therefore, the complexity classes contain only those prostate carcinomas that have the spatial distribution of cancer cell nuclei characterized by the similar values of the dimension D0 and hold a condition that the AUC in the ROC analysis equals 1.0. As expected, an identical re-stratification can be obtained with D1-values (compare Figure 3B). Each complexity class contains elements of the corresponding structural Gleason class and some carcinomas of the other structural Gleason classes with the D0-values within the range for that class. For example, Figure 3A presents each carcinoma of the test set as a point on the 2D-scatter plot with the D0 and the LFD as variables. Each carcinoma belongs to one out of seven classes of equivalence according to the D0-values. Table 3 shows the statistically significant differences between all complexity classes for the D0, and the D1 (p < 0.0001, t-Student test). This finding confirms a choice of the cut-off D0 or D1-values as a criterion for the re-stratification.
Figure 3. Two dimensional scatter plots showing all prostate carcinomas of the test set stratified according to a number of relationships between the parameters: (A) (D0, LFD), (B) (D1, H), (C) (D0, λ), and (D) (H, λ). Differences between the complexity classes defined by the D0 or the D1 are statistically significant at p < 0.0001 (Table 3). (A) The LFD divides the complexity classes into the five subsets (clusters 1–5 are marked by the numbers; compare Figure 4A). Those clusters contain carcinomas with different values of the LFD, and, therefore, with different intercellular connectivity. In particular, 10 carcinomas of the class C1 (black circles) and four carcinomas of the class C2 (red circles) compose a cluster two with the lowest values of both the D0 and the LFD (Figure 4A), and the lowest aggressiveness. The values are close to those of benign prostatic hyperplasia. The cluster five represents high-grade carcinomas with the largest values of both the D0 and the LFD, the largest complexity, and the largest tumor aggressiveness (red squares and blue diamonds). (B) The dimension D1 is as strong classifier as the D0, and divides the set of prostate carcinomas in the same classes of equivalence. Although the values of the Shannon entropy H overlap between the classes, and, therefore, cannot play a role of the independent classifier, they characterize the spatial distribution of cancer cell nuclei in each class unequivocally if coupled with the D0 or D1-dimension. The coefficient of correlation for this relationship is 0.950. (C) The values of lacunarity λ overlap between the classes of equivalence; however, a combination of this diversity measure with the D0 or D1-dimension allows the unequivocal classification of the underlying spatial distribution of cancer cell nuclei. The coefficient of correlation for this relationship is −0.902. The statistical linearity of this relationship confirms in the independent way the existence of fractal structure in the spatial distribution of cancer cell cell nuclei. (D) Neither the Shannon entropy H nor lacunarity λ are strong classifiers. The values overlap between the classes of equivalence. A small subset of carcinomas of the class C7 overlaps with carcinomas of the class C5. This results from a lower number of cells in cellular infiltrates. However, this relationship is tri-linear. The coefficients of correlation for this relationship are −0.911, −0.650, and −0.922 for H < 0.7, 0.7 < H < 0.9, and H > 0.9, respectively. The existence of tri-linearity suggests that some kind of qualitative transitions occur in natural history of prostate cancer between the class C2 and C3 as well as the class C5 and C6.
A Relationship between the Global Capacity Dimension D0 and Lacunarity λ or the Global Information Dimension D1 and Entropy H or is Linear
The relationship between the D0 and λ is given by linear equation (8):
The coefficient of correlation R for that relationship was -0.902. As it can be seen in Figure 3C, a given value of lacunarity λ may be associated with different values of the D0 dimension and different corresponding spatial distributions of cancer cell nuclei that belong to the neighbor classes of equivalence. Hence, the λ value is not a good classifier. However, if this parameter is coupled with the D0-value, it characterizes the spatial distribution of cancer cell nuclei unequivocally.
A similar capability possesses a pair (D1, H). As expected, there was a linear relationship between the global information fractal dimension D1 and entropy H given by the following statistical equation (9)
Figure 3B demonstrates a distribution of carcinomas around the linear curve defined by Equation (9). The coefficient of correlation R for that relationship was 0.950. The linear relationship between the global information dimension as a complexity measure of cell nuclei clustering and the Shannon entropy as a diversity measure of information content validates in the independent manner the choice of the classes of equivalence. The D1 is more accurate classifier than the Shannon entropy H (Figure 3B). Some carcinomas of the class C3, C4, C5, C6, or C7 have entropy values that overlap with the values of the cases belonging to the adjacent classes. In particular, carcinomas of the class C7 are composed of cellular infiltrates only. Some of them have, however, a much lower number of infiltrating cells. Entropy of those carcinomas is lower, respectively (see Figure 3B, Class C7, light blue diamond). The values of lacunarity for those cases overlap with those for carcinomas of the class C5 (see Table 3 and Figure 3D).
A statistical relationship between lacunarity λ and the Shannon entropy H is more complex. This relationship is described by three linear equations that are defined for different intervalls of the H-values. The first one is given by equation (10)
and holds for H = < 0.7. The coefficient of correlation R = −0.911.
The second one is given by Equation (11)
and holds for 0.7 < H < 0.9. The coefficient of correlation R = −0.650.
The third one is given by Equation (12)
and holds for H > 0.9. The coefficient of correlation R = −0.901. Those relationships characterize well only the spatial distribution of cancer cell nuclei in prostate carcinomas with very low or very high diversity as measured by the Shannon entropy. Prostate carcinomas with intermediate diversity cannot be stratified with sufficient accuracy (compare Figures 3D,4B). The values of the pairs (H, λ) overlap for many cases that belong to different classes of equivalence defined by the D0 or the D1 (Figure 3D). Figure 3A shows that the spatial distributions of cancer cell nuclei in many prostate carcinomas are unique as represented by different values of the pair (D0, LFD) (Figure 3A). Those values are not identical owing to multifractality. Simultaneously, many of those carcinomas have the same or almost identical values of the pair (D1, H) (Figure 3B) or (H, λ) (Figure 3D). Two different fractal dimensions denote the existence of two different spatial distributions of cell nuclei. However, if the difference is not large, the Shannon entropy H or lacunarity λ may be identical owing to variability in cell numbers or their clustering.
Figure 4. Results of cluster analysis for the test set of prostate carcinomas with the dimensions D0 and LFD (A) (ANOVA p < 0.001, FD0 45799, FLFD 31935) and the parameters D0, LFD, LCFD, H, λ (B) (ANOVA p < 0.001, FD0 75800, FLFD 25289, FLCFD 13202, FH 7925, Fλ 16102). The increment of the number of parameters from two to five reduces a number of clusters from seven to three. The objective analysis can identify in that way carcinomas with low-, intermediate-, and high complexity of the spatial distribution of cancer cell nuclei corresponding to low-, intermediate- or high tumor aggressiveness. The following conditions can be defined: for D0 < 1.5820, LFD < 1.3, LCFD > 1.5, H < 0.7 and λ > 0.8, the class C1 or C2 contains low-complexity, low-grade carcinomas exclusively; for D0 > 1.6980, LFD > 1.7644, LCFD > 1.7051, H > 0.9, and λ < 0.7, the class C6 or C7 contains high-complexity, high-grade carcinomas only.
The Co-Application of the Complexity and Diversity Measures in the Evaluation of the Spatial Distribution of Cancer Cell Nuclei in Prostate Carcinomas
The cluster analysis by the k-means algorithm shows that the co-application of both complexity and diversity measures, D0 (D1), LFD, LCFD, H, and λ leads to a reduction of a number of clusters from seven to three (Figure 4B). Since low grade prostate carcinomas possess well-preserved glandular structure, the values of all the above measures but λ are also low. The reduction of the number of clusters identifies some carcinomas of the complexity classes C1 and C2 as elements of the first cluster with the lowest values of those parameters, and, therefore, the lowest tumor aggressiveness. Similarly, one can define the quantitative criteria for carcinomas with the highest values of those parameters but λ, and the highest aggressiveness.
The mean LFD values do not increase along the complexity classes monotonically (see Figure 3A and Table 3). Rather, those values increase along each stratum, that is, a subset composed of the classes C1 and C2 or C3 and C4, or C5 and C6, or C7 (Table 3). The mean LFD values can be identical for carcinomas of different complexity classes.
The cluster analysis indicates five such subsets of prostate carcinomas characterized by the pair (D0, LFD), in which the D0 is a measure of the global spatial distribution of cancer cell nuclei and the local fractal dimension LFD is a measure of intercellular connectivity (see Figures 3A,4A).
In the first subset, characterized by the parameters D0 < 1.5820 and LFD < 1.35, fourteen prostate carcinomas of the class C1 or C2 (low grade prostate carcinomas with mean PSA < 6.1 ng/ml, standard deviation 1.2 ng/ml, pT1b, pT1c, or pT2a, Gleason score 6, relatively homogeneous well-preserved glandular structure with glands distributed in tissue sparsely, without a clustering or crowding, none of those patients developed recurrence or metastases, or died owing to prostate cancer during the 5 years of follow-up) revealed very high intercellular connectivity. Those features are reflected by the value of the local connected fractal dimensions LCFD close to 1.5 (the scatter plot not shown), low values of entropy H < 0.7 and high values of lacunarity λ > 0.8 (compare Figures 3A–D). The low values of the fractal dimensions indicate that both global and local complexity in those carcinomas is low. In addition, the values of those fractal dimensions are quite close to the values obtained for benign prostatic hyperplasia (see Table 3). All those local and global features of the spatial distribution of cancer cell nuclei are associated with a low potential for growth and aggressiveness.
The second subset is composed of the low risk adenocarcinomas (mean PSA 8.3 ng/ml, standard deviation 2.6 ng/ml, Gleason score 6, pT1b, pT1c, or pT2a) stratified also to the classes C1 or C2. However, the second subset has the same values of the D0, D1, almost the same values of entropy H < 0.78 as well as lacunarity λ > 0.78, and 1.300 < LFD < 1.750. The values of the first four parameters denote that number of cells and their spatial configurations between those two sets are comparable. However, the values of the LFD indicate the existence of more advanced alterations in intercellular connectivity among adenocarcinomas of the second subset in spite of the preserved glandular structure, and, therefore, denote a higher risk of progression. Those LFD values denote very slow, yet active growth owing to the lower intercellular connectivity, and, therefore, increased probability of metastasizing in the course of time in comparison with the first subset. Those cancer patients belong to a “gray” zone, in which a follow-up without any additional therapy may work fine for some cases of the lower classes of complexity, such as C1, C2, or even for some intermediate-risk cases of the class C3 and C4 having the low values of the LFD.
The third and the fourth subset comprise intermediate risk prostate cancers of the complexity class C3, C4, C5 (mean PSA 12.7 ng/ml, standard deviation 2.2 ng/ml, Gleason score 7 and 8, pT2b-2c). In the third subset, there are carcinomas with 1.5820 < D0 < 1.6980 and the LFD < 1.8. Prostate carcinomas of the fourth subset have the same values of the D0 as in the third set, but the LFD > 1.8 (see Figure 3A) at 0.65 < H < 0.85, and 0.7 < λ < 0.82 (Figures 3B,C). Surprisingly, there are two adenocarcinomas of the complexity class C2 in that subset. Their mean LFD value approximate the value of 1.9, that is typical of prostate carcinomas without glandular structure and with the increased potential to metastasize. In spite of the stratification to the low class of complexity owing to the D0-value, such carcinomas with very low intercellular connectivity may have a high metastatic potential.
The fifth subset comprises high risk prostate cancers (mean PSA 24.3 ng/ml, standard deviation 3.1 ng/ml, Gleason score 9–10, stage T3a, the largest ratio of prostate carcinoma related death or metastasis formation within 5 years from the surgical treatment 0.442) of the complexity class C6 or C7 with D0 > 1.6980 and the LFD > 1.8, H > 0.9, and λ < 0.7. The large values of the D0 and H indicate a significantly increased cell number owing to the intense cell proliferation and rapid growth. Since images are covered almost completely by cancer cells, and become almost symmetrical, the λ-values are low. The high values of the LFD indicate the existence of a very low intercellular connectivity. Cells fill up the entire available tissue space without any clustering, and, therefore, have the largest potential for metastasis formation.
The classes of equivalence defined by the cut-off values of the global capacity dimension D0 reported previously (Waliszewski et al., 2015) have been validated. The other measures but the D1 failed to define the same classes of equivalence. However, the co-application of all the measures allowed the formulation of two quantitative criteria identifying low or high aggressive prostate carcinomas.
A novel, quantitative strategy for the evaluation of tumor aggressiveness excludes the automated diagnosis of prostate carcinoma. The complexity or diversity measures were not able to discriminate between the spatial distribution of cell nuclei in benign prostatic hyperplasia, a common lesion in human prostate, and carcinomas with a well-preserved glandular structure (see Table 3). This approach implies a necessity to examine the entire prostate carefully, what does not seem to be the case in the routine practice (True, 1994). Indeed, the largest value of the fractal dimension D0 or D1 determines the stratification of carcinomas into a given class of equivalence.
Results of this study suggest that the number of the subjective Gleason grades of tumor aggressiveness could be decreased from seven to three without a loss of clinically important information, that is, to the low, intermediate, or high grade. The identical reduction could be done for a number of the objective classes of equivalence. Table 1 indicates that pathologists can distinguish quite well between the relatively homogeneous patterns with the Gleason score 3+3, 4+4, or 5+5 (see the reference set). However, heterogeneous patterns defined as the Gleason score 3+4, 4+3, 4+5, or 5+4 were classified with much greater inaccuracy (see the test set). The classes of equivalence defined by the cut-off D0-values enable the quantitative evaluation of tumor aggressiveness of borderline cases with mixed tissue architecture; a frequent problem in pathology of prostate carcinoma (Dong et al., 2012). It should be noticed that the discrepancy between the subjective Gleason classification and the classification according to complexity of the spatial distribution of cancer cell nuclei was on average about 53% in this study; a value that is in the range of the intra- and interobserver variability (McLean et al., 1997; Allsbrook et al., 2001; Nguyen et al., 2004; van der Kwast et al., 2010; McKenney et al., 2011; Netto et al., 2011; Dong et al., 2012; Egevad et al., 2013; Scott Lucia et al., 2013; Berney et al., 2014).
The classes of equivalence, i.e., the complexity classes defined by the cut-off values of the D0 dimension (or the D1 dimension) were validated by analyzing the test set of carcinomas with sensitivity, specificity, and AUC equal 1.0 (Waliszewski et al., 2015). This set was larger and independent from the set analyzed in the previously published study (Waliszewski et al., 2015). The values of the global capacity fractal dimension D0 are identical with those previously reported (Waliszewski et al., 2015) (see Tables 2, 3, Figure 3). Those two dimensions appear to be the best classifiers of prostate carcinomas according to the spatial distribution of cancer cell nuclei. There is a very good correlation between the values of the D0 and the D1 (correlation coefficient R = 0.962). Those results were to be expected since both sets comprised different numbers of carcinomas with the similar patterns of tumor growth. As expected, the D1-values correlate very well with the values of the Shannon entropy H; a diversity measure of both randomness in the spatial distribution of cell nuclei and information content in images of a given class (R = 0.950; see Figure 3B). There also is some correlation between the values H and λ that changes depending on the H intervals (compare Equations 10–12). The λ values overlap between classes (see Figure 3D), and, therefore, cannot be applied for stratification as the independent classifier. The overlapping results from the non-bijective relationship between some configurations of cell nuclei and the measures applied (Waliszewski et al., 2015). In particular, two borderline spatial distributions of cancer cell nuclei may have identical values of the complexity measures, but different values of the diversity measures. In spite of that weakness, the diversity measures help to identify carcinomas with low or high diversity corresponding to carcinomas with the low or high grade (see Figures 3B,C). Each of those measures characterizes quantitatively different aspects of the spatial distribution of cancer cell nuclei in images of prostate carcinomas. While the D0 or the D1 enables the stratification of carcinomas into the seven classes of equivalence, that remain unchanged if the LFD, the LCFD, the Shannon entropy H or λ are added to the class characteristics (see Figure 3), the co-application of all the above-studied measures allows the stratification of the same carcinomas into just three distinct clusters with low tumor aggressiveness (class C1, C2, and some carcinomas of the class C3), intermediate one (some carcinomas of the class C3, C4, and C5) or high one (class C6 and C7) (compare Figure 4). The appropriate quantitative criteria are defined in the Sections Results: A Relationship between the Global Capacity Dimension D0 and Lacunarity λ or the Global Information Dimension D1 and Entropy H or is Linear; The Co-Application of the Complexity and Diversity Measures in the Evaluation of the Spatial Distribution of Cancer Cell Nuclei in Prostate Carcinomas.
The presented data indicate some change in intercellular connectivity at the interface between classes C2 and C3 as well as C5 and C6 as measured by the LFD; a complexity measure of cellular connectivity, that is, a measure of strength of intercellular interactions (Table 3, Figure 3A). This change can also be seen in the scatter plot (H, λ) (Figure 3D). The relationship between the Shannon entropy H and lacunarity λ; a diversity measure of cell clustering also changes at the interface between the above-mentioned classes and depends on the H-values (R = −0.911 for H < 0.7, R = −0.650 for 0.7 < H < 0.9, and − 0.901 for H > 0.9). Those changes correspond with the enhanced dynamics of tumor growth at the interface between low and intermediate risk cancers as well as intermediate and high risk cancers. It is possible that all those changes result from some genetic instability that gains influence on growth dynamics at some specific check-points in natural history of prostate cancer. Second, results of the analysis point out that prostate carcinomas with the spatial distributions of cancer cell nuclei holding the conditions D0 < 1.5820, LFD < 1.3, and LCFD > 1.5 are true low grade carcinomas. Those carcinomas possess a well-preserved, relatively homogeneous glandular structure that resembles the structure of the normal prostatic glands. There is some intermediate cell clustering at the low values of the LFD and, therefore, high intercellular connectivity in cell populations of those adenocarcinomas (Figure 3A). Hence, there is no high growth potential in those cellular populations. If D0 > 1.6980, LFD > 1.7644, and LCFD > 1.7051, then tumor aggressiveness is large. Cancer cells infiltrate the entire available tissue space, have low intercellular connectivity, do not interact each other, and compose passively multiple cell clusters; the overcrowding effect. This constellation of dimensions indicates the existence of a large growth potential in a cellular population. The cases with intermediate conditions compose a gray zone. Some of those carcinomas may do well under follow-up, especially those with the lower values of the LFD and LCFD. The other cases will certainly require some kind of therapy in the course of time. It is interesting, however, that there is no carcinoma with the low values of both local fractal dimensions. If the LFD increases, then intercellular connectivity decreases, and cells lose their capability to interact with each other. Therefore, cells in some carcinomas including those of the class C1, C2, C3, and C4 do not create clusters (the lower values of the LCFD), while many others still do (the larger values of the LCFD). The loss of capability for cell clustering at low values of the LCFD may also denote both the increased growth potential and risk of metastasis formation.
The increasing values of the D0 dimension denote that cellular proliferation in a given carcinoma is very extensive, and more and more cancer cells fill up the available tissue space. If cancer cell gather in some kind of glands or gland-like objects, the D0-values are lower than if there dominate cellular infiltrates in the image (compare the values of the D0 for the classes C1, C2 or C3 with those for the classes C6 or C7 in Table 3). The larger is the value of the D0 the higher is the complexity class and the corresponding subjective grade. The higher is the complexity class the higher is a number of cancer cells present in the image of prostate carcinoma. Cancer cells lose in the course of natural tumor evolution their cluster organization in a form of glands (high values of the D1, the LFD, and the LCFD and low values of λ), become more randomly distributed (high values of H), and reveal lower connectivity (high values of the LFD and the LFCD) (compare Figures 3A–D). The risk of local progression or metastasis formation is increased in those spatial distributions of cancer cell nuclei.
The existence of both the global and local fractal dimensions for a given spatial distribution of cancer cell nuclei implies that that spatial distribution underlies a scale-dependent power law. The existence of the power law means that the spatial distribution of cancer cell nuclei undergoes somehow ordered rather than utterly random changes in the course of natural history of the prostatic tumorigenesis. That order of events appears to comprise the gradual loss of intercellular connectivity in carcinomas of the higher complexity classes (see Table 3 and Figure 3). This may occur owing to the loss of expression of adhesion molecules or some other molecules with biologically important functions leading to alterations of intercellular interactions (Busch et al., 2002; Huang et al., 2009; Mathieu et al., 2013; Schrecengost and Knudsen, 2013). For the majority of carcinomas, the values of both the global and local fractal dimension are not identical (see Table 3 and Figure 3). That difference between the values of the global and local fractal dimensions indicates the existence of multifractals in the spatial distributions of cancer cell nuclei, that is, fractals that scale with multiple scaling rules and hold the condition D0 > D1 > D2 (reviewed in Lopes and Betrouni, 2009). Multifractality suggests that the final spatial distribution of cancer cells depends on both a variety of molecular events and long-range phenomena such as diffusion. The interactions may occur at different levels of the hierarchical tissue dynamic system (Waliszewski, 1997; Waliszewski et al., 1998, 2001). Those complex phenomena compose a new world to be explored.
Two diversity measures CV-LFD and CV-LCDF do not play a role in the stratification process (see Table 3). Those diversity measures apply to differences in measurements of the local fractal dimensions within each complexity class (Page, 2011). There is no monotonic increment of the values of those parameters that might enable the stratification process. However, entropy, the other diversity measure, appears to be useful in the additional characterization of the spatial distribution of cell nuclei (Equations 2, 4, 5 and Figure 3B). Even though the local fractal dimensions in each complexity class that relate to single cancer cell nuclei are not diverse, cancer cells as a population form diverse dynamic spatial structures of the higher order, such as cell clusters (see Figure 3A). Those cell clusters can be characterized locally by the values of the local connected fractal dimension, entropy, and lacunarity. In addition, the LFD identifies cells with the lowest intercellular connectivity, and, therefore, having the maximal metastatic potential (see Figure 3A; Waliszewski, 1997, 2005, 2009, 2013; Waliszewski et al., 1998, 2001; Waliszewski and Konarski, 2001). Computer simulations of a model considering the evolution of cooperation in a spatial setting obtained by Nowak and May show clearly the existence of a relationship between emergent diversity and complexity (Nowak and May, 1992, 1993). The existence of that relationship in prostate carcinomas can be seen in equation 5 and Figure 3B.
To summarize, the following findings must be reiterated. First, the classes of equivalence for prostate carcinomas by the cut-off D0-values were validated. Second, no other measure studied but the D1 can define the same classes of equivalence. Third, carcinomas with low complexity of the spatial distribution of cancer cell nuclei reveal low biological aggressiveness, and vice versa. Some pairs of the parameters (D0, LFD), (D0, H) or (D0, λ) or (D1, LFD), (D1, H) or (D1, λ) characterize well complexity of the spatial distribution of cancer cell nuclei in each class of equivalence. The tri-linear course of the relationship (H, λ) suggests that among carcinomas of the classes C2 and C3 as well as C5 and C6 occurs some kind of a qualitative transition that implicates changes in tumor aggressiveness. Fourth, the co-application of all complexity and diversity measures reduces a number of clusters, but still enables the identification of carcinomas with low or high complexity of the spatial distribution of cancer cell nuclei. Those findings suggest the number of both the subjective Gleason grades and the objective classes of equivalence could be decreased from seven to three.
PW: all ideas, construction of the prostate data bank, digitalization of images, statistical analyses, writing of the manuscript.
Conflict of Interest Statement
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The author is grateful to Professor Michal Banaszak of Department of Physics, Adam-Mickiewicz University in Poznan, Poland and Professor John T. Isaacs, Department of Urology, The Johns Hopkins School of Medicine, Baltimore, USA for a discussion, minor suggestions, and editing.
Allsbrook, W. C. Jr., Mangold, K. A., Johnson, M. H., Lane, R. B., Lane, C. G., and Epstein, J. I. (2001). Interobserver reproducibility of Gleason grading of prostatic carcinoma: general pathologist. Hum. Pathol. 32, 81–88. doi: 10.1053/hupa.2001.21135
Berney, D. M., Algaba, F., Camparo, P., Compérat, E., Griffiths, D., Kristiansen, G., et al. (2014). The reasons behind variation in Gleason grading of prostatic biopsies: areas of agreement and misconception among 266 European pathologists. Histopathology 64, 405–411. doi: 10.1111/his.12284
Boccon-Gibod, L. M., Dumonceau, O., Toublanc, M., Ravery, V., and Boccon-Gibod, L. A. (2005). Micro-focal prostate cancer: a comparison of biopsy and radical prostatectomy specimen features. Eur. Urol. 48, 895–899. doi: 10.1016/j.eururo.2005.04.033
Busch, T. A., Hanssen, T. A., Wagener, C., and Oebrink, B. (2002). Down-regulation of CEACAM1 in human prostate cancer: correlation with loss of cell polarity, increased proliferation rate, and Gleason grade 3 to 4 transition. Hum. Pathol. 33, 290–298. doi: 10.1053/hupa.2002.32218
Cooperberg, M. R., Cowan, J. E., Hilton, J. F., Reese, A. C., Zaid, H. B., Porten, S. P., et al. (2011). Outcomes of active surveillance for men with intermediate-risk prostate cancer. J. Clin. Oncol. 29, 228–234. doi: 10.1200/JCO.2014.59.2329
Dall'Era, M. A., Konety, B. R., Cowan, J. E., Shinohara, K., Stauf, F., Cooperberg, M. R., et al. (2008). Active surveillance for the management of prostate cancer in a contemporary cohort. Cancer 112, 2664–2670. doi: 10.1002/cncr.23502
D'Amico, A. V., Chen, M. H., Roehl, K. A., and Catalona, W. J. (2004). Preoperative PSA velocity and the risk of death from prostate cancer after radical prostatectomy. N. Engl. J. Med. 351, 125–135. doi: 10.1056/NEJMoa032975
D'Amico, A. V., Chen, M. H., Roehl, K. A., and Catalona, W. J. (2005). Identifying patients at risk for significant versus clinically insignificant postoperative prostate-specific antigen failure. J. Clin. Oncol. 23, 4975–4979. doi: 10.1200/JCO.2005.08.904
D'Amico, A. V., Wu, Y., Chen, M. H., Nash, M., Renshaw, A. A., and Richie, J. P. (2000). Pathologic findings and prostate specific antigen outcome after radical prostatectomy for patients diagnosed on the basis of a single microscopic focus of prostate carcinoma with a Gleason score </=7. Cancer 89, 1810–1817. doi: 10.1002/1097-0142(20001015)89:8<1810::AID-CNCR22>3.0.CO;2-9
Daskivich, T. J., Lai, J., Dick, A. W., Setodji, C. M., Hanley, J. M., Litwin, M. S., et al. (2014). Variation in treatment associated with life expectancy in a population-based cohort of men with early-stage prostate cancer. Cancer 120, 3642–3650. doi: 10.1002/cncr.28926
Dong, F., Wang, C., Farris, A. B., Wu, S., Lee, H., Olumi, A. F., et al. (2012). Impact on the clinical outcome of prostate cancer by the 2005 international society of urological pathology modified Gleason grading system. Am. J. Surg. Pathol. 36, 838–843. doi: 10.1097/PAS.0b013e3182486faf
Doyle, S., Hwang, M., Shah, K., Madabhushi, A., Feldman, M., and Tomaszeweski, J. (2007). “Automated grading of prostate cancer using architectural and textural image features,” in Proceedings of IEEE ISBI (Washington, DC: IEEE), 1284–1287. doi: 10.1109/ISBI.2007.357094
Egevad, L., Ahmad, A. S., Algaba, F., Berney, D. M., Boccon-Gibod, L., Compérat, E., et al. (2013). Standardization of Gleason grading among 337 European pathologists. Histopathology 62, 247–256. doi: 10.1111/his.12008
ElShafei, A., Li, Y. H., Hatem, A., Moussa, A. S., Ethan, V., Krishnan, N., et al. (2013). The utility of PSA velocity in prediction of prostate cancer and high grade cancer after an initially negative prostate biopsy. Prostate 73, 1796–1802. doi: 10.1002/pros.22718
Epstein, J. I., Allsbrook, W. C. Jr., Amin, M. B., Egevad, L. L., and the ISUP Grading Committee (2005). The 2005 International Society of Urological Pathology (ISUP) consensus conference on Gleason grading of prostatic carcinoma. Am. J. Surg. Pathol. 29, 1228–1242. doi: 10.1097/01.pas.0000173646.99337.b1
Fine, S. W., Amin, M. B., Berney, D. M., Bjartell, A., Egevad, L., Epstein, J. I., et al. (2012). A contemporary update on pathology reporting for prostate cancer: biopsy and radical prostatectomy specimens. Eur. Urol. 62, 20–39. doi: 10.1016/j.eururo.2012.02.055
Heidenreich, A., Bastian, P. J., Bellmunt, J., Bolla, M., Joniau, S., van der Kwast, T., et al. (2014). European association of urology. EAU guidelines on prostate cancer. Part II: Treatment of advanced, relapsing, and castration-resistant prostate cancer. Eur. Urol. 65, 467–479. doi: 10.1016/j.eururo.2013.11.002
Hong-Jun, Y., Ching-Chung, L., Christudass, C., Veltri, R., Epstein, J. I., and Zhen, Z. (2011). “Cardinal multiridgelet-based prostate cancer histological image classification for Gleason grading,” in Proceedings IEEE BIBM (Atlanta, GA: IEEE), 315–320. doi: 10.1109/BIBM.2011.59
Klein, E. (2014). “Prostate cancer: risk stratification and choice of initial treatment,” in ed M. E. Ross. Available online at: http://www.uptodate.com/contents/prostate-cancer-risk-stratification-and-choice-of-initial-treatment (Accessed May 20, 2015).
Klotz, L., Vesprini, D., Sethukavalan, P., Jethava, V., Zhang, L., Jain, S., et al. (2014). Long-term follow-up of a large active surveillance cohort of patients with prostate cancer. J. Clin. Oncol. 33, 272–277. doi: 10.1200/JCO.2014.55.1192
Klotz, L., Zhang, L., Lam, A., Nam, R., Mamedov, A., and Loblaw, A. (2010). Clinical results of long-term follow-up of a large, active surveillance cohort with localized prostate cancer. J. Clin. Oncol. 28, 126–131. doi: 10.1200/JCO.2009.24.2180
Kristiansen, G., and Epstein, J. I. (2014). Immunohistochemistry in Prostate Pathology. DAKO. Available online at http://www.dako.com/de/ihc-prostate-pathology.pdf
Loeffler, M., Greulich, L., Scheibe, P., Kahl, P., Shaikhibrahim, Z., Braumann, U. D., et al. (2012). Classifying prostate cancer malignancy by quantitative histomorphometry. J. Urol. 187, 1867–1875. doi: 10.1016/j.juro.2011.12.054
Lopez, M., and Again, S. (2013). “A new set of wavelet- and fractals-based features for Gleason grading of prostate cancer histopathology images,” in Proceedings SPIE-IS&T Electronic Imaging (Burlingame, CA), 8655:865516.
Lu-Yao, G. L., Albertsen, P. C., Moore, D. F., Shih, W., Lin, Y., Di Paola, R. S., et al. (2014). Fifteen-year survival outcomes following primary androgen-deprivation therapy for localized prostate cancer. JAMA Int. Med. 174, 1460–1467. doi: 10.1001/jamainternmed.2014.3028
Mathieu, R., Evrard, B., Fromont, G., Rioux-Leclercq, N., Godet, J., Cathelineau, X., et al. (2013). Expression screening of cancer/testis genes in prostate cancer identifies nr6a1 as a novel marker of disease progression and aggressiveness. Prostate 73, 1103–1114. doi: 10.1002/pros.22659
McKenney, J. K., Simko, J., Bonham, M., True, L. D., Troyer, D., Hawley, S., et al. (2011). Canary/early detection research network prostate active surveillance study investigators. the potential impact of reproducibility of gleason grading in men with early stage prostate cancer managed by active surveillance; a multi-institutional study. J. Urol. 186, 465–469. doi: 10.1016/j.juro.2011.03.115
McLean, M., Srigley, J., Banerjee, D., Warde, P., and Hao, Y. (1997). Interobserver variation in prostate cancer Gleason scoring: are there implications for the design of clinical trials and treatment strategies? Clin. Oncol. 9, 222–225.
Miller, D. C., Gruber, S. B., Hollenbeck, B. K., Montie, J. E., and Wei, J. T. (2006). Incidence of initial local therapy among men with lower-risk prostate cancer in the United States. J. Natl. Cancer Inst. 98, 1134–1141. doi: 10.1093/jnci/djj308
Montironi, R., Lopez-Bertran, A., Cheng, L., Montorsi, F., and Scarpelli, M. (2013). Central prostate pathology review: should it be mandatory? Eur. Urol. 64, 199–203. doi: 10.1016/j.eururo.2013.04.002
Mottet, N., Bellmunt, J., Briers, E., van den Bergh, R. C. N., Bolla, M., van Casteren, N. J., et al. (2015). Guidelines on Prostate Cancer. European Association for Urology, 30–97. Available online at: http://uroweb.org/wp-content/uploads/09-Prostate-Cancer_LR.pdf (Accessed January 15, 2016).
Naik, S., Doyle, S., Agner, S., Madabhushi, A., Feldman, M., and Tomaszeweski, J. (2008). “Automated gland and nuclei segmentation for grading of prostate and breast cancer histopathology,” in Proceedings of IEEE ISBI (Paris: IEEE), 284–287. doi: 10.1109/ISBI.2008.4540988
National Comprehensive Cancer Network (NCCN) (2015). NCCN Clinical Practice Guidelines in Oncology. Available online at: http://www.nccn.org/professionals/physician_gls/f_guidelines.asp (Accessed May 20, 2015).
Netto, G. J., Eisenberger, M., and Epstein, J. I. (2011). Interobserver variability in histologic evaluation of radical prostatectomy between central and local pathologists: findings of TAX 3501 multinational clinical trial. Urology 77, 1155–1160. doi: 10.1016/j.urology.2010.08.031
Nguyen, P. L., Schultz, D., Renshaw, A. A., Volmer, R. T., Welch, W. R., Cote, K., et al. (2004). The impact of pathology review on treatment recommendations for patients with adenocarcinoma of the prostate. Urol. Oncol. Semin. Orginal Invest. 22, 295–299. doi: 10.1016/S1078-1439(03)00236-9
Partin, A. W., Kattan, M. W., Subong, E. N. P., Walsh, P. C., Wojno, K. J., Oesterling, J. E., et al. (1997). A combination of prostate-specific antigen, clinical stage, and Gleason score to predict pathological stage of localized prostate cancer. A multi-institutional Update. JAMA 277, 1445–1451.
Pin, E., Fredolini, C., and Petricoin, E. F. III. (2013). The role of proteomics in prostate cancer research: biomarker discovery and validation. Clin. Biochem. 46, 524–538. doi: 10.1016/j.clinbiochem.2012.12.012
Scott Lucia, M., Bostwick, D. G., Somerville, M. C., Fowler, I. L., and Rittmaster, R. S. (2013). Comparison of classic and international society of urological pathology 2005 modified gleason grading using needle biopsies from the reduction by dutasteride of prostate cancer events (REDUCE) Trial. Archiv. Pathol. Lab. Med. 137, 1740–1746. doi: 10.5858/arpa.2012-0447-OA
Stamey, T. A., Kabalin, J. N., McNeal, J. E., Johnstone, I. M., Freiha, F., Redwine, E. A., et al. (1989). Prostate specific antigen in the diagnosis and treatment of adenocarcinoma of the prostate: II. Radical prostatectomy treated patients. J. Urol. 141, 1076–1083.
Tabesh, A., Kumar, V. P., Pang, H. Y., Verbel, D., Kotsianti, A., Teverovskiy, M., et al. (2005). “Automated prostate cancer diagnosis and Gleason grading of tissue microarrays,” in Proceedings SPIE, Vol. 5747 (San Diego, CA), 58–70. doi: 10.1117/12.597250
Tabesh, A., Teverovskiy, M., Pang, H. Y., Kumar, V. P., Verbel, D., Kotsianti, A., et al. (2007). Multifeature Prostate cancer diagnosis and gleason grading of histological images. IEEE Trans. Med. Imaging 26, 1366–1378. doi: 10.1117/12.597250
Thompson, I. M., Ankerst, D. P., Chi, C., Goodman, P. J., Tangen, C. M., Lucia, M. S., et al. (2006). Assessing prostate cancer risk: Results from the prostate cancer prevention trial. J. Natl. Cancer Inst. 98, 529–534. doi: 10.1093/jnci/djj131
van der Kwast, T. H., Evans, A., Lockwood, G., Tkachuk, D., Bostwick, D. G., Epstein, J. I., et al. (2010). Variability in diagnostic opinion among pathologists for single small atypical foci in prostate biopsies. Am. J. Surg. Pathol. 34, 169–177. doi: 10.1097/PAS.0b013e3181c7997b
Waliszewski, P. (2009). A principle of fractal-stochastic dualism, couplings, complementarity and growth. CEAI 11, 45–52. Available online at: http://www.ceai.srait.ro
Waliszewski, P., Molski, M., and Konarski, J. (1998). On the holistic approach in cellular and cancer biology: nonlinearity, complexity, and quasi-determinism of the dynamic cellular network. J. Surg. Oncol. 68, 70–78.
Waliszewski, P., Molski, M., and Konarski, J. (2001). On the relationship between fractal geometry of space and time in which a system of interacting cells exists and dynamics of gene expression. Acta Biochim. Pol. 48, 209–220.
Waliszewski, P., Wagenlehner, F., Gattenloehner, S., and Weidner, W. (2015). On the Relationship between tumor structure and complexity of the spatial distribution of cancer cell nuclei; a fractal geometrical model of prostate carcinoma. Prostate 75, 399–414. doi: 10.1002/pros.22926
Waliszewski, P., Wagenlehner, F., Kribus, S., Schafhauser, W., Weidner, W., and Gattenlöhner, S. (2014). Objektivierung des tumorgradings bei prostatakarzinomen anhand der globalen und lokalen fraktalen dimensionen. Urol. A 53, 1504–1511. doi: 10.1007/s00120-014-3470-z
Womble, P. R., Montie, J. E., Ye, Z., Linsell, S. M., Lane, B. R., Miller, D. C., et al. (2015). A contemporary use of initial active surveillance among men in Michigan with low-risk prostate cancer. Eur. Urol. 67, 44–50. doi: 10.1016/j.eururo.2014.08.024
Wu, J. N., Fish, K. M., Evans, C. P., Devere White, R. W., and Dall'Era, M. A. (2014). No improvement noted in overall or cause-specific survival for men presenting with metastatic prostate cancer over a 20-year period. Cancer 120, 818–823. doi: 10.1002/cncr.28485
Keywords: fractals, complexity, Renyi dimensions, entropy, lacunarity, tumor aggressiveness, grading, prostate carcinoma
Citation: Waliszewski P (2016) The Quantitative Criteria Based on the Fractal Dimensions, Entropy, and Lacunarity for the Spatial Distribution of Cancer Cell Nuclei Enable Identification of Low or High Aggressive Prostate Carcinomas. Front. Physiol. 7:34. doi: 10.3389/fphys.2016.00034
Received: 13 October 2015; Accepted: 25 January 2016;
Published: 11 February 2016.
Edited by:John M. Beggs, Indiana University, USA
Reviewed by:Wentian Li, Feinstein Institute for Medical Research, North Shore LIJ Health System, USA
Chunhua Bian, Nanjing University, China
Copyright © 2016 Waliszewski. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Przemyslaw Waliszewski, email@example.com
This work was presented in part to the 15th Central European Meeting of the European Association of Urology in Budapest, Hungary on October 2–4, 2015.