eHealth Solutions for the Integrated HealthcareView this Special Issue
Research Article | Open Access
Andrius Lauraitis, Rytis Maskeliūnas, Robertas Damaševičius, "ANN and Fuzzy Logic Based Model to Evaluate Huntington Disease Symptoms", Journal of Healthcare Engineering, vol. 2018, Article ID 4581272, 10 pages, 2018. https://doi.org/10.1155/2018/4581272
ANN and Fuzzy Logic Based Model to Evaluate Huntington Disease Symptoms
We introduce an approach to predict deterioration of reaction state for people having neurological movement disorders such as hand tremors and nonvoluntary movements. These involuntary motor features are closely related to the symptoms occurring in patients suffering from Huntington’s disease (HD). We propose a hybrid (neurofuzzy) model that combines an artificial neural network (ANN) to predict the functional capacity level (FCL) of a person and a fuzzy logic system (FLS) to determine a stage of reaction. We analyzed our own dataset of 3032 records collected from 20 test subjects (both healthy and HD patients) using smart phones or tablets by asking a patient to locate circular objects on the device’s screen. We describe the preparation and labelling of data for the neural network, selection of training algorithms, modelling of the fuzzy logic controller, and construction and implementation of the hybrid model. The feed-forward backpropagation (FFBP) neural network achieved the regression R value of 0.98 and mean squared error (MSE) values of 0.08, while the FLS provides a final evaluation of subject’s reaction condition in terms of FCL.
Huntington disease (HD) is a progressive genetic neurodegenerative disorder causing involuntary movement and cognitive problems that significantly affect daily life of HD patients. HD affects about 1 in 10,000 to 20,000 people of European (Caucasian) descent , though in some isolated populations it is much higher. HD reduces life expectancy due to heart disease, pneumonia, physical injury from falls, and suicide. The most visible symptom of HD is chorea, which consists of jerky, involuntary movements of the upper and lower extremities, face or body, and occurs in about 90% of patients at some stage of their illness . Other symptoms include behavioural problems, cognitive impairment, psychiatric disorders, and dementia, which have a serious impact on daily living of a patient and often result in hospitalization. The societal and financial cost of HD on health and social care systems is significant and is estimated to be £195 million per year in the UK alone .
HD is currently incurable so most of the current research in this area focuses on identifying the deficits at the early stage of the disease, to benefit from future medical interventions that may help delaying the progress of the disease [4, 5]. This is also the case of the work presented in this paper. Traditional HD research often include magnetic resonance neuroimaging (MRI) measures of striatum and white matter volume, CAG repeat length in chromosome analysis, age, and striatal atrophy [6, 7]. Moreover, medical personnel and doctors who have experience in caring after HD patients and knowing that disease is cureless are not usually motivated to conduct scientific research themselves or to support multidiscipline (e.g., bioinformatics) investigations.
Any scientific result (device, technology, and theoretical model) that could contribute towards improvement of daily life of HD patient’s and help to monitor or predict the progress of the disease can be useful for both doctors and HD patients.
The problematics of data prediction evolved with the rise of artificial intelligence (AI) and machine learning (ML) methods and algorithms. Artificial neural networks (ANN) such as multilayer perceptron (MLP) can be used for classification of accelerometer-based tremor signals invoked by Parkinson patient’s involuntary movements . Prediction of Parkinson disease onset by adapting radial basis function neural network (RBFNN) for tremor activity data recorded via stimulation electrodes using electromyography (EMG) signals is described in . Dynamic neural network (DNN) is used to detect time-varying occurrences of tremor and dyskinesia from time series data acquired from EMG sensors and triaxial accelerometers worn by Parkinson patients . Another approach of designing a prediction model for Parkinson’s disease uses a decision tree and Iterative Dichotomiser (ID3) methods to analyze data collected from HD symptoms such as trembling in the legs, arms, hands, impaired speech articulation, and production difficulties . Hybrid models combine different AI and ML approaches for reproducing intelligent human reasoning process . By using information fusion, hybrid models combine heterogeneous ML approaches and improve quality of reasoning for complex regression and classification problems . Neurofuzzy systems combine neural network and fuzzy logic paradigms to avoid the limitations of neural network explanations to reach decision and limitations of fuzzy logic to automatically acquire the rules used for making those decisions . Fuzzy expert systems such as neurofuzzy system (ANFIS) can be applied in assessment of Parkinson’s disease with a noninvasive screening system for quantitative evaluation and analysis by using amplitude, frequency, spectral characteristics, and trembling localization parameters of input data . Hybrid model is adapted in designing a decision support system (DSS) for the intelligent identification of Alzheimer where neurofuzzy system explores approximation techniques from neural networks to find the parameter of a fuzzy system . Hybrid systems are also used as a classifier fusion strategy (Bayesian, SVM, k-nearest neighbours) in the prevalence of age-related diseases like Alzheimer’s and dementia , diagnostics and measurement  with wavelet transform (WT) and norm entropy feature extraction methods. The DSS that uses MLP and RBFNN is applied for monitoring patients with neurological disorders . The data is collected using noninvasive smart devices (modified mouse and 3-axis accelerometer sensor). Integration of neurofuzzy networks and information fusion for multimodal human cognitive state recognition is described in . Projection-based learning for metacognitive radial basis function network (PBL-McRBFN) is applied to predict Parkinson’s disease . Other hybrid systems and applications include nonlinear adaptive system, which fuses brain and gait information algorithmically using multistate Markov model . Accurate Parkinson disease diagnosis model based on cluster analysis uses random tree, classification and regression tree (C-RT), ID3, binary logistic regression, k-NN, partial least square regression (PLS), support vector machines (SVM) , and fuzzy c-means clustering (FCM) . Table 1 provides a summary of methods used by other authors.
Our previous work included the development of text input-based system for evaluating the condition of Huntington’s patients . The use of ANN for predicting the functional capacity of a Huntington’s patient was proposed in .
The aim of this paper is to create a computerized behavioural model, which predicts an impaired reaction condition for HD patients. We develop a mobile application to collect a dataset using finger touch coordinates and reaction time features extracted from test subjects (healthy and HD patients); create an ANN to predict the functional capacity level and fuzzy logic system (FLS) to determine the reaction condition (stage) for individual person; combine ANN with FLS into a hybrid model to predict the impaired reaction condition for HD patients; and simulate an experimental setup for test subjects to perform a provided exercise (test) at the different moments in time in order to predict a possibly impairing reaction condition with the help of the proposed hybrid model.
2. Materials and Methods
The study included ten (10) Huntington disease (HD) patients living in Lithuania. Each HD patient agreed to participate and allowed the data collected during the test to be used for scientific purposes. Every HD patient fall in the early clinical descriptor category of Huntington disease, that is, I and II stages according to Shoulson–Fahn evaluation system . Such HD patients have hand tremors, body movement distractions, but are capable to perform the test on a mobile application without extra help, for example, from medical personnel, nurses, or family members. Other ten (10) participants were healthy people with no signs of any neurological or neurodegenerative disorder.
The test can be performed using various mobile devices that support Android OS. The mobile application randomly generates circular shape objects (2, 3, and 5 circles at time) of particular color that are generated on the mobile device’s screen. Each circle is located in different positions of the screen, thus no possible collisions (overlapping) between two particular circles are possible. An active circle that needs to be touched is marked by a black contour so as to differ from other objects.
The subjects are instructed to touch every object, starting from first in sequence, by finger as close to center and as quickly as possible. When subject finishes the test, collected data is stored in external mobile device storage and sent to the database using the internet connection.
The collected dataset consists of 3032 data examples collected from 20 test subjects (10—healthy and 10—HD patients). The dataset (see a sample in Table 2) contains the ground truth coordinates of the generated object, the coordinates of subject’s touch, subject’s reaction time, subject’s label, and the marker of Huntington’s disease.
x,y: screen coordinates of the center of circular object to touch; xt, yt: screen coordinates of user touch; nC: number of circular objects rendered on the device screen; rt: user’s reaction time in seconds; delta: the Euclidean distance between object’s center and touch position. User: user ID; IsSick: indicates if test subject has Huntington disease (1 yes, 0 otherwise).
2.4. Feature Extraction and Class Labelling
The subject’s reaction time (rt) and the Euclidian distance between the two points of true and touched positions (delta) serve as features which are incorporated as input variables to ANN. We assume that smaller rt and delta values indicate better functional capacity level. The bigger delta value can show stronger hand tremoring, whereas higher rt value is an indicator of body stagnancy.
The statistical analysis of the rt and delta values has revealed that the values are not normally distributed, but after the applying the log transformation, which is commonly used in regression analysis of biological data with highly skewed distribution , the values become normal as confirmed by visual inspection in Figure 1 and skewness and kurtosis tests ( and ). For data samples greater than 300, values and are considered as acceptable for normality .
To analyze the power of rt and delta values to correctly predict the healthy or sick state of the subject, we have performed feature evaluation using the relative entropy (also known as the Kullback–Leibler distance or divergence) criterion, considering different number of objects presented at the screen. The results are presented in Figure 2. In all cases, delta feature has larger discriminative power than rt, and the features from 3 and 5 objects test are more statistically discriminative.
2.5. ANN for Functional Capacity Level Prediction
We have analyzed the following neural network models: (1) feed-forward backpropagation (FFBP); (2) feed-forward time delay neural network (FFTD); (3) cascade-forward backpropagation (CFBP); (4) nonlinear autoregressive exogenous model (NARX); (5) Elman neural network; (6) layer recurrent neural network (RNN); and (7) generalized regression neural network (GRNN).
FFBP is a simple neural network without any cycle connections between neurons . FFTD has no internal state and adds delayed copies as other inputs as an input signal to obtain time-shift invariance . In CFBP, the input values calculated after every hidden layer are backpropagated and the weights adjusted . NARX have a limited feedback, which comes only from the output neuron rather than from hidden layer . Elman network additionally has context units, which are connected to the hidden units, thus providing the network with memory . RNN represent an architecture where connections between units form a directed cycle . GRNN has only one (smoothness) parameter, and its convergence is guaranteed; fast and stable .
Each neural network has 2 inputs (rt, delta) and 1 output (Y). Neural network is composed of single neurons that are treated as a simple unit carrying signals (data) to each other or different layers via transfer functions, which correspond to sum of input signal. Training function is the optimization algorithm used for finding global minimum of a function. The outputs of ANN are class labels for determining the functional capacity of a person (the larger value indicates that a person is more capable to do motoric activities). Such scenario imitates the TFC scale measurement system for Huntington disease patients presented in Table 3 .
Table 4 illustrates the setup for analyzed ANN models with their parameters.
2.6. Training and Testing
The dataset was randomly divided into 3 sets: training, validation, and testing. Training set uses all samples from 70% of users. Validation set (15%) is used to measure network generalization and to stop training when necessary. Testing set (15%) provides independent performance of the network afterwards. We also analyzed a different partition of the dataset (40% for training, 30% for validation, and 30% for testing); however, there were no significant differences in the performance of ANN.
Overfitting was prevented by using the early stopping technique, which controls error on the validation set which is monitored during training process: when error increases for a specified number of iterations then the training is stopped and the weights and biases at the minimum of the validation error are returned.
For each neural network model, we have repeated the training and testing process for 20 times in order to allow calculation of statistical characteristics (mean, standard deviation) of ANN performance measures and to perform statistical comparison.
2.7. Reaction Stage Determination Using Fuzzy Logic System (FLS)
The aim of the FLS system is to determine the reaction stage of a patient (test subject) according to some predefined parameters. The FLS consists of three main parts: fuzzification block, inference mechanism, and defuzzification block. Membership functions, linguistic variables are created in fuzzification module. Inference engine is responsible for applying logical rules (fuzzy rule base) to the knowledge base and deduce new knowledge. Defuzzification module converts all the fuzzy terms created by the rule base of the controller to crisp terms (numerical values). The FLS uses triangular membership Mamdani-type functions with fuzzy set inference mechanism (minimum implication, maximum aggregation, minimum AND operator, maximum OR operator) and centroid defuzzification method.
The parameters of the FLS are derived from the ANN output corresponding to the functional capacity level, so in the FLS design process, the model input and output values need to be considered accordingly. There are three input and one output variable in the FLS. The input parameters are AVG1, AVG2, and AVG3, which correspond to the average of ANN output values when test subject is working with two, three, and five objects, respectively. All three inputs can have values in range [0; 10].
The linguistic variables (terms) for AVG1, AVG2, and AVG3 are (1)LOW [0 2 4];(2)AVERAGE [3.6 5.5 7];(3)HIGH [6.6; 8.5 10].
The model has one output parameter ReactionStage can have one of five values: close to peaks 1, 3, 5, 7, or 9, that is, each peak corresponds to particular linguistic variable of ReactionStage. The terms for output parameter ReactionStage are (1)ADVANCED [0 1 2];(2)LATE [1.5 3 4];(3)AVERAGE [3.5 5 6];(4)EARLY [5.5 7 8];(5)HEALTHY/PRECLINICAL [7.5 9 10].
The FLS rule base is formed from 27 fuzzy rules. Table 5 illustrates the principles of constructing fuzzy rule base. These can be interpreted as general fuzzy IF-THEN rules containing only fuzzy logical AND operators, for example,
IF AVG1 is LOW AND AVG2 is LOW AND AVG3 is LOW
THEN ReactionStage is ADVANCED.
3. Proposed Hybrid Model
The hybrid model (see Figure 3) is composed of four sub models: (1) dataset formation; (2) ANN prediction model; (3) fuzzy logic expert system (FLS); and (4) decision module for determination of person’s condition.
During dataset formation, test subjects (under the supervision of a healthy person—a medical doctor or a nurse) use smart devices to perform reaction and accuracy test experiments with their fingers. The collected data is stored in the database. The ANN submodel predicts the functional capacity level of a person using the data from the database. The network is trained by observing regression (R), that is, correlation measurement between outputs and targets and mean squared error (MSE) values. Once the network is trained, it can make predictions on new sample data. Finally, to evaluate the reaction condition of a test subject, the test session is repeated at a different time and the ANN predictions are aggregated, and the reaction stage of a person is evaluated using a fuzzy rules system.
4. Experimental Results
The hybrid model was implemented with MATLAB Neural Network and Fuzzy Logic Toolbox software (MathWorks Inc.). The results of regression and comparison of the prediction results of the analyzed ANN models is presented in Table 6, whereas the performance of neural networks in terms of means and 95% confidence intervals of R and MSE is given in Table 7. The “TFC” field indicates the ground truth evaluation of the patient state provided by a medical neurologist expert according to the TFC scale. The R metric measures the correlation between output and targets, whereas the MSE metric is the average squared difference between outputs and targets.
Nonparametric Friedman test was conducted to compare the performance results (MSE) among ANN models. Results show that there is a significant difference in performance among all ANN models (chi-square = 133.15; ). Posthoc Nemenyi tests further reveal that the performance of FFBP is the best among all ANN models (Figure 4).
Figure 5 shows an example of FFBP best performance equal to R = 0.993 and MSE = 0.094 on the validation set.
Table 8 illustrates impaired reaction condition simulation example on a single test subject using the FLS system. In order to make comparison, data samples were collected at different time moments. Feature (rt1, delta1, rt2, and delta2) values are presented in all three modes (10 attempts), thus giving two separate ANN (in the example provided, FFBP model was used) prediction outputs, which are used to calculate average values and evaluate the reaction condition in the FLS.
We have presented an actual experimental framework to assess finger-tapping tests performed by patients suffering from the Huntington’s disease (HD). The proposed model was validated using a dataset of 3032 data records collected from 20 test subjects (both healthy and HD patients). The reaction condition was determined using the developed Mamdani Type-1 fuzzy logic expert system (FLS) with 3 input (3 linguistic variables), 1 output (5 linguistic variables), triangular membership functions, and 27 fuzzy rules base.
We describe an architecture that combines several artificial neural networks (ANN) of different type (FFBP, FFDT, CFBP, NARX, Elman, RNN, and GRNN) to create a hybrid (neurofuzzy) model, which integrates feature extraction, prediction, and classification routines to forecast the impaired reaction condition for HD patients. The best results were achieved using the feed-forward backpropagation (FFBP) neural network model, which predicts the total functionality capability (TFC score) with high performance results, that is, it has obtained regression R value not less than 0.98 and mean squared error (MSE) values of 0.08, while FLS evaluates several measurements taken time apart to provide a final evaluation of the subject’s reaction condition.
Future work will focus on the validation of the proposed system using a larger dataset, which includes the data collected from the Parkinson’s and Alzheimer’s patients as well, the analysis and use of more sophisticated finger-tapping features, and the comparison of the ANN results with those of SVM regression.
Human Studies. Research on human subjects was approved by the Institutional Review Board of the Faculty of Informatics of Kaunas University of Technology.
Conflicts of Interest
The authors declare that there is no conflict of interest regarding the publication of this article.
Authors would like to thank the president of Lithuania Huntington disease association, Dr. Zivile Navikiene, for contacting family members of HD patients to help carry out experiments for the investigation described in this paper, as well as for practical support and advice.
- R. A. C. Roos, “Huntington’s disease: a clinical review,” Orphanet Journal of Rare Diseases, vol. 5, no. 1, p. 40, 2010.
- M. S. Haddad and J. L. Cummings, “Huntington’s disease,” The Psychiatric Clinics of North America, vol. 20, no. 4, pp. 791–807, 1997.
- C. Jones, M. Busse, L. Quinn et al., “The societal cost of Huntington’s disease: are we underestimating the burden?” European Journal of Neurology, vol. 23, no. 10, pp. 1588–1590, 2016.
- M. Bennasar, Y. Hicks, S. Clinch et al., “Huntington’s disease assessment using tri axis accelerometers,” Procedia Computer Science, vol. 96, pp. 1193–1201, 2016.
- T. V. Wiecki, C. A. Antoniades, A. Stevenson et al., “A computational cognitive biomarker for early-stage Huntington’s disease,” PLoS One, vol. 11, no. 2, article e0148409, 2016.
- J. S. Paulsen, D. R. Langbehn, J. C. Stout et al., “Detection of Huntington’s disease decades before diagnosis: the Predict-HD study,” Journal of Neurology, Neurosurgery & Psychiatry, vol. 79, no. 8, pp. 874–880, 2008.
- E. H. Aylward, B. F. Sparks, K. M. Field et al., “Onset and rate of striatal atrophy in preclinical Huntington disease,” Neurology, vol. 63, no. 1, pp. 66–72, 2004.
- M. Engin, S. Demirag, E. Engin et al., “The classification of human tremor signals using artificial neural network,” Expert Systems with Applications, vol. 33, no. 3, pp. 754–761, 2007.
- D. Wu, K. Warwick, Z. Ma et al., “Prediction of Parkinson’s disease tremor onset using a radial basis function neural network based on particle swarm optimization,” International Journal of Neural Systems, vol. 20, no. 02, pp. 109–116, 2010.
- B. T. Cole, S. H. Roy, C. J. De Luca, and S. H. Nawab, “Dynamic neural network detection of tremor and dyskinesia from wearable sensor data,” in 2010 Annual International Conference of the IEEE Engineering in Medicine and Biology, vol. 2010, pp. 6062–6065, Buenos Aires, Argentina, 2010.
- A. Chandrashekhar, S. Jain, and V. Kumar Jha, “Design and analysis of data mining based prediction model for Parkinson’s disease,” International Journal of Computer Science Engineering (IJCSE), vol. 3, no. 03, 2014.
- T. Furuhashi, S. Tano, and H.-A. Jacobsen, Deep Fusion of Computational and Symbolic Processing, vol. 59, Publish in Physica, Fuzziness and Soft Computing, New York, NY, USA, 2012.
- P. Kazienko, E. Lughofer, and B. Trawinski, “Hybrid and ensemble methods in machine learning,” Published in Journal of Universal Computer Science, vol. 19, no. 4, pp. 457–461, 2013.
- R. Fuller, “Introduction to neuro-fuzzy systems,” in Advances in Intelligent and Soft Computing, vol. 2, Springer Science & Business Media (2013), 2013.
- O. Geman, “Parkinson’s disease assessment using fuzzy expert system and nonlinear dynamics,” Advances in Electrical and Computer Engineering, vol. 13, no. 1, pp. 41–46, 2013.
- J. C. Obi and A. A. Imainvan, “Decision support system for the intelligent identification of Alzheimer using neuro fuzzy logic,” International Journal on Soft Computing (IJSC), vol. 2, no. 2, pp. 25–38, 2011.
- S. Iram, P. Fergus, D. Al–Jumeily, A. Hussain, and M. Randles, “A classifier fusion strategy to improve the early detection of neurodegenerative diseases,” International Journal of Artificial Intelligence and Soft Computing (IJAISC), vol. 5, no. 1, pp. 23–44, 2015.
- Z. J. Viharos and K. B. Kis, “Survey on neuro-fuzzy systems and their applications in technical diagnostics and measurement,” Measurement, vol. 67, pp. 126–136, 2015.
- I. Chiuchisan and O. Geman, “An approach of a decision support and home monitoring system for patients with neurological disorders using internet of things concepts,” WSEAS Transactions on Systems, vol. 13, pp. 460–469, 2014.
- G. Yang, Y. Lin, and P. Bhattacharya, “Multimodality inferring of human cognitive states based on integration of neuro-fuzzy network and information fusion techniques,” EURASIP Journal on Advances in Signal Processing, vol. 2008, no. 1, article 371621, 2007.
- G. Sateesh Babu and S. Suresh, “Parkinson’s disease prediction using gene expression – a projection based learning meta-cognitive neural classifier approach,” Expert Systems with Applications, vol. 40, no. 5, pp. 1519–1529, 2013.
- O. Geman, “Towards an inclusive Parkinson’s screening system,” in 2014 18th International Conference on System Theory, Control and Computing (ICSTCC), pp. 470–475, Sinaia, Romania, 2014.
- G. R. Ramani and G. Sivagami, “Parkinson disease classification using data mining algorithms,” International Journal of Computer Applications, vol. 32, no. 9, pp. 17–22, 2011.
- I. Rustempasic and M. Can, “Diagnosis of Parkinson’s disease using fuzzy c-means clustering and pattern recognition,” Southeast Europe Journal of Soft Computing, vol. 2, pp. 42–49, 2013.
- J. Gelsvartas, R. Simutis, and R. Maskeliunas, “User adaptive text predictor for mentally disabled Huntington’s patients,” Computational Intelligence and Neuroscience, vol. 2016, Article ID 3054258, 6 pages, 2016.
- A. Lauraitis and R. Maskeliūnas, “Investigation of predicting functional capacity level for Huntington disease patients,” in Information and Software Technologies. ICIST 2017. Communications in Computer and Information Science, R. Damaševičius and V. Mikašytė, Eds., vol. 756, pp. 142–149, Springer, Cham, 2017.
- I. Shoulson and S. Fahn, “Huntington disease: clinical care and evaluation,” Neurology, vol. 29, no. 1, pp. 1–3, 1979.
- S. Gustavsson, B. Fagerberg, G. Sallsten, and E. M. Andersson, “Regression models for log-normal data: comparing different methods for quantifying the association between abdominal adiposity and biomarkers of inflammation and insulin resistance,” International Journal of Environmental Research and Public Health, vol. 11, no. 12, pp. 3521–3539, 2014.
- H.-Y. Kim, “Statistical notes for clinical researchers: assessing normal distribution (2) using skewness and kurtosis,” Restorative Dentistry & Endodontics, vol. 38, no. 1, pp. 52–54, 2013.
- D. E. Rumelhart, G. E. Hinton, and R. J. Williams, “Learning representations by back-propagating errors,” Nature, vol. 323, no. 6088, pp. 533–536, 1986.
- A. Waibel, T. Hanazawa, G. Hinton, K. Shikano, and K. J. Lang, “Phoneme recognition using time-delay neural networks,” IEEE Transactions on Acoustics Speech and Signal Processing, vol. 37, no. 3, pp. 328–339, 1989.
- S. Goyal and G. K. Goyal, “Cascade and feedforward backpropagation artificial neural network models for prediction of sensory quality of instant coffee flavored sterilized drink,” Canadian Journal on Artificial Intelligence, Machine Learning and Pattern Recognition, vol. 2, no. 6, pp. 78–82, 2011.
- S. Billings, Nonlinear System Identification: NARMAX Methods in the Time, Frequency, and Spatio-Temporal Domains, Wiley, Chichester, UK, 2013.
- J. L. Elman, “Distributed representations, simple recurrent networks, and grammatical structure,” Machine Learning, vol. 7, no. 2-3, pp. 195–225, 1991.
- C. Goller and A. Küchler, “Learning task-dependent distributed representations by backpropagation through structure,” in Neural Networks, 1996, IEEE International Conference on, vol. 1, pp. 347–352, Washington, DC, USA, 1996.
- D. F. Specht, “A general regression neural network,” IEEE Transactions on Neural Networks, vol. 2, no. 6, pp. 568–576, 1991.
Copyright © 2018 Andrius Lauraitis et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.