1. Introduction
The Hong Kong population is aging rapidly such that the elderly population (those aged 65 years or older) is estimated to increase from approximately 13% in 2009 to 28% by 2039, there will be 2.5 million elder people in Hong Kong [1]. Since dementia is an age-related condition, the prevalence of dementia is expected to rise rapidly, and the consequent has an extensive impact on individuals, families, health and social welfare systems [2]. The number of dementia cases is projected to increase from 85,012 elder people in 2009 to 271,320 people in 2039. This is an important issue to be tackled and a dementia strategy for Hong Kong needs to be formulated. In order to cope with it, one direction is to develop accurate diagnostic tools, particular at early stage when treatment outcomes are most effective. Nuclear imaging and biomarkers diagnosis are the promising approach.
1.1. Definition of Dementia
Dementia is a group of diseases that primarily affect older people; but it is not a normal part of aging. The World Health Organization’s International Classification of Disease (1992-1994) defines dementia as:
“Dementia is a syndrome due to disease of the brain, usually of a chronic or progressive nature, in which there is disturbance of multiple higher cortical functions, including memory, thinking, orientation, comprehension, calculation, learning capacity, language and judgment. The impairments of cognitive function are commonly accompanied, and occasionally preceded by deterioration in emotional control, social behavior, or motivation”.
1.2. Types and Causes of Dementia [3]
There are numerous different forms of dementia and each has its own causes. Dementia can also be a symptom that develops in the later stages of other illnesses. The common types and causes of dementia include:
・ Alzheimer’s disease (AD) [4] accounts for 50% - 70% of all dementia. It is a progressive degeneration disease that shrinks brain cells, builds up tangles in soma and senile plaques, a neuritic process around an amyloid core, in extracellular space. These damage neurons connections disrupt messages within the brain.
・ Vascular dementia (VaD) [5] is the second (15% - 20%) most common form of dementia, which is caused by cerebrovascular conditions, usually caused by a series of small stroke. As the result, symptoms may not be noticeable or may only be temporary.
・ Mixed dementia (MD) [6] is a condition of which pathologies of AD and VaD coexist. Around 45% of AD patients exhibit significant vascular pathology.
・ Lewy bodies dementia (LBD) [7] is caused by abnormal microscopic deposits of spherical structuring Lewy bodies in the cerebral cortex.
・ Frontotemporal Lobar Degeneration (FTLD) [8] damages the cells in the frontal lobes, anterior temporal lobes or both, giving rise to a spectrum of dementing disorders. It has a high familial incidence with approximately 40% of patients having a similar disorder history in a first degree-relative.
・ Mild Cognitive Impairment (MCI) [9 , 10] is considered as a prodromal stage of dementia, with noticeable cognitive changes, but no significant disability to interfere with daily life.
1.3. Assessments of Dementia [11]
There is no single diagnostic test for all dementias and a comprehensive assessment (such as clinical history taking, cognitive and physical examination, medication review, brain scanning, etc.) is required in order to exclude treatable causes of cognitive impairment.
Common assessments include:-
・ Mini-mental state examination (MMSE), using in initial checkup, to analyze a subject’s cognitive status. It uses a set of questions to provide measures of orientation, registration, short-term memory, motor and language functioning.
・ Routine heametology and biochemistry tests to detect untreated syphilis, abnormality of thyroid function and decline of B12.
・ MRI/CT brain scan to show the level of any atrophy and any vascular disease or infarcts.
However, absolute confirmation of the diagnosis can only be made by a biopsy and an autopsy of brain tissues.
2. 18F-FDG PET-CT Brian Scan and Cortex-ID
Numerous studies have found that AD and other dementias could alter brain metabolisim significantly. AD patients usually present the posterior cingulate and parietotemporal cortices hypometabolisim and spread into the frontal lobes in advanced stage. In contrast, FTLD patients show manifestly hypometabolism in the frontal and anterior temporal cortices, while DLB patients present hypometabolism in the posterior brain comprising primarily the parietoocipital regions. Theoretically, 18F-FDG PET scan can help in early diagnosis of AD and other dementias by highlighting these decreased FDG uptake cortex regions before MRI or CT scans can detect any structural brain damage.
To quantify the cortex metabolic activity, the reconstructed FDG PET images of the patient are processed with CortexID, developed by Dr. Satoshi Minoshima [12 , 13]. The software provides automated brain assessment through quantification of the comparison of local peak activity values, at standardized anatomical locations, with the corresponding reference normal peak activity in age-stratified control subjects. 3D stereotactic surface projections of the brain are generated and 31 regional Z-scores are calculated. A positive Z-score is given in a hypometabolic CortexID area.
By using these 31 regional Z-scores with MRI cerebrovascular condition, two ANNs were designed for dementias differentiation and MMSE score mapping. The details will be presented in section 4.
3. Multilayer perceptron (MLP) [14 - 16]
ANN is a mathematical model to simulate the structure and functionalities of biological neural networks. The basic building block, artificial neuron, consists of multiplying every input value with individual weight, summing all weighted values and activating through the activation function (Figure 1 and Figure 2). Its mathematical description is:
(1)
where
,
and
are input, weight and output values in discrete time k respectively,
b is bias,
g is a activation function.
A multilayer perceptron (Figure 3) is a feedforward artificial neural network and consists at least of three layers of nodes (artifical neuron): an input layer, a hidden layer and an output layer. It can divide the multi-dimensional feature space into arbitrartly complex decision regions and represent any boolean function.
Figure 1. Working principle of an artifical neuron.
Figure 2. Different types of activation functions: (a) Threshold, (b) Piecewise linear, (c) Sigmoid, and (d) Gaussian.
Therefore MLP is widely and successfully used for pattern classification.
MLP utilizes the backpropagation rule for training. The rule repetitively calculates an error function for each input and backpropagates the error from one layer to the previous one. The weights are adjusted in direct proportional to the error in the nodes to which it is connected. The back-propagation algorithm is described below:
Step:
1) Initialize the weights randomly.
2) Arbitrarily select an input X.
3) Feed the signal through the network.
4) Compute
in the output layer
(2)
where
is the net input to the ith neuron in the Lth layer,
is the derivative of the activation function,
and
are the desired and feed forward output respectively.
5) Back-propagate the errors and compute the preceding layers deltas
(3)
6) Update weights
where
(4)
7) Repeat step 2 until the error in the output layer is below a prespecified values.
4. Building NNs of using CortexID to Model Physician’s Diagnosis and MMSE Score Mapping
MLPs are universal function approximators so they can be used to create mathematical models for function mapping as well as classifier. Therefore, MLPs were chosen in our application. Since too much hidden layers results in memorizing and not learning, one hidden layer is sufficient for the tasks. Except linear output neurons were used in MMSE score mapping, all neurons employ sigmoid activition function and trained with back-propagation algorithm.
4.1. Data Preparation
108 PET-CT examination radiological reports from University of Hong Kong at Queen Mary Hospital were screened. 87 reports are dementia cases and 17 of them were excluded because of incomplete CortexID data or unclassified or only having a single incidence. Moreover, 55 reports have MMSE scores and the frequency distribution of dementia is tabulated in Table 1.
It is observed that the distribution is severe skewed which will provide poor accuracies in most ANNs. To cope with the problems, class distributions need to be normalized by creating synthetic samples with Z-scores perturbed randomly [17]. Since VaD can only be accurately detected by MRI/CT imaging, the cerebrovascular condition becomes one of the input attributes of neural network in addition with 31 CortextID Z-scores.
Table 1. The dementia distribution in the dataset for classifer training.
4.2. Training Procedure
We used Matlab 2010 neural toolbox to build the ANNs. The data was divided randomly into the training set, the validation set and the test set. The training set is used for computing the gradient and updating the network weights. The validation set is used to improve generalization and to stop training when the network starts to overfit the data. The test set error is for different models comparison. For the dementia classifier training, the data is partitioned with a 60% training set, 20% validation set and 20% test set. In MMSE score mapping, the division portions were adjusted to 70%, 10% and 20% respectively.
4.3. Trained NN Models
In order to determine the optimum number of hidden neurons, the amount of neurons in the hidden layer is progressive increasing until the errors of the training, validation and test data sets are below the prespecified threshold.
The obtained NN models for dementia classifier (Figure 4) and MMSE score mapper (Figure 5) are shown.
5. Results
5.1. Dementia Classifier
The performance of the classifier can be visualized in the confusion matrices (Figure 6(a)). 7.4% of AD and VaD cases were misclassified as MCI and MD respectively. Moreover, since the ROC curves (Figure 6(b))
(a)(b)(c)
Figure 6. (a) The confusion matrices [1-MCI, 2-AD, 3-VaD, 4-MD, 5-LBD]; (b) the receiver-operating characteristics and (c) mean squared error of the dementia classifier.
are very close to the upper left corner, this implies the overall accuracy of the classifier is very high. Figure 6(c) reflects that the network was continued adjusting the weights to classify the data pattern and the best validation performance was obtained at epoch 45 with mean square error 0.021872.
5.2. MMSE Score Mapping
The correlations between the outputs of mapping function network and desired MMSE target values were plotted in Figure 7. The root means squared errors between output and target of training, validation and testing datasets are 1.11, 1.74 and 6.62 respectively. An outlier occurred in the test-set upshifts the regression line significantly. Moreover, the overall regression R-value is 0.87.
6. Performance Analysis
To prevent biased, untrained data sets are used to calculate the sensitivities and specificities of the dementia classifier. Other related performance indices are also computed and tabulated in Table 2.
Based on the analysis table, it is observed that the dementia classifier was learnt to categorise different types of dementia very accurately as a radiologist by using 31 regional CortexID Z-scores with cerebrovascular condition.
In MMSE score mapping study, the sample sizes of validation and testing data sets are 6 and 11; their regression values are 0.983 and 0.646. Therefore, the combined regression value is
Figure 7. Regression Plots between MMSE Mapping Function outputs and targets.
Table 2. Performance Analysis of Dementia Classifier on the Untrained Data Set (Validation + Test).
Since the
is very close to 1, the mapping network correlates highly with MMSE score.
7. Conclusions
Two multilayer perceptrons, based on 31 regional CortexID Z-scores, were designed and trained to classify different types of dementia and map with MMSE score. Comparing with physicians diagnosis, the dementia classifier achieved 96% accuracy. The mapper network can correctly predict the MMSE score with 0.782 regression value.
However, due to the limit number of records, we do not know whether the network can be trained to differentiate other types of dementias, such as frontotemporal dementia, Parkinsonian syndromes, etc. Further study and verification are required.
AcknowledgementS
This is a retrospective chart review study and is supported by the Department of Diagnostic Radiology, University of Hong Kong and the Hospital Authority, Hong Kong. We also wish to acknowledge the clinical data support of Prof. Leung-Wing Chu of the Department of Medicine of the University of Hong Kong.