Translator Disclaimer
19 October 2017 Tissue multifractality and hidden Markov model based integrated framework for optimum precancer detection
Author Affiliations +
Abstract
We report the application of a hidden Markov model (HMM) on multifractal tissue optical properties derived via the Born approximation-based inverse light scattering method for effective discrimination of precancerous human cervical tissue sites from the normal ones. Two global fractal parameters, generalized Hurst exponent and the corresponding singularity spectrum width, computed by multifractal detrended fluctuation analysis (MFDFA), are used here as potential biomarkers. We develop a methodology that makes use of these multifractal parameters by integrating with different statistical classifiers like the HMM and support vector machine (SVM). It is shown that the MFDFA-HMM integrated model achieves significantly better discrimination between normal and different grades of cancer as compared to the MFDFA-SVM integrated model.

1.

Introduction

Disease diagnosis through optical methods is an area of considerable research interest.1 Optical tools are sensitive and are hence potentially capable of discriminating different stages of disease progression.17 However, tissue being a complex medium, with several fluorophores, scatterers, and absorption domains, makes it difficult for proper diagnosis through optical means.1 Hence, identifying reliable markers for accurately depicting the tissue condition through noninvasive optical methods has received significant attention.1 For this purpose, recent approaches have focused on extracting intrinsic fluorescence2 and tissue multifractality, characterizing the morphological changes by multifractal detrended fluctuation analysis (MFDFA).3 Other approaches make use of principal component analysis4 for identification of underlying spectral correlation and other image processing tools like wavelets5 for pin pointing parameters that faithfully capture the disease progression. Clinical application of this approach, not only depends on these biomarkers but also crucially depends on the validation of the diagnosis outcome through a suitable diagnostic algorithm, which can accurately classify the measured spectra from an unknown tissue, using the stored database of spectra of tissues of known histopathologic classification. This will supplement and augment the histopathological approach, the current industry gold standard. Over the last few decades, a variety of diagnostic algorithms have been developed for optical diagnosis of cancer.7 Classification schemes like artificial neural network6 and support vector machine (SVM)7 have been found promising in binary classification, e.g., normal versus cancer.

However, in a clinical situation, it is often required to classify the tissue site as normal and different grades of cancer. Hence, there is strong interest on statistical classifiers to extract the information content of the entire spectral data in order to get the best diagnostic features and enhance accuracy in classification of tissues into corresponding histopathologic categories.8 Total principal component regression developed by Tan et al.8 has classified various cancers, based on gene expression profiles and provided the optimized results compared to other methods for multiclass classifications. SVM has been deployed for multiclass cancer classification with classwise optimized gene.9 Its usefulness for optical diagnosis is yet to be explored. In recent times, hidden Markov models (HMMs) are being widely used in biological sequence analysis as a robust method.10 An HMM has been deployed to analyze hyperspectral images and a new HMM-based spectral measure has been referred to as the HMM information divergence in order to characterize the spectral properties.11

Here, we demonstrate the efficacy of MFDFA-HMM integrated framework for optical diagnosis of cancer. More specifically, it is found that HMM on the multifractal light scattering properties of the tissues shows remarkable efficiency in differentiating normal and different stages of precancer. It is particularly effective when applied on global fractal parameters like generalized Hurst exponent and the corresponding singularity spectrum width/strength of multifractality, characterizing the global morphological conditions of the tissues for multiclass cancer classification, as compared to the MFDFA-SVM integrated framework under same application.

2.

Methods and Materials

2.1.

Sample Preparation

Biopsied cervical precancer tissue slices (the histopathologically characterized grade I, grade II, and grade III precancer tissue) and normal tissues were collected from Ganesh Shankar Vidyarthi Memorial (GSVM) Medical College, Kanpur, India (age of patients between 35 and 60 years; ntotal=35, with ngradeI=14, ngradeII=6, ngradeIII=9; and six biopsies from the normal counterparts, nnormal=6). The standardized histological preparation of the excised tissues involving fixation, dehydration, embedding in wax, sectioning under a rotary microtome with thickness 5  μm, and lateral dimension 4  mm×6  mm, is followed by performing subsequent dewaxing. The consent for the use of all the intact tissue (human cervix with cancer and normal) samples in our study was obtained from the Ethical Committee, GSVM Medical College and Hospital, Kanpur, India. The sample preparation methods follow approved guidelines in our study.

2.2.

Experimental System

The spatial distribution of tissue refractive index (RI) was recorded by a differential interference contrast (DIC) microscope (Olympus IX-81, United States). At a magnification of 60×, these DIC images were recorded by a CCD camera (ORCA-ERG, Hamamatsu, 1344×1024  pixel dimension 6.45  μm). The elastic scattering spectra from the multiple sites of the biopsied tissue sections were recorded by the angle resolved spectral light scattering measurements (Fig. 1). In brief, light emitted from a Xe-lamp (HPX-2000, Ocean Optics, United States) was collimated by a combination of lenses and illuminated the tissue sample at the center of a goniometric arrangement (spot size 1-mm-diameter). The collimated scattered light from the sample was focused into a collecting fiber probe coupled to a spectrometer (USB4000FL, Ocean Optics) for wavelength resolved signal detection. The recordings of spectra were performed (360 to 800 nm) with a spectral resolution of 2.05 nm, where the angular range was kept at 10 deg to 150 deg with an interval of 10 deg. For the inverse multifractal study, the spectra were recorded at backscattering angle θ=150  deg (Fig. 3).

Fig. 1

Schematic of the spectral light scattering measurement. Xe lamp: excitation source; A: aperture; L1: collimating lens; L2: illuminating lens; L3 and L4: collecting lenses; f1: focal length of collimating lens L1; f2: focal length of illuminating lens L2; and f3 and f4: focal lengths of collecting lenses L3 and L4, respectively.

JBO_22_10_105005_f001.png

2.3.

Background Analysis

In the first step, the elastic scattering spectra data were processed through Fourier domain preprocessing via the Born approximation, followed by the multifractal analysis. For analyzing the dataset, 6 normal, 14 grade-I, 6 grade-II, and 9 grade-III samples were taken.

2.3.1.

Light scattering-based inverse analysis in Born approximation

At first, we extract the fractal parameters from the spatial variation of RI, as manifested through corresponding light scattering data.12 We have followed the inverse analysis strategy12 for quantification of multifractal signature from the scattered light signal through MFDFA, before applying an HMM on them. The normalized RI fluctuations of a weakly fluctuating scattering medium are given by Δn(r)n(r)n0n0, where n0 and r are the average RI of the medium and location within the volume, respectively. The fluctuation part Δn(r) is responsible for phase distortion and scattering. It is known that the elastic scattering field for scalar excitation can be related to Δn(r) in first-order Born approximation via Fourier transform.12 For continuous random media, such as tissues, the expression for scattered intensity is given by

Eq. (1)

I(β)k4σ2|η(r)ei(βr)d3r|2,
where λ is the wavelength; k=2π/λ; θ is the scattering angle; η(r) is the spatial inhomogeneity distribution of index; and σ=n0δn is the strength of fluctuation index, where the amplitude of fluctuating index is δn. The scattering vector β has the modulus β=2ksin(θ/2), where β is related to the spatial frequency ν via β=2πν. For a medium exhibiting self-similarity in index fluctuations, the information on multifractality can be extracted from the scattering signal as follows:

Eq. (2)

η/(ρ)k2I(β=2π)ei(βr)d3β.

The parameter η/(ρ) exhibits the index inhomogeneity distribution with spatial scale ρ=|rr/|, where ρ is the distance between any two points in the medium. It represents the randomness of the medium in statistical sense and embodies the essential multifractal features of index fluctuations in complex systems, such as tissues. This is subsequently analyzed using MFDFA13 to yield the multifractal tissue optical properties.

2.3.2.

Multifractal detrended fluctuation analysis

MFDA is an effective tool for quantitative estimation of the generalized Hurst exponents. In MFDFA, one first generates a profile function Y(i) (spatial series of length N, i=1,,N) from the one-dimensional spatial index fluctuations. Here, length of the series N=256×256  pixels. Subsequently, the profile is divided into Ns=int(N/s) nonoverlapping segments b of equal segment length s. The segment length s varies from 16 to 128. The local trend of the series [yb(i)] is determined for each segment b by least square polynomial fitting and then subtracted from the segmented profiles to yield detrended fluctuations. The resulting variance of the detrended fluctuation is determined for each segment as follows:

Eq. (3)

F2(b,s)=1si=1s{Y[(b1)s+i]yb(i)}2.

The moment (q) dependent fluctuation function can be obtained by taking the average over all the segments:

Eq. (4)

Fq(s)={12Nsb=12Ns[F2(b,s)]q2}1q,
where q may vary from negative to positive values with fraction or integer values. Since length N of the series is often not an integer multiple of variable segment length s, in order to take into account all the data points, the procedure was performed twice on the series, starting from either end of the series.13 The scaling behavior can be obtained by analyzing the variations of   Fq(s) versus s for each value q, assuming the general scaling function as follows:

Eq. (5)

Fq(s)sh(q).

The relation between the generalized Hurst exponent h(q) and the multifractal scaling exponent τ(q) can be demonstrated as follows:

Eq. (6)

τ(q)=qh(q)1.

The Hurst exponent (H) is defined as h(q=2) and the values H>0.5, H=0.5, and H<0.5 correspond to long range correlation, uncorrelated random fluctuations, and anticorrelated behavior, respectively.11 The Hurst exponent h(q) and the scaling exponent τ(q), along with the singularity spectrum f(α) completely, characterize any nonstationary multifractal fluctuation series. Here, f(α) is related to τ(q) via a Legendre transformation:

Eq. (7)

α=dτdq,f(α)=qατ(q),
where α is the singularity strength and the width (Δα=|α1α2|) [considered at f(α)=0] of f(α) is a quantitative measure of multifractality.13

2.4.

Data Analysis

The obtained two global fractal parameters, generalized Hurst exponent and the corresponding singularity spectrum width, were subjected to (a) SVM-based multiclass classification and (b) HMM-based multiclass classification. For analyzing the dataset, 6 normal, 14 grade-I, 6 grade-II, and 9 grade-III samples were taken.

2.4.1.

Support vector machine

SVMs are powerful statistical classifiers under the supervised learning scheme. The central idea behind SVM operation is to separate classes with a surface that maximizes margin between them by avoiding overfitting to form an optimal separating hyperplane (OSH). Hence, by following structural risk minimization (SRM) of statistical learning makes prediction on a function f(x) as f(x)=i=1Nwik(x,xi)+w0, where k(x,xi) is the kernel function defined on a basis function, {wi} is the corresponding model weights, and w0 is the bias weight.

The training data points lie far away from the OSH, does not participate in the specification and hence receives zero weight. Data point that lies close to decision boundary receives nonzero weights. These training data points are “support vectors.”14,15 If we remove these points, it will change the boundary location. Unlike relevance vector machine, there are restrictions while choosing of kernels in SVM. An appropriate selection of kernel function is an important aspect as it defines the accuracy level of SVM-based operation while determining training data classification. The kernel function will produce optimum results in classification as long as it obeys the Mercer’s theorem.14,15 Figure 2 displays the simplified workflow of SVM-based multiclass classification on extracted multifractal parameters from tissue samples.

Fig. 2

SVM-based multiclass classification on extracted multifractal parameters from tissue samples.

JBO_22_10_105005_f002.png

2.4.2.

Hidden Markov model

HMM16 is a statistical Markov model with hidden states and is also the simplest dynamic Bayesian network. An HMM is closely related with mixture models, which are statistically independent. An HMM can be efficiently employed for a time series data, where actual parameters are unknown and only observational information are known. From this series of observations, the probabilities of parameters giving such observations and the transmission probabilities can be found by the Baum–Welch algorithm.17 The basic principle of an HMM can be described as follows. For s0,,st states as input to the model, the st+1’th state can be predicted by the traditional Markov model, where given the present input, the future is independent of the past:

Eq. (8)

p(st+1|st,,s0)=p(st+1|st).

If qn’th observation can be made on the basis of qn1’th observation, it is a first-order assumption, which generally is used for Bayesian modeling. Using the Markov assumption, we also can write it as follows:

Eq. (9)

p(q1,,qn)=i=1np(qi|qi1).

According to the Bayes’ formula:

Eq. (10)

p(qi/xi)=p(xi/qi)p(qi)p(xi).

In a more general way, the above Eq. (8) can be rewritten as follows:

Eq. (11)

p(q1,,qnx1,,xn)=p(x1,,xnq1,,qn)p(q1,,qn)p(x1,,xn).

The measure of probability can be achieved by likelihood parameter L, which is proportional to the probability:

Eq. (12)

p(q1,,qnx1,,xn)L(q1,,qnx1,,xn)=p(x1,,xnq1,,qn)p(q1,,qn).

In the case of an HMM, the model is prepared with the training data Θ={π,A,B} and a sequence of N states is S={s1,,sN}. Here, π is defined as the prior probability, A is the transition probability, and B is the emission probability. The probability of a state sequence Q={q1,,qn} obtained from an HMM with parameters Θ can be expressed as a product of the transition probabilities:

Eq. (13)

p(QΘ)=πq1n=1N1aqn,qn+1=πq1aq1,q2aq2,q3aqN1,qN.

For an observational sequence X={x1,,xN}, a (hidden) state sequence Q={q1,,qN} can be determined from an HMM with parameter Θ. Hence, the likelihood of X along the path Q takes the form:

Eq. (14)

p(XQ,Θ)=n=1Np(xnqn,Θ)=bq1,x1.bq2,x2bqN,xN.

It can be expressed as the product of the emission probabilities computed along the considered path.

With the likelihood of an observed sequence X={x1,,xN} and the parameter Θ defined by an HMM, the probability p can be expanded as follows:

Eq. (15)

p(XΘ)=p(X,QΘ).

Using the Baum–Welch algorithm, the hidden parameters in an HMM are found. This algorithm utilizes the expectation maximization (EM) algorithm for finding the maximum likelihood estimation of the parameters of an HMM, given a set of observed feature vectors.

As we mentioned earlier, hidden Markov chain can be represented as θ=(A,B,π). Here, the stochastic transition matrix A={aij}=p(Xt=j|Xt1=i), where Xt is the discrete variable. The emission probability B=[bj(yt)]=p(Yt=yt|Xt=j), where Yt is the observational sequence. The EM algorithm defines a local maximum for θ*=argmax[p(Y/X)]. After defining the initial condition, the Baum–Welch algorithm follows the forward and backward procedure to find the proper estimation of the predicted results.

In the case of forward procedure, let the probability of viewing y1,,yt at state i in time t is αi(t)=p(Y1=y,1,Yt=yt|Xt=i,θ).

Under recursion procedure, αi(1)=πibi(y1) and αj(t+1)=bj(yt+1)i=1Nαi(t)aij.

Similarly, let the probability of viewing yt+1,,yT at state i in time t is βi(t)=p(Yt+1=yt+1,,YT=yT|Xt=i,θ). Under recursion procedure, βi(T)=l and βi(t)=j=1Nβj(t+1)aijbj(yt+1).

In the final step, according to the Bayes’ theorem, the probability of the observed sequence Y and the parameters θ in state i at time t given as γi(t)=p(Xt=i|Y,θ)=αi(t)βi(t)j=1Nαj(t)βj(t).

The probability of being in state i and j at times t and t+1, respectively, given the observed sequence Y and parameters θ:

ϵij(t)=p(Xt=i,Xt+1=j|Y,θ)=αi(t)aijβj(t+1)bj(yt+1)k=1Nαk(T).

Hence, θ can be updated at expected frequency spent in state i at time l as πi*=γi(l). The expected number of transitions from state i to state j compared to the expected total number of transitions away from state i is aij*=t=1T1ϵij(t)t=1T1γi(t).

If bi*(vk) is the expected number of times, the output observations have been equal to vk while in state i over the expected total number of times in state i, then bi*(vk)=t=1Tlyt=vkγi(t)t=1Tγi(t). Here, the indicator function lyt=vk=1 exists only for yt=vk.

These steps are iterated until a desired level of convergence is achieved.

In this paper, an HMM is applied on the known multifractal fractal parameters, which leads to a significant improvement on the prediction accuracy. For experimentation, we first train the multifractal parameters in an HMM for each of the categories. Prior probabilities are first selected as a random function. A and B are modeled as Gaussian densities, with mean 0 and variance 1. Subsequently, a representative data is trained on the model iteratively to fit and modify the model using EM algorithm. The model is optimized using Lagrange multipliers. We use forward and backward algorithms to compute a set of sufficient statistics for our EM step tractably. Once the model is sufficiently trained for a given sequence of data we calculate the likelihood of sequence for each category, i.e., as P(X/θi), when the sum of the joint likelihoods of the sequence over all possible state sequences Q, allowed by the model for each category. The maximum likelihood gives the prediction for the sequence data.

The proposed HMM-based model on light scattering derived from multifractal tissue optical properties has been demonstrated in Fig. 3. The HMM-based data analysis steps for normal and different cancerous grades have been shown here in detail for the ease of understanding.

Fig. 3

The flowchart of an HMM-based model on multifractal tissue optical properties derived from light scattering spectra. First row depicts the singularity spectrum versus generalized Hurst exponent for testing and training datasets. Second row shows the processing of test data and training data through emission probability density, initial state probability, and transition probability (model trained using the training data) has been performed on the normal and different grades of cancer, respectively. The output has been processed by the maximum likelihood classifier method to perform the HMM-based classification.

JBO_22_10_105005_f003.png

3.

Results and Discussion

The DIC images of different pathological grades have been presented in Fig. 4 for comparisons. The histopathologically characterized tissue samples were provided by the pathologists of GSVM Medical College and Hospital, Kanpur, where cancer stages were defined by the pathologists.

Fig. 4

The representative DIC images of (a) healthy (normal) and (b–d) different stages (histopathologically characterized grade I, II, III) of precancer tissues, respectively.

JBO_22_10_105005_f004.png

The results of the inverse analysis on the light scattering spectra recorded (using the spectral light scattering measurement system in Fig. 1) from a grade-1 dysplastic cervical tissue (corresponding to Fig. 1) are displayed in Fig. 5. The large variation of the slopes of log Fq(s) versus logs [Fig. 5(c)] is the evidence of strong multifractality. The derived h(q) spectrum and the singularity spectrum f(α) [Figs. 5(d) and 5(e)] demonstrate the strong multifractality in the spatial variations of tissue refractive indices. The parameter η/(ρ), obtained by Fourier preprocessing, contains submicron level spatial index fluctuations information [as evident from Fig. 5(b)]. The observed small-scale fluctuations of h(q) [Fig. 5(d)] and the resulting width of f(α) [Fig. 5(e)] possibly originate from the microarchitecture of the fibrous network of connective tissue. To summarize Fourier preprocessing of light scattering in Born approximation and its subsequent analysis through MFDFA documents, there are small spectral variations as signature of subtle or hidden changes in the refractive indices spatial distribution via multifractal parameters.

Fig. 5

Results of the inverse analysis performed on the light scattering spectra recorded from a grade-1 precancer human cervical tissue slice (corresponding to Fig. 1). (a) The recorded light scattering intensity [I(v)versusv,v=2λsin(θ2), θ=150  deg, λ=360 to 740 nm, shown here]. (b) The representative index fluctuations with spatial scale ρ [η/(ρ)] extracted via Fourier domain preprocessing [using Eq. (2)] on I(v) (shown following polynomial detrending). (c)–(e) Results of the MFDFA inverse analysis on the detrended fluctuations η/(ρ). (c) The variation of logFq(s) versus logs for different moments q (=4 to +4 shown here). Multifractality in η/(ρ) is evident from significant variations in the slopes with varying q. (d) The MFDFA-derived moment dependence of generalized Hurst exponent h(q). (e) The resulting singularity spectrum f(α).

JBO_22_10_105005_f005.png

After MFDFA analysis, the observed trends are h(q=2)=0.63±0.02, 0.56±0.05, 0.48±0.03, 0.41±0.04, and Δα=0.86±0.01, 0.90±0.03, 0.96±0.04, 0.99±0.01 for normal and different grades of cancer, respectively. Here, it can be observed that there exists an overlapping of multifractal parameters like Hurst exponent and singularity spectrum width among normal and different grade of precancerous tissues. Hence, the supervised classifiers like SVM and HMM have been applied here for multiclass classification purpose.

Our initial dataset consists of 35 samples. We have used Monte Carlo cross-validation, where we randomly split the dataset into training and testing dataset. This process has been repeated 100 times. The size of training and testing dataset varies for each split; we just ensured that a minimum two number of samples have been included in both training and testing dataset for each split. For each such split, the model has been fit to the training data, and predictive accuracy has been assessed using just the validation data (testing data). Then, the results of the testing data have been averaged over the splits. The advantage of this method is that the proportion of the training/validation split is not dependent on the number of iterations, and the predictive accuracy is more or less independent of the samples used for training dataset. We additionally used nine unknown samples taken at different time than the dataset and the prediction done by our model has been compared by manual verification. Results of the unknown samples also have been averaged and presented in the results. The above process has been repeated for SVM and HMM.

The data in Table 1 depicts the mean and variance (mean±variance) for each MFDA parameter of each grade for the entire dataset. The same dataset has been used for model generation in SVM and HMM.

Table 1

Summary of the multifractal tissue optical properties derived from light scattering spectra.

NormalGrade-IGrade-IIGrade-III
Hurst exponent [h(q=2)]0.63±0.020.56±0.050.48±0.030.41±0.04
Singularity spectrum width (Δα)0.86±0.010.90±0.030.96±0.040.99±0.01

There are nine unknown samples, which are considered as test samples for SVM and HMM classification purpose.

SVM creates an optimum manifold barrier with radial basis function kernel between healthy and different grades of cancer depending on MFDFA parameters. Figure 6 displays the prediction analysis carried out by training data in SVM.

Fig. 6

SVM-based tissue classification based on light scattering-derived multifractal tissue optical properties. The horizontal surfaces of these figures have been sectioned into 4×4 rectangles. The accurate and inaccurate prediction of each stage (normal, grade I, grade II, and grade III) have been represented by diagonal and off diagonal rectangles, respectively.

JBO_22_10_105005_f006.png

From Fig. 6, it is clearly visible that SVM works very well in binary classification, i.e., between normal and grade III with 98.5% and 100% accuracy, respectively. While trying out SVM classification for multiclass classification, the MFDFA-SVM integrated framework performs poorly by degrading the overall performance of the system with 57.14% and 55% accuracy, respectively. This error in prediction has occurred due to wrongly predicting normal tissues as grade I tissues (1.5%), grade I as normal (14.29%), and grade II tissues (28.57%) as well as grade II tissues as normal (10%) and grade I tissues (35%).

An HMM creates abstract Markov models for the classification purpose on multifractal parameters. Figure 7 displays the prediction analysis carried out by training data in an HMM.

Fig. 7

HMM-based tissue classification based on light scattering –derived multifractal tissue optical properties. The horizontal surfaces of these figures have been sectioned into 4×4 rectangles. The accurate and inaccurate prediction of each stage (normal, grade I, grade II, and grade III) have been represented by diagonal and off diagonal rectangles, respectively.

JBO_22_10_105005_f007.png

As can be seen from the graph, the parameters clearly show distinction between normal and different grades of precancer with singularity spectrum width (Δα) and Hurst exponent h(q). The normal and grade III tissues were correctly predicted. However, there is error while predicting grade I and II tissues (correct prediction rates are 85.45% and 78.57%, respectively). This error in prediction is also limited to wrongly predicting grade I tissues as grade II tissues (14.54%) as well as grade II tissues as grade I tissues (21.43%).

The results demonstrate that binary classification between normal and cancerous tissues (grade III) is very good both in SVM and HMM. Meanwhile, in multiclass classification cases, when precancerous grades (grade I, grade II) are to be classified along with normal and precancerous tissues (grade III), abstract parameters achieved using an HMM that performs better than the SVM. The presence of noise in the obtained signal damages the SVM performance as SVM clearly classifies based on the kernel formed after considering all the multifractal parameters. While in the case of HMM, the Markov model finds abstract parameters by controlling the actual multifractal parameters and produces a prediction based on the derived abstract parameters. As a consequence, an HMM avoids the noise added to the signal and able to produce better multiclass classification results than SVM.

In our current manuscript, we proceeded with Fourier domain preprocessing as our main focus was on global behavior analysis of time series through MFDFA-HMM and MFDFA-SVM integrated model classifications in early-stage cancer detection. It is known to us that wavelet domain analysis is more appreciable than Fourier domain (short-time Fourier transform) analysis in the case of local behavior analysis of time series.18 In wavelet domain preprocessing, the performance of the SVM model can perform better than an HMM model as SVM has better generalization due to the principle of SRM than an HMM for system abnormality detection.19 A comparative study between MFDFA-SVM and MFDFA-HMM integrated model for local behavior analysis of time series through wavelet domain preprocessing in early-stage cancer diagnosis will be a part of our future study.

4.

Conclusions

We have explored an integrated framework of light scattering-derived multifractal tissue optical properties (generalized Hurst exponent and width of singularity spectrum) along with a robust HMM for multiclass classification of different precancerous grades of human uterine cervix. The results clearly demonstrate that the use of HMM on the multifractal properties leads to significantly improved classification as compared to MFDFA-SVM based integrated model for multiclass classification. These MFDFA-HMM based classification results show considerable promise by exploring multifractal tissue optical properties as a biomarker for precancer detection. We are currently expanding our investigations toward in-vivo deployment of this integrated approach for precancer detection using tissue light scattering spectra. In general, the use of this MFDFA-HMM integrated model on elastic scattering spectroscopic data may lead to a diagnostic modality for the detection of other types of cancer.

Disclosures

No conflicts of interest, financial or otherwise, are declared by the authors.

Acknowledgments

The authors thank Dr. Asha Agarwal, GSVM Medical College and Hospital, Kanpur, for providing the histopathologically characterized tissue samples.

References

1. 

P. N. Prasad, Introduction to Biophotonics, Wiley-Interscience, Hoboken, New Jersey (2003). Google Scholar

2. 

K. Pandey et al., “Fluorescence spectroscopy: a new approach in cervical cancer,” J. Obstet. Gynaecol. India, 62 (4), 432 –436 (2012). https://doi.org/http://dx.doi.org/10.1007/s13224-012-0298-6 JOBYA4 Google Scholar

3. 

N. Das et al., “Probing multifractality in tissue refractive index: prospects for precancer detection,” Opt. Lett., 38 211 –213 (2013). https://doi.org/http://dx.doi.org/10.1364/OL.38.000211 OPLEDP 0146-9592 Google Scholar

4. 

S. Devi, P. K. Panigrahi and A. Pradhan, “Detecting cervical cancer progression through extracted intrinsic fluorescence and principal component analysis,” J. Biomed. Opt., 19 127003 (2014). https://doi.org/http://dx.doi.org/10.1117/1.JBO.19.12.127003 JBOPFO 1083-3668 Google Scholar

5. 

A. H. Gharekhan et al., “Distinguishing cancer and normal breast tissue autofluorescence using continuous wavelet transform,” IEEE J. Sel. Top. Quantum Electron., 16 893 –899 (2010). https://doi.org/http://dx.doi.org/10.1109/JSTQE.2009.2033018 IJSQEN 1077-260X Google Scholar

6. 

G. A. Rovithakis et al., “Artificial neural networks for discriminating pathologic from normal peripheral vascular tissue,” IEEE Trans. Biomed. Eng., 48 (10), 1088 –1097 (2001). https://doi.org/http://dx.doi.org/10.1109/10.951511 Google Scholar

7. 

S. K. Majumder, N. Ghosh and P. K. Gupta, “Support vector machine for optical diagnosis of cancer,” J. Biomed. Opt., 10 (2), 024034 (2005). https://doi.org/http://dx.doi.org/10.1117/1.1897396 JBOPFO 1083-3668 Google Scholar

8. 

Y. Tan et al., “Multi-class cancer classification by total principal component regression (TPCR) using microarray gene expression data,” Nucleic Acids Res., 33 (1), 56 –65 (2005). https://doi.org/http://dx.doi.org/10.1093/nar/gki144 Google Scholar

9. 

A. Anand and P. N. Suganthan, “Multiclass cancer classification by support vector machines with class-wise optimized genes and probability estimates,” J. Theor. Biol., 259 (3), 533 –540 (2009). https://doi.org/http://dx.doi.org/10.1016/j.jtbi.2009.04.013 JTBIAP 0022-5193 Google Scholar

10. 

B. Yoon, “Hidden Markov models and their applications in biological sequence analysis,” Curr. Genomics, 10 402 –415 (2009). https://doi.org/http://dx.doi.org/10.2174/138920209789177575 CGUEA8 1389-2029 Google Scholar

11. 

Q. Du and C. I. Chang, “Hidden Markov model approach to spectral analysis for hyperspectral imagery,” Opt. Eng., 40 (10), 2277 –2284 (2001). https://doi.org/http://dx.doi.org/10.1117/1.1404430 Google Scholar

12. 

N. Das et al., “Tissue multifractality and Born approximation in analysis of light scattering: a novel approach for precancers detection,” Sci. Rep., 4 6129 (2014). https://doi.org/http://dx.doi.org/10.1038/srep06129 SRCEC3 2045-2322 Google Scholar

13. 

J. W. Kantelhardt et al., “Multifractal detrended fluctuation analysis of nonstationary time series,” Phys. A, 316 87 –114 (2002). https://doi.org/http://dx.doi.org/10.1016/S0378-4371(02)01383-3 Google Scholar

14. 

V. N. Vapnik, Statistical Learning Theory, Wiley, New York (1998). Google Scholar

15. 

C. J. C. Burges, “A tutorial on support vector machines for pattern recognition,” Data Min. Knowl. Discovery, 2 (2), 121 –167 (1998). https://doi.org/http://dx.doi.org/10.1023/A:1009715923555 Google Scholar

16. 

L. E. Baum and T. Perie, “Statistical inference for probabilistic functions of finite state Markov chains,” Ann. Math. Stat., 37 (6), 1554 –1563 (1966). https://doi.org/http://dx.doi.org/10.1214/aoms/1177699147 AASTAD 0003-4851 Google Scholar

17. 

R. Hogg, J. McKean and A. Craig, Introduction to Mathematical Statistics, 7th ed.Pearson Prentice Hall, Boston (2005). Google Scholar

18. 

F. Auger and P. Flandrin, “Improving the readability of time–frequency and time-scale representations by the reassignment method,” IEEE Trans. Signal Process., 43 (5), 1068 –1089 (1995). https://doi.org/http://dx.doi.org/10.1109/78.382394 ITPRED 1053-587X Google Scholar

19. 

Q. Miao, H. Huang and X. Fang, “A comparison study of support vector machines and hidden Markov models in machinery condition monitoring,” J. Mech. Sci. Technol., 21 607 –615 (2007). https://doi.org/http://dx.doi.org/10.1007/BF03026965 1738-494X Google Scholar

Biographies for the authors are not available.

© 2017 Society of Photo-Optical Instrumentation Engineers (SPIE) 1083-3668/2017/$25.00 © 2017 SPIE
Sabyasachi Mukhopadhyay, Nandan K. Das, Indrajit Kurmi, Nirmalya Ghosh, Asima Pradhan, and Prasanta K. Panigrahi "Tissue multifractality and hidden Markov model based integrated framework for optimum precancer detection," Journal of Biomedical Optics 22(10), 105005 (19 October 2017). https://doi.org/10.1117/1.JBO.22.10.105005
Received: 12 May 2017; Accepted: 29 September 2017; Published: 19 October 2017
JOURNAL ARTICLE
8 PAGES


SHARE
Advertisement
Advertisement
Back to Top