Abstract
Keywords
Introduction
Mild Cognitive Impairment (MCI) is a prodromal stage of dementia, characterized by cognitive decline beyond that expected for age but without significantly impacting daily activities. 1 It is considered as an intermediate state between normal cognitive aging and dementia. Early diagnosis of MCI is crucial for planning potential interventions and management strategies. 2 The prevalence of MCI increases with age in both genders. Patel et al (2018) has reported that there is a prevalence of 15% in elderly people of age greater than 80 years, to 2.6% in people of age less than 60 years in India. 3
There is sufficient evidence for the association between the Heart Rate Variability (HRV) and MCI.4,5 The HRV refers to the variation in time elapsed between successive heartbeats. It is a feasible, non-invasive measure of influence of the Autonomic Nervous system (ANS) on heart rate regulation and it reflects the dynamic balance between sympathetic and parasympathetic nervous system activity. 6 Derby et al (2023), has proved that HRV is associated with global cognition and cognitive impairment. 7 Arakaki et al (2023) showed the evidence that changes in Heart Rate (HR) and HRV may serve as indicators of cognitive decline and the progression of neurodegenerative disorders. 8 Longitudinal studies have found associations between baseline HRV measures and subsequent cognitive decline or dementia diagnosis. 9 Altered or reduced HRV has been associated with different factors like autonomic dysfunction, 10 cardiovascular risk factors, inflammation, and oxidative stress, which ultimately affects the brain-heart axis.11,12
Liu et al (2023) have examined that AD patients with agitation (presence) have a lower HRV and/or a steeper fall in HRV over time than those without agitation. 13 Another study has illustrated that there is an increased risk of dementia when the RHR is greater than or equal to 80 beats per minute, regardless of vascular risk factors and underlying cardio-vascular diseases (CVDs). 14 Similarly, the participants with RHR greater than or equal to 70 bpm experienced a faster deterioration in general cognitive function than those with RHR of 60-69 bpm.
Deng et al (2022) conducted a comprehensive logistic regression analysis, revealing a significant association between RHR and cognitive function. Elevated RHR values, predominantly exceeding 80 beats per minute, were frequently observed among individuals experiencing cognitive decline. 15 The study included 4177 participants, encompassing 2354 cases of Alzheimer’s Disease (AD) and 989 instances of vascular dementia (VaD). Similarly, a meta-analytic study reported an association between increased RHR and all combined types of dementia. They suggest that, in order to ascertain whether elevated RHR is a direct cause of cognitive deterioration, large-scale prospective cohort studies are necessary. 16 Molloy et al (2023) found an association between higher RHR and lower Mini-Mental State Examination (MMSE) scores among individuals with pathological amyloid/tau ratios in the cerebrospinal fluid. 17
There exist a cross-sectional association between HRV and cognitive performance which was reported in a study conducted among community-dwelling elderly Mexican Americans, an ethnic group at high risk for cardiovascular disease risk factors. 18 Reduced HRV was reported to be associated with worse performance on the global test of cognitive function. MCI patients are in general found to have significantly lower HRV compared to the healthy controls (HC) at resting state. 19
Although there have been certain studies that unwrap the correlation between HR, HRV and cognitive performance at global scale, 20 the applicability of the association on Indian elderly population with MCI is unknown. The objective of this study is to investigate the correlation between HR, ultra-short HRV (usHRV), and cognitive decline, with a specific emphasis on those diagnosed with MCI within the elderly Indian population. Our aim is to provide a comprehensive analysis of the potential implications of HR and usHRV as biomarkers for early identification of cognitive impairment. We develop a signal processing pipeline that can detect R-wave peaks and calculate HR and usHRV from brief ECG recordings, for point-of-care/bedside applications.
Methods
Participants
The present study was conducted as part of the Centre for Brain Research-Tata Longitudinal Study for Ageing (CBR-TLSA) and approved by the Institutional Ethics Committee. The study protocol follows all necessary ethical guidelines including participants’ safety, confidentiality, informed consent and voluntary participation. The cohort underwent comprehensive clinical, biochemical, and radiological assessments at regular intervals. The included participants in the study were 45 years or older belonging to urban, working middle-class backgrounds, with reasonable literacy levels and having undergone significant lifestyle changes over the past few decades. The recruitment process for the cohort involved stringent exclusion criteria to ensure the robustness of the study. Participants with terminal medical illnesses, psychiatric disorders, severe visual or hearing impairments, or locomotor disabilities likely to interfere with study assessments were excluded from the cohort.
The study participants were recruited from the CBR-TLSA cohort. A total of 600 participants were initially enrolled in this investigation. However, data from 536 participants were only utilised for the current study, as the remaining participants’ data lacked adequate ECG waveforms or had missing cognitive scores. Among the cohort of 536 participants, 239 individuals were subsequently excluded from analysis on the basis of their concurrent use of antihypertensive medications, given the potential influence of such medications on both HR and usHRV measures which led to a final cohort size of 297 (Figure 1). The mean age of the participants in the cohort was 61.42 ± 9.75 years. Out of the 297 individuals, 48.48% (144) were male, and 51.51% (153) were female. Flowchart representing the selection of participants.
Cognitive Assessment
Addenbrooke’s Cognitive Examination (ACE-III) is used to assess cognitive function of the participants and identify the MCI cases. Generally, the ACE is used to aid in the detection and evaluation of cognitive impairment, including dementia and MCI. The ACE test is commonly used in clinical practice and research settings to assess cognitive function quickly and reliably. The ACE-III is a revised version of the original ACE test, incorporating additional subtests and enhancements to improve its sensitivity and specificity in detecting cognitive impairment. The ACE-III provides a total score as well as scores for individual cognitive domains, allowing for a comprehensive assessment of cognitive function. 21 ACE-III was administered in English as the participants belonged to an urban community and were proficient in English. For those who were not comfortable in English, a translated version of ACE-III in Kannada was administered as it is the local native language. 22
As already mentioned, in the current study, the cognitive assessment and group classification are performed using the ACE-III test, consisting of 24 diverse questions encompassing various cognitive domains, such as attention, memory, fluency, language, and visuospatial abilities. Participants scoring below 88 for the ACE-III test were classified as MCI, whereas those scoring greater than or equal to 88 were classified as healthy controls. 21 The comprehensive use of the ACE-III test facilitated the identification of individuals experiencing cognitive decline, contributing valuable insights towards the association between HR dynamics, usHRV, and cognitive impairments. Among the subjects, 19.19% of the participants (57 in number) showed MCI, leaving the remaining 80.8% (240) classified as cognitively normal. The age spectrum ranged from 45 to 92 years, with a calculated mean age of 61.42 ± 9.75 years. Specifically, the mean age among individuals with MCI was 61.75 ± 10.19 years, and that for non-MCI individuals was 61.34 ± 9.66 years.
Computation of HR and usHRV Features
Classical 12-lead ECG recordings are collected from the participants. The ECG recordings are obtained using the BPL Cardiart 9108 12-channel ECG recorder, operating at a high sampling frequency of 16 kHz. During the recording session, the participants were in a supine position. ECG data were collected for a standardized duration of 10 s. Compared to other studies, this method utilizes a short duration of ECG acquisition. From the acquired signal, the Lead II was used further as the R wave peak is more prominent and isolated in the lead II. Subsequently, preprocessing is done. In the preprocessing stage, the signal undergoes detrending. The detrending eliminates low-frequency trends or baseline drifts. After detrending, a band-pass filter (0.5 to 100 Hz) is applied to extract the spectral components of diagnostic importance and avoid radio frequency (RF) interferences. A notch filter is used in addition to the band pass filter to avoid powerline interference.
After the preprocessing stage, a custom signal processing pipeline implemented in Matlab® to compute the usHRV features was used. The process begins with the detection of prominent peaks in the lead II of the ECG signal. At each sample with index
Only the peaks with relatively higher amplitude are detected as QRS complexes. An arbitrary sample
The mean NN interval is.
In (1),
The standard deviation of the NN intervals is,
The root mean square of the successive differences between NN intervals is,
Finally, Wilcoxon rank-sum test is used to assesses the statistical significance of the features. Architecture of the proposed signal processing pipeline for computing HR & usHRV features and validation discussed above is shown in Figure 2. Architecture of the proposed signal processing pipeline for computing HR & HRV features and their validation.
Statistical Validation of HR and usHRV Features
We used Wilcoxon rank sum test to measure the extent to which the HR, and usHRV features are different among MCI and non-MCI cases. 23 The Wilcoxon rank sum test is a non-parametric hypothesis test used to compare the features extracted from two independent groups when the data do not meet the assumptions of normality required by parametric tests like the t test. It is commonly used when the data are ordinal, interval, or ratio scale but not normally distributed. The Wilcoxon rank sum test offers lower P-values significantly lesser than the level of significance when the distribution parameters of feature sets of different classes are unlike. We use wrapper methods features 24 also to assess the ‘goodness’ of HR feature (RMS). The RMS feature is assessed in terms of the accuracy offered by three ML models, namely SVM, DA and NB in distinguishing MCI, and non-MCI cases. A 10-fold cross validation strategy is adopted in the ML-based validation. 25 In the wrapper method, 70% of the feature data is assigned as the training set and 30 % as the testing set.
Results
P-Values Obtained From Wilcoxon Rank-Sum Test on HR, and HRV Features of MCI Patients and Healthy Individuals.

A bar graph illustrating the P-values of four different features: mean of NN intervals, RMS of NN intervals, SDNN, and RMSSD.
Figure 4 shows the box whisker plots demonstrating the separability between the features, mean, and RMS of the NN intervals, SDNN and RMSSD. We can clearly see the nonoverlapping blocks of non-MCI and MCI groups for mean and RMS of NN intervals. Nonoverlapping blocks, and unequal median lines of the boxes of mean and RMS of the NN intervals of MCI and non-MCI classes justify that HR features are significantly different between the two groups. Box Whisker plots showing the separability between features: (a) mean of NN Intervals, (b) RMS of NN intervals, (c) SD of NN intervals, and d) RMSSD of NN intervals.
Figure 5 shows the consecutive RR interval time series of 5 participants from MCI and non-MCI groups. For MCI patients, the RR interval is less, demonstrating that they have higher heart rate. Whereas, for non-MCI patients, the RR interval is more, indicating that the heart rate is comparatively less. It is also seen that the RR intervals of the non-MCI subjects exhibit greater variability compared to the RR interval values of the MCI patients. This indicates that the non-MCI subjects exhibit increased usHRV. RR intervals of (a) 5 non-MCI and (b) 5 MCI participants.
The k-fold validation accuracies of mean for three different machine learning algorithms, namely, SVM, DA, and NB. The RMS feature produced a very high accuracy of 80.80%, regardless of the classifier. The high accuracy values justify the potential of HR feature to be used as one of the biomarkers for detecting MCI.
Discussion
Our research highlights the potential of HR and usHRV as a biomarker for MCI detection among the elderly Indian population. Notably, our study is one of a kind on its focus on examining HR and its ultra- short variability for MCI detection within an Indian population. It is also the first study done in an elderly Indian population cohort, utilizing an ECG recording of only 10 s. Our findings align significantly with prior investigations that have explored the role of HR (Molloy et al 2023), evaluated by the measures of central tendencies of NN intervals (mean and RMS), which can act as a biomarker for cognitive decline.
The literature demonstrates a connection between HR, HRV and cognition through the rank sum test and RE measure. These two methods displayed significant values as they showed a greater dissimilarity between the features, making it easier to distinguish the MCI group from the non-MCI group. This can be interpreted from the box whisker plots.
A comprehensive systematic review exploring the correlation between HRV and MCI found a prevailing positive relationship across 19 out of 20 studies, signifying those individuals with higher HRV generally exhibited higher cognitive performance. 7 This suggests a potential link between HRV and MCI. Notably, the variance in the outcomes of our study might be attributed to the duration of the ECG signal considered; we opted for shorter durations of ECG signal, of 10 s, while most of the other studies employed a minimum of 5-minute ECG signal duration. 26 This temporal distinction could contribute to the observed disparities in our study outcomes compared to prior research.
In a study conducted by Schaich et al (2020), higher SDNN was associated with higher scores on the Cognitive Abilities Screening Instrument (CASI) and Digit Symbol Coding test (DSC) with P = .018 and P = .013 respectively. But there was no connection between HRV change and any of the cognitive scores. 27 Here, we are showing a result which suggests that MCI has an influence on HR and usHRV. The mean and RMS values in non-MCI are very high as seen from the box whisker plots. This signifies that the HR in such individuals without MCI will be low, whereas the HR in MCI patients will be high. Our observations are in accordance with the results of the study conducted by Imahori et al (2022). 14 This proves that HR is a major criterion that will distinguish individuals with MCI from those without MCI. Similarly, there is a greater usHRV in patients with MCI than those without MCI. So, both HR and usHRV can be used as biomarkers for the detection of MCI.
Based on the P-value obtained from the rank-sum test, RMS of the RR intervals is the most significant feature. In general, the detection of MCI from a single feature is a thresholding-based binary classification problem. Hence, we used a simple thresholding binary classifier that operates on the RMS feature to detect the MCI. For any typical thresholding operation, selection of the apt threshold is the most important step. Percentage classification error vs threshold graph of the binary classifier is shown in Figure 6A. In Figure 6A, minimum percentage error of 18.18% (at 81.82% Accuracy) is observed at a threshold of 0.629933. It is evident in the boxplots in Figure that the MCI patients have RMS less than the threshold. Percentage classification error vs threshold graph of the binary classifier after excluding tachycardia and bradycardia cases is shown Figure 6B. In Figure 6B minimum percentage error of 20.46% (at 79.54% Accuracy) is observed at a threshold of 0.0629971. Despite, the higher overall accuracy, and specificity, the binary classifier does not offer appreciable sensitivity. The class imbalance is also contributing to this uncertainty. Integration of more features to the classification framework may help improving the sensitivity also. There are certain subjects with their heart rate are above 100 and below 60. Exclusion of the tachycardia and bradycardia cases causes an increase in the percentage classification error offered by the binary classifier from 18.18% to 20.46%. Percentage classification error versus threshold under two different contexts. (a) Including tachycardia and bradycardia cases, (b) excluding tachycardia and bradycardia cases.
The major strengths of this study are: Here, we have used the ACE-III Cognitive Battery which is culturally validated in the Indian population to detect MCI. Demented patients were not included in this study; so, we demonstrated the ability of HR and usHRV features to detect MCI in its pre-clinical stage itself. The use of short-term ECG recordings makes it convenient for clinical utility and point-of-care applications. In addition to the usual hypothesis tests, we used ML models also, to validate the HR and usHRV features.
The current study did not include the transform-domain features like the power-ratio were not included since, we used ECG signals of short duration. The future work can focus on computing the transform-domain features without compromising the frequency resolution and signal quality. Also, combining the clinical parameters along with HR and usHRV features could improve the diagnostic accuracy of MCI.
Conclusion
We investigated the effect of MCI on HR and usHRV in the Indian population, using very short ECG recordings of individuals with MCI and those without MCI, for clinical utility and point-of-care applications. We designed a signal processing pipeline to detect the R-peaks from the pre-processed ECG signal and extract measures of central tendency (Mean and RMS) and dispersion (SDNN and RMSSD) of beat-to-beat intervals.
We obtained very small P-values for the mean, RMS, SDNN and RMSSD, which are very less than the level of significance, conveying the effect of MCI on HR and usHRV. This makes it easier for the identification of features that distinguish MCI class from non-MCI class. People with MCI have increased HR and lower usHRV when compared with healthy individuals. We assessed the discriminating capability of RMS feature with the help of various machine learning classifiers namely, SVM, DA and NB. The k-fold cross-validation accuracies were high for all three classifiers, proving the significance for MCI detection. The suitability of the proposed method for the early detection of MCI and point-of-care applications along with the use of ACE-III battery are some of the major advantages of our study.
Footnotes
Acknowledgments
We acknowledge all the participants for their valuable time and consent for data collection. We also thank TATA trusts for financial support. Last but not the least we thank our director for his unwavering support in conducting research.
Declaration of Conflicting Interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: this work was supported by the TATA Trusts.
Ethical Statement
Data Availability Statement
The datasets generated and/or analysed during the current study are not publicly available as the study is a longitudinal cohort study and is currently ongoing, the data is still being collected and curated and being monitored by the Institutional Ethics Committee (IEC) and Technical Advisory Committee (TAC). Therefore, it is not made public at this point of time. Data requests can be directed to the corresponding author Dr Thomas Gregor Issac who is the PI of the TLSA study and data will be shared if approved by the IEC and TAC.
