Decreased complexity of glucose dynamics preceding the onset of diabetes in mice and rats

RespirAnalyzer: an R package for analyzing data from continuous monitoring of respiratory signals

Motivation

The analysis of data obtained from continuous monitoring of respiratory signals (CMRS) holds significant importance in improving patient care, optimizing sports performance, and advancing scientific understanding in the field of respiratory health.

Results

The R package RespirAnalyzer provides an analytic tool specifically for feature extraction, fractal and complexity analysis for CMRS data. The package covers a wide and comprehensive range of data analysis methods including obtaining inter-breath intervals (IBI) series, plotting time series, obtaining summary statistics of IBI series, conducting power spectral density, multifractal detrended fluctuation analysis (MFDFA) and multiscale sample entropy analysis, fitting the MFDFA results with the extended binomial multifractal model, displaying results using various plots, etc. This package has been developed from our work in directly analyzing CMRS data and is anticipated to assist fellow researchers in computing the related features of their CMRS data, enabling them to delve into the clinical significance inherent in these features.

Availability and implementation

The package for Windows is available from both Comprehensive R Archive Network (CRAN): https://cran.r-project.org/web/packages/RespirAnalyzer/index.html and GitHub: https://github.com/dongxinzheng/RespirAnalyzer.

Decreasing complexity of glucose time series derived from continuous glucose monitoring is correlated with deteriorating glucose regulation

Most information used to evaluate diabetic statuses is collected at a special time-point, such as taking fasting plasma glucose test and providing a limited view of individual's health and disease risk. As a new parameter for continuously evaluating personal clinical statuses, the newly developed technique "continuous glucose monitoring" (CGM) can characterize glucose dynamics. By calculating the complexity of glucose time series index (CGI) with refined composite multi-scale entropy analysis of the CGM data, the study showed for the first time that the complexity of glucose time series in subjects decreased gradually from normal glucose tolerance to impaired glucose regulation and then to type 2 diabetes (P for trend < 0.01). Furthermore, CGI was significantly associated with various parameters such as insulin sensitivity/secretion (all P < 0.01), and multiple linear stepwise regression showed that the disposition index, which reflects β-cell function after adjusting for insulin sensitivity, was the only independent factor correlated with CGI (P < 0.01). Our findings indicate that the CGI derived from the CGM data may serve as a novel marker to evaluate glucose homeostasis.

Dynamic properties of glucose complexity during the course of critical illness: a pilot study

Methods to control the blood glucose (BG) levels of patients in intensive care units (ICU) improve the outcomes. The development of continuous BG levels monitoring devices has also permitted to optimize these processes. Recently it was shown that a complexity loss of the BG signal is linked to poor clinical outcomes. Thus, it becomes essential to decipher this relation to design efficient BG level control methods. In previous studies the BG signal complexity was calculated as a single index for the whole ICU stay. Although, these approaches did not grasp the potential variability of the BG signal complexity. Therefore, we setup this pilot study using a continuous monitoring of central venous BG levels in ten critically ill patients (EIRUS platform, Maquet Critical CARE AB, Solna, Sweden). Data were processed and the complexity was assessed by the detrended fluctuation analysis and multiscale entropy (MSE) methods. Finally, recordings were split into 24 h overlapping intervals and a MSE analysis was applied to each of them. The MSE analysis on time intervals revealed an entropy variation and allowed periodic BG signal complexity assessments. To highlight differences of MSE between each time interval we calculated the MSE complexity index defined as the area under the curve. This new approach could pave the way to future studies exploring new strategies aimed at restoring blood glucose complexity during the ICU stay.

A comprehensive comparison and overview of R packages for calculating sample entropy

Sample entropy is a powerful tool for analyzing the complexity and irregularity of physiology signals which may be associated with human health. Nevertheless, the sophistication of its calculation hinders its universal application. As of today, the R language provides multiple open-source packages for calculating sample entropy. All of which, however, are designed for different scenarios. Therefore, when searching for a proper package, the investigators would be confused on the parameter setting and selection of algorithms. To ease their selection, we have explored the functions of five existing R packages for calculating sample entropy and have compared their computing capability in several dimensions. We used four published datasets on respiratory and heart rate to study their input parameters, types of entropy, and program running time. In summary, NonlinearTseries and CGManalyzer can provide the analysis of sample entropy with different embedding dimensions and similarity thresholds. CGManalyzer is a good choice for calculating multiscale sample entropy of physiological signal because it not only shows sample entropy of all scales simultaneously but also provides various visualization plots. MSMVSampEn is the only package that can calculate multivariate multiscale entropies. In terms of computing time, NonlinearTseries, CGManalyzer, and MSMVSampEn run significantly faster than the other two packages. Moreover, we identify the issues in MVMSampEn package. This article provides guidelines for researchers to find a suitable R package for their analysis and applications using sample entropy.

Analyzing Complexity and Fractality of Glucose Dynamics in a Pregnant Woman with Type 2 Diabetes under Treatment

Currently, the rapid development of continuous glucose monitoring (CGM) device brings new insights into the treatment of diabetic patients including those during pregnancy. Complexity and fractality have recently under fast development for extracting information embodied in glucose dynamics measured using CGM. Although scientists have investigated the difference of complexity in glucose dynamics between diabetes and non-diabetes in order to discover better approaches for diabetes care, no one has analyzed the complexity and fractality of glucose dynamics during the process of adopting CGM to successfully treat pregnant women with type 2 diabetes. Thus, we analyzed the complexity and fractality using power spectral density (PSD), multi-scale sample entropy (MSE) and multifractal detrended fluctuation analysis (MF-DFA) in a clinical case. Our results show that (i) there exists multifractal behavior in blood glucose dynamics; (ii) the alpha stable distribution fits to the glucose increment data better than the Gaussian distribution; and (iii) the "global" complexity indicated by multiscale entropy, spectrum exponent and Hurst exponent increase and the "local" complexity indicated by multifractal spectrum decrease after the successful therapy. Our results offer findings that may bring value to health care providers for managing glucose levels of pregnant women with type 2 diabetes as well as provide scientists a reference on applying complexity and fractality in the clinical practice of treating diabetes.

An Improved Method of Handling Missing Values in the Analysis of Sample Entropy for Continuous Monitoring of Physiological Signals

Medical devices generate huge amounts of continuous time series data. However, missing values commonly found in these data can prevent us from directly using analytic methods such as sample entropy to reveal the information contained in these data. To minimize the influence of missing points on the calculation of sample entropy, we propose a new method to handle missing values in continuous time series data. We use both experimental and simulated datasets to compare the performance (in percentage error) of our proposed method with three currently used methods: skipping the missing values, linear interpolation, and bootstrapping. Unlike the methods that involve modifying the input data, our method modifies the calculation process. This keeps the data unchanged which is less intrusive to the structure of the data. The results demonstrate that our method has a consistent lower average percentage error than other three commonly used methods in multiple common physiological signals. For missing values in common physiological signal type, different data size and generating mechanism, our method can more accurately extract the information contained in continuously monitored data than traditional methods. So it may serve as an effective tool for handling missing values and may have broad utility in analyzing sample entropy for common physiological signals. This could help develop new tools for disease diagnosis and evaluation of treatment effects.

Decreased complexity of glucose dynamics in diabetes in rhesus monkeys

Until recently, preclinical and clinical work on diabetes has focused on the understanding of blood glucose elevation and its detrimental metabolic sequelae. The advent of continuous glucose monitoring (CGM) technology now allows real time monitoring of blood glucose levels as a time series, and thus the exploration of glucose dynamics at short time scales. Previous work has shown decreases in the complexity of glucose dynamics, as measured by multiscale entropy (MSE) analysis, in diabetes in humans, mice, and rats. Analyses for non-human primates (NHP) have not been reported, nor is it known if anti-diabetes compounds affect complexity of glucose dynamics. We instrumented four healthy and six diabetic rhesus monkeys with CGM probes in the carotid artery and collected glucose values at a frequency of one data point per second for the duration of the sensors' life span. Sensors lasted between 45 and 78 days. Five of the diabetic rhesus monkeys were also administered the anti-diabetic drug liraglutide daily beginning at day 39 of the CGM monitoring period. Glucose levels fluctuated during the day in both healthy and diabetic rhesus monkeys, peaking between 12 noon - 6 pm. MSE analysis showed reduced complexity of glucose dynamics in diabetic monkeys compared to healthy animals. Although liraglutide decreased glucose levels, it did not restore complexity in diabetic monkeys consistently. Complexity varied by time of day, more strongly for healthy animals than for diabetic animals. And by dividing the monitoring period into 3-day or 1-week subperiods, we were able to estimate within-animal variability of MSE curves. Our data reveal that decreased complexity of glucose dynamics is a conserved feature of diabetes from rodents to NHPs to man.

The Difference δ-Cells Make in Glucose Control

The role of beta and α-cells to glucose control are established, but the physiological role of δ-cells is poorly understood. Delta-cells are ideally positioned within pancreatic islets to modulate insulin and glucagon secretion at their source. We review the evidence for a negative feedback loop between delta and β-cells that determines the blood glucose set point and suggest that local δ-cell-mediated feedback stabilizes glycemic control.

Heart Rate Fragmentation as a Novel Biomarker of Adverse Cardiovascular Events: The Multi-Ethnic Study of Atherosclerosis

Background: A major objective of precision medicine is the elucidation of non-invasive biomarkers of cardiovascular (CV) risk. Recently, we introduced a new dynamical marker of sino-atrial instability, termed heart rate fragmentation (HRF), which outperformed traditional and nonlinear heart rate variability metrics in separating ostensibly healthy subjects from patients with coronary artery disease. Accordingly, we hypothesized that HRF may be a dynamical biomarker of adverse cardiovascular events (CVEs). Methods: This study employed data from a cohort of participants in the Multi-Ethnic Study of Atherosclerosis (MESA), a prospective study of sub-clinical heart disease. Interbeat interval time series (n = 1963), derived from the electrocardiographic channel of the polysomnogram study, were analyzed using the newly introduced metrics of fragmentation, as well as traditional heart rate variability (HRV) indices and the short-term detrended fluctuation analysis exponent. Cox regression analysis was used to assess the association between HR dynamic indices and CV outcomes in unadjusted and adjusted models. Results: The mean (± SD) follow-up time was 2.97 ± 0.63 years. In adjusted models, higher fragmentation was significantly associated with incident CVEs (number of events; hazard ratio [95% confidence interval]: n = 72, 1.43 [1.16-1.76]) and CV death (n = 21; 1.65 [1.15-2.36]). The traditional HRV and the fractal indices were not associated with CVEs or CV death. The most discriminatory fragmentation indices added significant value to Framingham and MESA CV risk indices in all analyses. Conclusion: Our findings show that HRF has promise as a non-invasive, automatable biomarker of CV risk. The basic mechanisms underlying fragmentation remain to be delineated. Its association with incident outcomes raises the possibility of connections to degenerative changes in the multisystem network controlling SAN function.

Demonstrating the Potential of Using Transcutaneous Oxygen and Carbon Dioxide Tensions to Assess the Risk of Pressure Injuries

Pressure injuries have a high incidence in elderly and critically ill patients, and can endanger lives in severe cases. The key to reducing the incidence of pressure injuries is to find an objective, noninvasive, automatic and consistent scientific method for assessing pressure injuries. To serve this need, we conducted a clinical study to investigate the potential of using transcutaneous oxygen tension (TcPO2) and transcutaneous carbon dioxide tension (TcPCO2) for assessing pressure injuries. From the results of the study we found that first, the values of TcPO2 and TcPCO2 are sensitive to the change of pressure imposed on the measured region and to the risk status of a pressure injury when a pressure is imposed. Second, the magnitude of change in TcPO2 and TcPCO2 is higher in patients with a high risk of a pressure injury compared with those who have a low risk. Third, TcPO2 and TcPCO2 are both significantly correlated with the Braden score, the widely used score for assessing the risk of a pressure injury. Therefore, TcPO2 and TcPCO2 have a potential to be an effective and convenient scientific tool for assessing the risk of pressure injuries.

Glucose Variability: A Review of Clinical Applications and Research Developments

Glycemic variability (GV) is a major consideration when evaluating quality of glycemic control. GV increases progressively from prediabetes through advanced T2D and is still higher in T1D. GV is correlated with risk of hypoglycemia. The most popular metrics for GV are the %Coefficient of Variation (%CV) and standard deviation (SD). The %CV is correlated with risk of hypoglycemia. Graphical display of glucose by date, time of day, and day of the week, and display of simplified glucose distributions showing % of time in several ranges, provide clinically useful indicators of GV. SD is highly correlated with most other measures of GV, including interquartile range, mean amplitude of glycemic excursion, mean of daily differences, and average daily risk range. Some metrics are sensitive to the frequency, periodicity, and complexity of glycemic fluctuations, including Fourier analysis, periodograms, frequency spectrum, multiscale entropy (MSE), and Glucose Variability Percentage (GVP). Fourier analysis indicates progressive changes from normal subjects to children and adults with T1D, and from prediabetes to T2D. The GVP identifies novel characteristics for children, adolescents, and adults with type 1 diabetes and for adults with type 2. GVP also demonstrated small rapid glycemic fluctuations in people with T1D when using a dual-hormone closed-loop control. MSE demonstrated systematic changes from normal subjects to people with T2D at various stages of duration, intensity of therapy, and quality of glycemic control. We describe new metrics to characterize postprandial excursions, day-to-day stability of glucose patterns, and systematic changes of patterns by day of the week. Metrics for GV should be interpreted in terms of percentiles and z-scores relative to identified reference populations. There is a need for large accessible databases for reference populations to provide a basis for automated interpretation of GV and other features of continuous glucose monitoring records.

Detection of sputum by interpreting the time-frequency distribution of respiratory sound signal using image processing techniques

Motivation

Sputum in the trachea is hard to expectorate and detect directly for the patients who are unconscious, especially those in Intensive Care Unit. Medical staff should always check the condition of sputum in the trachea. This is time-consuming and the necessary skills are difficult to acquire. Currently, there are few automatic approaches to serve as alternatives to this manual approach.

Results

We develop an automatic approach to diagnose the condition of the sputum. Our approach utilizes a system involving a medical device and quantitative analytic methods. In this approach, the time-frequency distribution of respiratory sound signals, determined from the spectrum, is treated as an image. The sputum detection is performed by interpreting the patterns in the image through the procedure of preprocessing and feature extraction. In this study, 272 respiratory sound samples (145 sputum sound and 127 non-sputum sound samples) are collected from 12 patients. We apply the method of leave-one out cross-validation to the 12 patients to assess the performance of our approach. That is, out of the 12 patients, 11 are randomly selected and their sound samples are used to predict the sound samples in the remaining one patient. The results show that our automatic approach can classify the sputum condition at an accuracy rate of 83.5%.

Availability and implementation

The matlab codes and examples of datasets explored in this work are available at Bioinformatics online.

Contact

yesoyou@gmail.com or douglaszhang@umac.mo.

Supplementary information

Supplementary data are available at Bioinformatics online.