Dynamic Functional Connectivity Analysis Using Network-Based Brain State Identification, Application on Temporal Lobe Epilepsy.

Epilepsy is a brain network disorder caused by discharges of interconnected groups of neurons and resulting brain dysfunction. The brain network can be characterized by intra- and inter-regional functional connectivity (FC). However, since the BOLD signal is inherently non-stationary, the FC is evidenced to be varying over time. Considering the dynamic characteristics of the functional network, we aimed to obtain dynamic brain states and their properties using network-based analyses for the comparison of healthy control and temporal lobe epilepsy (TLE) groups and also lateralization of TLE patients. We used dwelling time, transition time, and brain network connection in each state as the dynamic features for this purpose. Results showed a significant difference in dwelling time and transition time between the healthy control group and both left TLE and right TLE groups and also a significant difference in brain network connections between the left and right TLE groups.

Comparison of multimodal findings on epileptogenic side in temporal lobe epilepsy using self-organizing maps.

OBJECTIVE: To develop a decision-making tool to evaluate and compare the performance of neuroimaging markers with clinical findings and the significance of attributes for presurgical lateralization of mesial temporal lobe epilepsy (mTLE). METHODS: Thirty-five unilateral mTLE patients who qualified as candidates for surgical resection were studied. Seizure semiology, ictal EEG, ictal epileptogenic zone, interictal-irritative zone, and MRI findings were used as clinical markers. Hippocampal T1 volumetry and FLAIR intensity, DTI estimated; mean diffusivity (MD) in the hippocampus and fractional anisotropy (FA) in posteroinferior cingulum and crus of fornix, and the output of logistic regression method on volumetrics of the hippocampus, amygdala, and thalamus were adopted as neuroimaging markers. The self-organizing map (SOM) method was applied to markers to provide predictive methods for mTLE lateralization. RESULTS: The SOM clustered all clinical attributes correctly with 100% accuracy and sensitivity for both the left and right mTLE. Among the clinical markers, seizure semiology and interictal-irrelative zone are the most sensitive attribute for the left-mTLE group lateralization. The accuracy achieved by applying the SOM method to the neuroimaging attributes was 94%, while the sensitivity was achieved 90% for left and 100% for right mTLE. SOM evidence indicated that the hippocampal volume is the most sensitive attribute for the prediction of the laterality in left-mTLE groups. CONCLUSION: The proposed SOM method showed that neuroimaging markers may not replace with clinical findings. Nevertheless, multimodal neuroimaging can play an effective role in preoperative lateralization to reduce the costs and risks of surgical resection.

Identification of effective features of LFP signal for making closed-loop deep brain stimulation in parkinsonian rats.

Traditional deep brain stimulation (DBS) is one of the acceptable methods to relieve the clinical symptoms of Parkinson's disease in its advanced stages. Today, the use of closed-loop DBS to increase stimulation efficiency and patient satisfaction is one of the most important issues under investigation. The present study was aimed to find local field potential (LFP) features of parkinsonian rats, which can determine the timing of stimulation with high accuracy. The LFP signals from rats were recorded in three groups of parkinsonian rat models receiving stimulation (stimulation), without getting stimulation (off-stimulation), and sham-controlled group. The frequency domain and chaotic features of signals were extracted for classifying three classes by support vector machine (SVM) and neural networks. The best combination of features was selected using the genetic algorithm (GA). Finally, the effective features were introduced to determine the on/off stimulation time, and the optimal stimulation parameters were identified. It was found that a combination of frequency domain and chaotic features with an accuracy of about 99% was able to determine the time the DBS must switch on. In about 80.67% of the 1861 different stimulation parameters, the brain was able to maintain its state for about 3 min after stimulation discontinuation.

Dynamic functional connectivity in temporal lobe epilepsy: a graph theoretical and machine learning approach.

PURPOSE: Functional magnetic resonance imaging (fMRI) in resting state can be used to evaluate the functional organization of the human brain in the absence of any task or stimulus. The functional connectivity (FC) has non-stationary nature and consented to be varying over time. By considering the dynamic characteristics of the FC and using graph theoretical analysis and a machine learning approach, we aim to identify the laterality in cases of temporal lobe epilepsy (TLE). METHODS: Six global graph measures are extracted from static and dynamic functional connectivity matrices using fMRI data of 35 unilateral TLE subjects. Alterations in the time trend of the graph measures are quantified. The random forest (RF) method is used for the determination of feature importance and selection of dynamic graph features including mean, variance, skewness, kurtosis, and Shannon entropy. The selected features are used in the support vector machine (SVM) classifier to identify the left and right epileptogenic sides in patients with TLE. RESULTS: Our results for the performance of SVM demonstrate that the utility of dynamic features improves the classification outcome in terms of accuracy (88.5% for dynamic features compared with 82% for static features). Selecting the best dynamic features also elevates the accuracy to 91.5%. CONCLUSION: Accounting for the non-stationary characteristics of functional connectivity, dynamic connectivity analysis of graph measures along with machine learning approach can identify the temporal trend of some specific network features. These network features may be used as potential imaging markers in determining the epileptogenic hemisphere in patients with TLE.

Toward a closed-loop deep brain stimulation in Parkinson's disease using local field potential in parkinsonian rat model.

Deep brain stimulation (DBS) is an invasive method used for treating Parkinson's disease in its advanced stages. Nowadays, the initial adjustment of DBS parameters and their automatic matching proportion to the progression of the disease is viewed as one of the research areas discussed by the researchers, which is called closed-loop DBS. Various studies were conducted regarding finding the signal(s) which reflects different symptoms of the disease. Local Field Potential (LFP) is one of the signals that is suitable for using as feedback, because it can be recorded by the same implemented electrodes for stimulation. The present study aimed to identify the distinguishing features of patients from healthy individuals using LFP signals. METHODS: In the present study, LFP was recorded from the rats in sham and parkinsonian model groups. After evaluating the signals in the frequency domain, sixty-six features were extracted from power spectral density of LFPs. The features were classified by Support Vector Machine (SVM) to determine the ability of features for separating parkinsonian rats from healthy ones. Finally, the most effective features were selected for distinguishing between the sham and parkinsonian model groups using a genetic algorithm. RESULTS: The results indicated that the frequency domain features of LFP signals from rats have capacity of using them as a feedback for closed-loop DBS. The accuracy of the Support Vector Machine classification using all 66 features was 80.42% which increased to 84.41% using 38 features selected by genetic algorithm. The proposed method not only increase the accuracy, but it also reduce computation by decreasing the number of the effective features. The results indicate the significant capacity of the proposed method for identifying the effective high-frequency features to control the closed-loop DBS. CONCLUSIONS: The ability of using LFP signals as feedback in closed-loop DBS was shown by extracting useful information in frequency bands below and above 100 Hz regarding LFP signals of parkinsonian rats and sham ones. Based on the results, features at frequencies above 100 Hz were more powerful and robust than below 100 Hz. The genetic algorithm was used for optimizing the classification problem.

Evolution of Graph Theory in Dynamic Functional Connectivity for Lateralization of Temporal Lobe Epilepsy.

Resting-state functional magnetic resonance imaging (rsfMRI) has described the functional architecture of the human brain in the absence of any task or stimulus. Since the functional connectivity (FC), has non-stationary nature, it is evidenced to be varying over time. Using dynamic functional connectivity, six graph theoretical characteristics were measured and compared between left and right temporal lobe epilepsy (TLE). We also obtain a trend for each characteristic in the time course of experiments. The results demonstrated that the static connectivity analysis failed to fully separate the left and right TLE patients for some characteristics, whereby the dynamic analysis has been shown capable of identifying the laterality. Furthermore, the results suggest that the temporal trend of some graph theoretical characteristics can be exploited as a novel marker for TLE laterality.

SSVEP recognition by modeling brain activity using system identification based on Box-Jenkins model.

The steady-state visual evoked potential (SSVEP) based brain-computer interface (BCI) has received increasing attention in recent years. The present study proposes a new method for recognition based on system identification. The method relies on modeling the electroencephalogram (EEG) signals using the Box-Jenkins model. In this approach, the recorded EEG signal is considered as a combination of an SSVEP signal evoked by periodic visual stimulation and a background EEG signal whose components are modeled by a moving average (MA) process and an auto-regressive moving average (ARMA) process, respectively. Then, the target frequency is determined by comparing the modeled SSVEP signals for all stimulation frequencies. The experimental results of the proposed method for recorded EEG signals from five subjects (each subject with four stimulation frequencies) demonstrated a significant improvement in the accuracy of the SSVEP recognition in contrast to canonical correlation analysis, least absolute shrinkage and selection operator, and multivariate linear regression methods. The proposed method exhibits enhanced accuracy especially for short data length and a small number of channels. This superiority suggests that the proposed method is an appropriate choice for the implementation of real-time SSVEP based BCI systems.

A novel approach to the extraction of fetal electrocardiogram based on empirical mode decomposition and correlation analysis.

Fetal heart rate monitoring is the process of checking the condition of the fetus during pregnancy and it would allow doctors and nurses to detect early signs of trouble during labor and delivery. The fetal ECG (FECG) signal is so weak and also is corrupted by other signals and noises, mainly by maternal ECG signal. It is so hard to acquire a noise-free, precise and reliable FECG using the conventional methods. In this study, a combination of empirical mode decomposition (EMD) algorithms, correlation and match filtering is used for extracting FECG from maternal abdominal ECG signals. The proposed method benefits from match filtering ability to detect fetal signal and QRS complex to detect weak QRS peaks. The combined method, has been applied successfully on different signal qualities, even for signals that their analysis was hard and complicated for other methods. This method is able to detect R-R intervals with high accuracy. It was proved that the complete ensemble empirical mode decomposition method provides a better frequency resolution of modes and also requires less iterations that leads to a considerably less computational cost than EMD and ensemble empirical mode decomposition and can reconstruct the FECG completely from the calculated modes. We believe that this method opens a new field in non-invasive maternal abdominal signal processing so the FECG signal could be extracted with high speed and accuracy.

Providing an Efficient Algorithm for Finding R Peaks in ECG Signals and Detecting Ventricular Abnormalities With Morphological Features.

Ventricular arrhythmias are one of the most important causes of annual deaths in the world, which may lead to sudden cardiac deaths. Accurate and early diagnosis of ventricular arrhythmias in heart diseases is essential for preventing mortality in cardiac patients. Ventricular activity on the electrocardiogram (ECG) signal is in the interval from the beginning of QRS complex to T wave end. Variations in the ECG signal and its features may indicate heart condition of patients. The first step to extract features of ECG in time domain is finding R peaks. In this paper, a combination of two algorithms of Pan-Tompkins and state logic machine has been used to find R peaks in heart signals for normal sinus signals and ventricular abnormalities. Then, a healthy or sick beat may be realized by comparing the difference between R peaks obtained from two algorithms in each beat. The morphological features of the ECG signal in the range of QRS complex are evaluated. Ventricular tachycardia (VT), ventricular flutter (VFL), ventricular fibrillation (VFI), ventricular escape beat (VEB), and premature ventricular contractions (PVCs) are abnormalities studied in this paper. In the classification step, the support vector machine (SVM) classifier with Gaussian kernel (one in front of everyone) is used. Accuracy percentages of ventricular abnormalities mentioned above and normal sinus rhythm are respectively obtained as 95.8%, 92.8%, 94.5, 98.9%, 91.5%, and 100%. The database of this paper has been taken from normal sinus rhythm and MIT-SCD banks available on Physionet.org.

Arrangement and Applying of Movement Patterns in the Cerebellum Based on Semi-supervised Learning.

Biological control systems have long been studied as a possible inspiration for the construction of robotic controllers. The cerebellum is known to be involved in the production and learning of smooth, coordinated movements. Therefore, highly regular structure of the cerebellum has been in the core of attention in theoretical and computational modeling. However, most of these models reflect some special features of the cerebellum without regarding the whole motor command computational process. In this paper, we try to make a logical relation between the most significant models of the cerebellum and introduce a new learning strategy to arrange the movement patterns: cerebellar modular arrangement and applying of movement patterns based on semi-supervised learning (CMAPS). We assume here the cerebellum like a big archive of patterns that has an efficient organization to classify and recall them. The main idea is to achieve an optimal use of memory locations by more than just a supervised learning and classification algorithm. Surely, more experimental and physiological researches are needed to confirm our hypothesis.

Recognition of amyotrophic lateral sclerosis disease using factorial hidden Markov model.

Amyotrophic lateral sclerosis (ALS) is a common disease among neurological disorders that can change the pattern of gait in human. One of the effective methods for recognition and analysis of gait patterns in ALS patients is utilizing stride interval time series. With proper preprocessing for removing unwanted artifacts from the raw stride interval times and then extracting meaningful features from these data, the factorial hidden Markov model (FHMM) was used to distinguish ALS patients from healthy subjects. The results of classification accuracy evaluated using the leave-one-out (LOO) cross-validation algorithm showed that the FHMM method provides better recognition of ALS and healthy subjects compared to standard HMM. Moreover, comparing our method with a state-of-the art method named least square support vector machine (LS-SVM) showed the efficiency of the FHMM in distinguishing ALS subjects from healthy ones.

Right ventricle functional parameters estimation in arrhythmogenic right ventricular dysplasia using a robust shape based deformable model.

Assessment of cardiac right-ventricle functions plays an essential role in diagnosis of arrhythmogenic right ventricular dysplasia (ARVD). Among clinical tests, cardiac magnetic resonance imaging (MRI) is now becoming the most valid imaging technique to diagnose ARVD. Fatty infiltration of the right ventricular free wall can be visible on cardiac MRI. Finding right-ventricle functional parameters from cardiac MRI images contains segmentation of right-ventricle in each slice of end diastole and end systole phases of cardiac cycle and calculation of end diastolic and end systolic volume and furthermore other functional parameters. The main problem of this task is the segmentation part. We used a robust method based on deformable model that uses shape information for segmentation of right-ventricle in short axis MRI images. After segmentation of right-ventricle from base to apex in end diastole and end systole phases of cardiac cycle, volume of right-ventricle in these phases calculated and then, ejection fraction calculated. We performed a quantitative evaluation of clinical cardiac parameters derived from the automatic segmentation by comparison against a manual delineation of the ventricles. The manually and automatically determined quantitative clinical parameters were statistically compared by means of linear regression. This fits a line to the data such that the root-mean-square error (RMSE) of the residuals is minimized. The results show low RMSE for Right Ventricle Ejection Fraction and Volume (≤ 0.06 for RV EF, and ≤ 10 mL for RV volume). Evaluation of segmentation results is also done by means of four statistical measures including sensitivity, specificity, similarity index and Jaccard index. The average value of similarity index is 86.87%. The Jaccard index mean value is 83.85% which shows a good accuracy of segmentation. The average of sensitivity is 93.9% and mean value of the specificity is 89.45%. These results show the reliability of proposed method in these cases that manual segmentation is inapplicable. Huge shape variety of right-ventricle led us to use a shape prior based method and this work can develop by four-dimensional processing for determining the first ventricular slices.

A novel approach to predict sudden cardiac death (SCD) using nonlinear and time-frequency analyses from HRV signals.

Investigations show that millions of people all around the world die as the result of sudden cardiac death (SCD). These deaths can be reduced by using medical equipment, such as defibrillators, after detection. We need to propose suitable ways to assist doctors to predict sudden cardiac death with a high level of accuracy. To do this, Linear, Time-Frequency (TF) and Nonlinear features have been extracted from HRV of ECG signal. Finally, healthy people and people at risk of SCD are classified by k-Nearest Neighbor (k-NN) and Multilayer Perceptron Neural Network (MLP). To evaluate, we have compared the classification rates for both separate and combined Nonlinear and TF features. The results show that HRV signals have special features in the vicinity of the occurrence of SCD that have the ability to distinguish between patients prone to SCD and normal people. We found that the combination of Time-Frequency and Nonlinear features have a better ability to achieve higher accuracy. The experimental results show that the combination of features can predict SCD by the accuracy of 99.73%, 96.52%, 90.37% and 83.96% for the first, second, third and forth one-minute intervals, respectively, before SCD occurrence.

How can selection of biologically inspired features improve the performance of a robust object recognition model?

Humans can effectively and swiftly recognize objects in complex natural scenes. This outstanding ability has motivated many computational object recognition models. Most of these models try to emulate the behavior of this remarkable system. The human visual system hierarchically recognizes objects in several processing stages. Along these stages a set of features with increasing complexity is extracted by different parts of visual system. Elementary features like bars and edges are processed in earlier levels of visual pathway and as far as one goes upper in this pathway more complex features will be spotted. It is an important interrogation in the field of visual processing to see which features of an object are selected and represented by the visual cortex. To address this issue, we extended a hierarchical model, which is motivated by biology, for different object recognition tasks. In this model, a set of object parts, named patches, extracted in the intermediate stages. These object parts are used for training procedure in the model and have an important role in object recognition. These patches are selected indiscriminately from different positions of an image and this can lead to the extraction of non-discriminating patches which eventually may reduce the performance. In the proposed model we used an evolutionary algorithm approach to select a set of informative patches. Our reported results indicate that these patches are more informative than usual random patches. We demonstrate the strength of the proposed model on a range of object recognition tasks. The proposed model outperforms the original model in diverse object recognition tasks. It can be seen from the experiments that selected features are generally particular parts of target images. Our results suggest that selected features which are parts of target objects provide an efficient set for robust object recognition.

Modeling the gait of normal and Parkinsonian persons for improving the diagnosis.

In this study, we present a model for the gait of normal and Parkinson's disease (PD) persons. Gait is semi-periodic and has fractal properties. Sine circle map (SCM) relation has a sinusoidal term and can show chaotic behaviour. Therefore, we used SCM as a basis for our model structure. Moreover, some similarities exist between the parameters of this relation and basal ganglia (BG) structure. This relation can explain the complex behaviours and the complex structure of BG. The presented model can simulate the BG behaviour globally. A model parameter, Ω, has a key role in the model response. We showed that when Ω is between 0.6 and 0.8, the model simulates the behaviour of normal persons; the amounts greater or less than this range correspond to PD persons. Our statistical tests show that there is a significant difference between the Ω of normal and PD patients. We conclude that Ω can be introduced as a parameter to distinguish normal and PD persons. Additionally, our results showed that Spearman correlation between the Ω and the severity of PD is 0.586. This parameter may be a good index of PD severity.

Identification of voice disorders using long-time features and support vector machine with different feature reduction methods.

Identification of voice disorders has a fundamental role in our life nowadays. Therefore, many of these diseases must be diagnosed at early stages of occurrence before they lead to a critical condition. Acoustic analysis can be used to identify voice disorders as a complementary technique with other traditional invasive methods, such as laryngoscopy. In this article, we followed an extensive study in the diagnosis of voice disorders using the statistical pattern recognition techniques. Finally, we proposed a combined scheme of feature reduction methods followed by pattern recognition methods to classify voice disorders. Six classifiers are used to evaluate feature vectors obtained by principal component analysis or linear discriminant analysis (LDA) as feature reduction methods. Furthermore, individual, forward, backward, and branch-and-bound methods are examined as feature selection methods. The performance of each combined scheme is evaluated in terms of the accuracy, sensitivity, specificity, and area under the receiver operating characteristic curve (AUC). The experimental results denote that LDA along with support vector machine (SVM) has the best performance, with a recognition rate of 94.26% and AUC of 97.94%. Additionally, this structure has the lowest complexity in comparison with other architectures. Among feature selection methods, individual feature selection followed by SVM classifier shows the best recognition rate of 91.55% and AUC of 95.80%.

Uterine segmentation and volume measurement in uterine fibroid patients' MRI using fuzzy C-mean algorithm and morphological operations.

BACKGROUND: Uterine fibroids are common benign tumors of the female pelvis. Uterine artery embolization (UAE) is an effective treatment of symptomatic uterine fibroids by shrinkage of the size of these tumors. Segmentation of the uterine region is essential for an accurate treatment strategy. OBJECTIVES: In this paper, we will introduce a new method for uterine segmentation in T1W and enhanced T1W magnetic resonance (MR) images in a group of fibroid patients candidated for UAE in order to make a reliable tool for uterine volumetry. PATIENTS AND METHODS: Uterine was initially segmented using Fuzzy C-Mean (FCM) method in T1W-enhanced images and some morphological operations were then applied to refine the initial segmentation. Finally redundant parts were removed by masking the segmented region in T1W-enhanced image over the registered T1W image and using histogram thresholding. This method was evaluated using a dataset with ten patients' images (sagittal, axial and coronal views). RESULTS: We compared manually segmented images with the output of our system and obtained a mean similarity of 80%, mean sensitivity of 75.32% and a mean specificity of 89.5%. The Pearson correlation coefficient between the areas measured by the manual method and the automated method was 0.99. CONCLUSIONS: The quantitative results illustrate good performance of this method. By uterine segmentation, fibroids in the uterine may be segmented and their properties may be analyzed.