Functional brain-heart interplay extends to the multifractal domain.

The study of functional brain-heart interplay has provided meaningful insights in cardiology and neuroscience. Regarding biosignal processing, this interplay involves predominantly neural and heartbeat linear dynamics expressed via time and frequency domain-related features. However, the dynamics of central and autonomous nervous systems show nonlinear and multifractal behaviours, and the extent to which this behaviour influences brain-heart interactions is currently unknown. Here, we report a novel signal processing framework aimed at quantifying nonlinear functional brain-heart interplay in the non-Gaussian and multifractal domains that combines electroencephalography (EEG) and heart rate variability series. This framework relies on a maximal information coefficient analysis between nonlinear multiscale features derived from EEG spectra and from an inhomogeneous point-process model for heartbeat dynamics. Experimental results were gathered from 24 healthy volunteers during a resting state and a cold pressor test, revealing that synchronous changes between brain and heartbeat multifractal spectra occur at higher EEG frequency bands and through nonlinear/complex cardiovascular control. We conclude that significant bodily, sympathovagal changes such as those elicited by cold-pressure stimuli affect the functional brain-heart interplay beyond second-order statistics, thus extending it to multifractal dynamics. These results provide a platform to define novel nervous-system-targeted biomarkers. This article is part of the theme issue 'Advanced computation in cardiovascular physiology: new challenges and opportunities'.

Drowsiness detection from polysomnographic data using multivariate selfsimilarity and eigen-wavelet analysis.

Because drowsiness is a major cause in vehicle accidents, its automated detection is critical. Scale-free temporal dynamics is known to be typical of physiological and body rhythms. The present work quantifies the benefits of applying a recent and original multivariate selfsimilarity analysis to several modalities of polysomnographic measurements (heart rate, blood pressure, electroencephalogram and respiration), from the MIT-BIH Polysomnographic Database, to better classify drowsiness-related sleep stages. Clinical relevance- This study shows that probing jointly temporal dynamics amongst polysomnographic measurements, with a proposed original multivariate multiscale approach, yields a gain of above 5% in the Area-under-Curve quanti-fying drowsiness-related sleep stage classification performance compared to univariate analysis.

Learning Event Representations for Temporal Segmentation of Image Sequences by Dynamic Graph Embedding.

Recently, self-supervised learning has proved to be effective to learn representations of events suitable for temporal segmentation in image sequences, where events are understood as sets of temporally adjacent images that are semantically perceived as a whole. However, although this approach does not require expensive manual annotations, it is data hungry and suffers from domain adaptation problems. As an alternative, in this work, we propose a novel approach for learning event representations named Dynamic Graph Embedding (DGE). The assumption underlying our model is that a sequence of images can be represented by a graph that encodes both semantic and temporal similarity. The key novelty of DGE is to learn jointly the graph and its graph embedding. At its core, DGE works by iterating over two steps: 1) updating the graph representing the semantic and temporal similarity of the data based on the current data representation, and 2) updating the data representation to take into account the current data graph structure. The main advantage of DGE over state-of-the-art self-supervised approaches is that it does not require any training set, but instead learns iteratively from the data itself a low-dimensional embedding that reflects their temporal and semantic similarity. Experimental results on two benchmark datasets of real image sequences captured at regular time intervals demonstrate that the proposed DGE leads to event representations effective for temporal segmentation. In particular, it achieves robust temporal segmentation on the EDUBSeg and EDUBSeg-Desc benchmark datasets, outperforming the state of the art. Additional experiments on two Human Motion Segmentation benchmark datasets demonstrate the generalization capabilities of the proposed DGE.

Quantifying Functional Links between Brain and Heartbeat Dynamics in the Multifractal Domain: a Preliminary Analysis.

Quantification of brain-heart interplay (BHI) has mainly been performed in the time and frequency domains. However, such functional interactions are likely to involve nonlinear dynamics associated with the two systems. To this extent, in this preliminary study we investigate the functional coupling between multifractal properties of Electroencephalography (EEG) and Heart Rate Variability (HRV) series using a channel- and time scale-wise maximal information coefficient analysis. Experimental results were gathered from 24 healthy volunteers undergoing a resting state and a cold-pressure test, and suggest that significant changes between the two experimental conditions might be associated with nonlinear quantifiers of the multifractal spectrum. Particularly, major brain-heart functional coupling was associated with the secondorder cumulant of the multifractal spectrum. We conclude that a functional nonlinear relationship between brain- and heartbeat-related multifractal sprectra exist, with higher values associated with the resting state.

Revisiting Functional Connectivity for Infraslow Scale-Free Brain Dynamics Using Complex Wavelets.

The analysis of human brain functional networks is achieved by computing functional connectivity indices reflecting phase coupling and interactions between remote brain regions. In magneto- and electroencephalography, the most frequently used functional connectivity indices are constructed based on Fourier-based cross-spectral estimation applied to specific fast and band-limited oscillatory regimes. Recently, infraslow arrhythmic fluctuations (below the 1 Hz) were recognized as playing a leading role in spontaneous brain activity. The present work aims to propose to assess functional connectivity from fractal dynamics, thus extending the assessment of functional connectivity to the infraslow arrhythmic or scale-free temporal dynamics of M/EEG-quantified brain activity. Instead of being based on Fourier analysis, new Imaginary Coherence and weighted Phase Lag indices are constructed from complex-wavelet representations. Their performances are first assessed on synthetic data by means of Monte-Carlo simulations, and they are then compared favorably against the classical Fourier-based indices. These new assessments of functional connectivity indices are also applied to MEG data collected on 36 individuals both at rest and during the learning of a visual motion discrimination task. They demonstrate a higher statistical sensitivity, compared to their Fourier counterparts, in capturing significant and relevant functional interactions in the infraslow regime and modulations from rest to task. Notably, the consistent overall increase in functional connectivity assessed from fractal dynamics from rest to task correlated with a change in temporal dynamics as well as with improved performance in task completion, which suggests that the complex-wavelet weighted Phase Lag index is the sole index is able to capture brain plasticity in the infraslow scale-free regime.

Multifractal formalisms for multivariate analysis.

Multifractal analysis, that quantifies the fluctuations of regularities in time series or textures, has become a standard signal/image processing tool. It has been successfully used in a large variety of applicative contexts. Yet, successes are confined to the analysis of one signal or image at a time (univariate analysis). This is because multivariate (or joint) multifractal analysis remains so far rarely used in practice and has barely been studied theoretically. In view of the myriad of modern real-world applications that rely on the joint (multivariate) analysis of collections of signals or images, univariate analysis constitutes a major limitation. The goal of the present work is to theoretically ground multivariate multifractal analysis by studying the properties and limitations of the most natural extension of the univariate formalism to a multivariate formulation. It is notably shown that while performing well for a class of model processes, this natural extension is not valid in general. Based on the theoretical study of the mechanisms leading to failure, we propose alternative formulations and examine their mathematical properties.

Wavelet p-Leader Non Gaussian Multiscale Expansions for Heart Rate Variability Analysis in Congestive Heart Failure Patients.

OBJECTIVE: Numerous indices were devised for the statistical characterization of temporal dynamics of heart rate variability (HRV) with the aim to discriminate between healthy subjects and nonhealthy patients. Elaborating on the concepts of (multi)fractal and nonlinear analyses, the present contribution defines and studies formally novel non Gaussian multiscale representations. METHODS: A methodological framework for non Gaussian multiscale representations constructed on wavelet p-leaders is developed, relying a priori neither on exact scale-free dynamics nor on predefined forms of departure from Gaussianity. Its versatility in quantifying the strength and nature of departure from Gaussian is analyzed theoretically and numerically. The ability of the representations to discriminate between healthy subjects and congestive heart failure (CHF) patients, and between survivors and nonsurvivor CHF patients, is assessed on a large cohort of 198 subjects. RESULTS: The analysis leads to conclude that i) scale-free and multifractal dynamics are observed, both for healthy subjects and CHF patients, for time scales shorter than [Formula: see text]; ii) a circadian evolution of multifractal and non Gaussian properties of HRV is evidenced for healthy subjects, but not for CHF patients; iii) non Gaussian multiscale indices possess high discriminative abilities between survivor and nonsurvivor CHF patients, at specific time scales ([Formula: see text] and [Formula: see text]). CONCLUSIONS: The non Gaussian multiscale representations provide evidence for the existence of short-term cascade-type multifractal mechanisms underlying HRV for both healthy and CHF subjects. A circadian evolution of this mechanism is only evidenced for the healthy group, suggesting an alteration of the sympathetic-parasympathetic balance for CHF patients. SIGNIFICANCE: Results obtained for a large cohort of subjects suggest that the novel non Gaussian indices might robustly quantify crucial information for clinical risk stratification in CHF patients.

Mortality Prediction in Severe Congestive Heart Failure Patients with Multifractal Point-Process Modeling of Heartbeat Dynamics.

Multifractal analysis of human heartbeat dynamics has been demonstrated to provide promising markers of Congestive Heart Failure (CHF). Yet, it crucially builds on the interpolation of RR intervals series, which has been generically performed with limited links to CHF pathophysiology. We devise a novel methodology estimating multifractal autonomic dynamics from heartbeat-derived series defined in the continuous time. We hypothesize that markers estimated from our novel framework are also effective for mortality prediction in severe CHF. We merge multifractal analysis within a methodological framework based on inhomogeneous point process models of heartbeat dynamics. Specifically, wavelet coefficients and wavelet leaders are computed over measures extracted from instantaneous statistics of probability density functions characterizing and predicting the time until the next heartbeat event occurs. The proposed approach is tested on data from 94 CHF patients, aiming at predicting survivor and non-survivor individuals as determined after a 4 years follow up. Instantaneous markers of vagal and sympatho-vagal dynamics display power-law scaling for a large range of scales, from s to s. Using standard SVM algorithms, the proposed inhomogeneous point-process representation based multifractal analysis achieved the best CHF mortality prediction accuracy of 79.11 % (sensitivity 90.48%, specificity 67.74%). Our results suggest that heartbeat scaling and multifractal properties in CHF patients are not generated at the sinus-node level, but rather by the intrinsic action of vagal short-term control and of sympatho-vagal fluctuations associated with circadian cardiovascular control, especially within the VLF band. These markers might provide critical information in devising a clinical tool for individualized prediction of survivor and non-survivor CHF patients.

Scattering Transform of Heart Rate Variability for the Prediction of Ischemic Stroke in Patients with Atrial Fibrillation.

BACKGROUND: Atrial fibrillation (AF) is an identified risk factor for ischemic strokes (IS). AF causes a loss in atrial contractile function that favors the formation of thrombi, and thus increases the risk of stroke. Also, AF produces highly irregular and complex temporal dynamics in ventricular response RR intervals. Thus, it is hypothesized that the analysis of RR dynamics could provide predictors for IS. However, these complex and nonlinear dynamics call for the use of advanced multiscale nonlinear signal processing tools. OBJECTIVES: The global aim is to investigate the performance of a recently-proposed multiscale and nonlinear signal processing tool, the scattering transform, in predicting IS for patients suffering from AF. METHODS: The heart rate of a cohort of 173 patients from Fujita Health University Hospital in Japan was analyzed with the scattering transform. First, p-values of Wilcoxon rank sum tests were used to identify scattering coefficients achieving significant (univariate) discrimination between patients with and without IS. Second, a multivariate procedure for feature selection and classification, the Sparse Support Vector Machine (S-SVM), was applied to predict IS. RESULTS: Groups of scattering coefficients, located at several time-scales, were identified as significantly higher (p-value < 0.05) in patients who developed IS than in those who did not. Though the overall predictive power of these indices remained moderate (around 60 %), it was found to be much higher when analysis was restricted to patients not taking antithrombotic treatment (around 80 %). Further, S-SVM showed that multivariate classification improves IS prediction, and also indicated that coefficients involved in classification differ for patients with and without antithrombotic treatment. CONCLUSIONS: Scattering coefficients were found to play a significant role in predicting IS, notably for patients not receiving antithrombotic treatment. S-SVM improves IS detection performance and also provides insight on which features are important. Notably, it shows that AF patients not taking antithrombotic treatment are characterized by a slow modulation of RR dynamics in the ULF range and a faster modulation in the HF range. These modulations are significantly decreased in patients with IS, and hence have a good discriminant ability.

Spatially regularized multifractal analysis for fMRI data.

Scale-free dynamics is nowadays a massively used paradigm to model infraslow macroscopic brain activity. Multifractal analysis is becoming the standard tool to characterize scale-free dynamics. It is commonly used on various modalities of neuroimaging data to evaluate whether arrhythmic fluctuations in ongoing or evoked brain activity are related to pathologies (Alzheimer, epilepsy) or task performance. The success of multifractal analysis in neurosciences remains however so far contrasted: While it lead to relevant findings on M/EEG data, less clear impact was shown when applied to fMRI data. This is mostly due to their poor time resolution and very short duration as well as to the fact that analysis remains performed voxelwise. To take advantage of the large amount of voxels recorded jointly in fMRI, the present contribution proposes the use of a recently introduced Bayesian formalism for multifractal analysis, that regularizes the estimation of the multifractality parameter of a given voxel using information from neighbor voxels. The benefits of this regularized multifractal analysis are illustrated by comparison against classical multifractal analysis on fMRI data collected on one subject, at rest and during a working memory task: Though not yet statistically significant, increased multifractality is observed in task-negative and task-positive networks, respectively.

Bayesian estimation of the multifractality parameter for image texture using a whittle approximation.

Texture characterization is a central element in many image processing applications. Multifractal analysis is a useful signal and image processing tool, yet, the accurate estimation of multifractal parameters for image texture remains a challenge. This is due in the main to the fact that current estimation procedures consist of performing linear regressions across frequency scales of the 2D dyadic wavelet transform, for which only a few such scales are computable for images. The strongly non-Gaussian nature of multifractal processes, combined with their complicated dependence structure, makes it difficult to develop suitable models for parameter estimation. Here, we propose a Bayesian procedure that addresses the difficulties in the estimation of the multifractality parameter. The originality of the procedure is threefold. The construction of a generic semiparametric statistical model for the logarithm of wavelet leaders; the formulation of Bayesian estimators that are associated with this model and the set of parameter values admitted by multifractal theory; the exploitation of a suitable Whittle approximation within the Bayesian model which enables the otherwise infeasible evaluation of the posterior distribution associated with the model. Performance is assessed numerically for several 2D multifractal processes, for several image sizes and a large range of process parameters. The procedure yields significant benefits over current benchmark estimators in terms of estimation performance and ability to discriminate between the two most commonly used classes of multifractal process models. The gains in performance are particularly pronounced for small image sizes, notably enabling for the first time the analysis of image patches as small as 64 × 64 pixels.