Every brain is unique, having its structural and functional organization shaped by both genetic and environmental factors over the course of its development. Brain image studies tend to produce ...results by averaging across a group of subjects, under the common assumption that it is possible to subdivide the cortex into homogeneous areas while maintaining a correspondence across subjects. We investigate this assumption: can the structural properties of a specific region of an atlas be assumed to be the same across subjects? This question is addressed by looking at the network representation of the brain, with nodes corresponding to brain regions and edges to their structural relationships. Using an unsupervised graph matching strategy, we align the structural connectomes of a set of healthy subjects, considering parcellations of different granularity, to understand the connectivity misalignment between regions. First, we compare the obtained permutations with four different algorithm initializations: Spatial Adjacency, Identity, Barycenter, and Random. Our results suggest that applying an alignment strategy improves the similarity across subjects when the number of parcels is above 100 and when using Spatial Adjacency and Identity initialization (the most plausible priors). Second, we characterize the obtained permutations, revealing that the majority of permutations happens between neighbors parcels. Lastly, we study the spatial distribution of the permutations. By visualizing the results on the cortex, we observe no clear spatial patterns on the permutations and all the regions across the context are mostly permuted with first and second order neighbors.
Can the structural properties of a specific region of an atlas be assumed to be the same across subjects? Using an unsupervised graph matching strategy, we align the structural connectomes of a set of healthy subjects and with parcellations of different granularity to understand which is the connectivity misalignment between regions.
Muscle electrophysiology has emerged as a powerful tool to drive human machine interfaces, with many new recent applications outside the traditional clinical domains, such as robotics and virtual ...reality. However, more sophisticated, functional, and robust decoding algorithms are required to meet the fine control requirements of these applications. Deep learning has shown high potential in meeting these demands, but requires a large amount of high-quality annotated data, which is expensive and time-consuming to acquire. Data augmentation using simulations, a strategy applied in other deep learning applications, has never been attempted in electromyography due to the absence of computationally efficient models. We introduce a concept of Myoelectric Digital Twin - highly realistic and fast computational model tailored for the training of deep learning algorithms. It enables simulation of arbitrary large and perfectly annotated datasets of realistic electromyography signals, allowing new approaches to muscular signal decoding, accelerating the development of human-machine interfaces.
Neuromechanical studies investigate how the nervous system interacts with the musculoskeletal (MSK) system to generate volitional movements. Such studies have been supported by simulation models that ...provide insights into variables that cannot be measured experimentally and allow a large number of conditions to be tested before the experimental analysis. However, current simulation models of electromyography (EMG), a core physiological signal in neuromechanical analyses, remain either limited in accuracy and conditions or are computationally heavy to apply. Here, we provide a computational platform to enable future work to overcome these limitations by presenting NeuroMotion, an open-source simulator that can modularly test a variety of approaches to the full-spectrum synthesis of EMG signals during voluntary movements. We demonstrate NeuroMotion using three sample modules. The first module is an upper-limb MSK model with OpenSim API to estimate the muscle fibre lengths and muscle activations during movements. The second module is BioMime, a deep neural network-based EMG generator that receives nonstationary physiological parameter inputs, like the afore-estimated muscle fibre lengths, and efficiently outputs motor unit action potentials (MUAPs). The third module is a motor unit pool model that transforms the muscle activations into discharge timings of motor units. The discharge timings are convolved with the output of BioMime to simulate EMG signals during the movement. We first show how MUAP waveforms change during different levels of physiological parameter variations and different movements. We then show that the synthetic EMG signals during two-degree-of-freedom hand and wrist movements can be used to augment experimental data for regressing joint angles. Ridge regressors trained on the synthetic dataset were directly used to predict joint angles from experimental data. In this way, NeuroMotion was able to generate full-spectrum EMG for the first use-case of human forearm electrophysiology during voluntary hand, wrist, and forearm movements. All intermediate variables are available, which allows the user to study cause-effect relationships in the complex neuromechanical system, fast iterate algorithms before collecting experimental data, and validate algorithms that estimate non-measurable parameters in experiments. We expect this modular platform will enable validation of generative EMG models, complement experimental approaches and empower neuromechanical research.
•A new computational framework for the analytical generation of a complete set of algebraically independent Rotation Invariant Features (RIF) for spherical functions.•These new invariants can be ...linked to statistical and geometrical measures of spherical functions e.g.: the mean, the variance and the volume.•We apply our new RIF to diffusion MRI in particular to the Apparent Diffusion Coefficient, the fiber Orientation Distribution Function, and the diffusion signal itself.•Using both synthetic and real data, we assess the sensitivity of our invariants to brain tissue microstructure related changes in the diffusion signal.
Display omitted
In this work, we present a novel computational framework for analytically generating a complete set of algebraically independent Rotation Invariant Features (RIF) given the Laplace-series expansion of a spherical function. Our computational framework provides a closed-form solution for these new invariants, which are the natural expansion of the well known spherical mean, power-spectrum and bispectrum invariants. We highlight the maximal number of algebraically independent invariants which can be obtained from a truncated Spherical Harmonic (SH) representation of a spherical function and show that most of these new invariants can be linked to statistical and geometrical measures of spherical functions, such as the mean, the variance and the volume of the spherical signal. Moreover, we demonstrate their application to dMRI signal modeling including the Apparent Diffusion Coefficient (ADC), the diffusion signal and the fiber Orientation Distribution Function (fODF). In addition, using both synthetic and real data, we test the ability of our invariants to estimate brain tissue microstructure in healthy subjects and show that our framework provides more flexibility and open up new opportunities for innovative development in the domain of microstructure recovery from diffusion MRI.
•Structure-function mappings based on eigenmodes are unified in a general framework.•Two new mappings are proposed and their performance is compared to existing mappings.•Recently published results ...are reproduced on 50 subjects of the HCP.•A glass ceiling on the prediction of the functional connectivity is obtained.
Display omitted
Characterizing the connection between brain structure and brain function is essential for understanding how behaviour emerges from the underlying anatomy. A number of studies have shown that the network structure of the white matter shapes functional connectivity. Therefore, it should be possible to predict, at least partially, functional connectivity given the structural network. Many structure–function mappings have been proposed in the literature, including several direct mappings between the structural and functional connectivity matrices. However, the current literature is fragmented and does not provide a uniform treatment of current methods based on eigendecompositions. In particular, existing methods have never been compared to each other and their relationship explicitly derived in the context of brain structure–function mapping. In this work, we propose a unified computational framework that generalizes recently proposed structure–function mappings based on eigenmodes. Using this unified framework, we highlight the link between existing models and show how they can be obtained by specific choices of the parameters of our framework. By applying our framework to 50 subjects of the Human Connectome Project, we reproduce 6 recently published results, devise two new models and provide a direct comparison between all mappings. Finally, we show that a glass ceiling on the performance of mappings based on eigenmodes seems to be reached and conclude with possible approaches to break this performance limit.
The state of the art in sampling theory now contains several theorems for signals that are non-bandlimited. For signals on the sphere however, most theorems still require the assumptions of ...bandlimitedness. In this work we show that a particular class of non-bandlimited signals, which have a finite rate of innovation, can be exactly recovered using a finite number of samples. We consider a sampling scheme where K weighted Diracs are convolved with a kernel on the rotation group. We prove that if the sampling kernel has a bandlimit L = 2K then (2K - 1)(4K - 1) + 1 equiangular samples are sufficient for exact reconstruction. If the samples are uniformly distributed on the sphere, we argue that the signal can be accurately reconstructed using 4K 2 samples and validate our claim through numerical simulations. To further reduce the number of samples required, we design an optimal sampling kernel that achieves accurate reconstruction of the signal using only 3K samples, the number of parameters of the weighted Diracs. In addition to weighted Diracs, we show that our results can be extended to sample Diracs integrated along the azimuth. Finally, we consider kernels with antipodal symmetry which are common in applications such as diffusion magnetic resonance imaging.
High angular resolution diffusion imaging (HARDI)-based tractography has been increasingly used in longitudinal studies on white matter macro- and micro-structural changes in the language network ...during language acquisition and in language impairments. However, test-retest reliability measurements are essential to ascertain that the longitudinal variations observed are not related to data processing. The aims of this study were to determine the reproducibility of the reconstruction of major white matter fiber bundles of the language network using anatomically constrained probabilistic tractography with constrained spherical deconvolution based on HARDI data, as well as to assess the test-retest reliability of diffusion measures extracted along them. Eighteen right-handed participants were scanned twice, one week apart. The arcuate, inferior longitudinal, inferior fronto-occipital, and uncinate fasciculi were reconstructed in the left and right hemispheres and the following diffusion measures were extracted along each tract: fractional anisotropy, mean, axial, and radial diffusivity, number of fiber orientations, mean length of streamlines, and volume. All fiber bundles showed good morphological overlap between the two scanning timepoints and the test-retest reliability of all diffusion measures in most fiber bundles was good to excellent. We thus propose a fairly simple, but robust, HARDI-based tractography pipeline reliable for the longitudinal study of white matter language fiber bundles, which increases its potential applicability to research on the neurobiological mechanisms supporting language.
Numerical models of electromyography (EMG) signals have provided a huge contribution to our fundamental understanding of human neurophysiology and remain a central pillar of motor neuroscience and ...the development of human-machine interfaces. However, while modern biophysical simulations based on finite element methods (FEMs) are highly accurate, they are extremely computationally expensive and thus are generally limited to modeling static systems such as isometrically contracting limbs. As a solution to this problem, we propose to use a conditional generative model to mimic the output of an advanced numerical model. To this end, we present BioMime, a conditional generative neural network trained adversarially to generate motor unit (MU) activation potential waveforms under a wide variety of volume conductor parameters. We demonstrate the ability of such a model to predictively interpolate between a much smaller number of numerical model's outputs with a high accuracy. Consequently, the computational load is dramatically reduced, which allows the rapid simulation of EMG signals during truly dynamic and naturalistic movements.
The interactions between different brain regions can be modeled as a graph,
called connectome, whose nodes correspond to parcels from a predefined brain
atlas. The edges of the graph encode the ...strength of the axonal connectivity
between regions of the atlas that can be estimated via diffusion magnetic
resonance imaging (MRI) tractography. Herein, we aim to provide a novel
perspective on the problem of choosing a suitable atlas for structural
connectivity studies by assessing how robustly an atlas captures the network
topology across different subjects in a homogeneous cohort. We measure this
robustness by assessing the alignability of the connectomes, namely the
possibility to retrieve graph matchings that provide highly similar graphs. We
introduce two novel concepts. First, the graph Jaccard index (GJI), a graph
similarity measure based on the well-established Jaccard index between sets; the
GJI exhibits natural mathematical properties that are not satisfied by previous
approaches. Second, we devise WL-align, a new technique for aligning connectomes
obtained by adapting the Weisfeiler-Leman (WL) graph-isomorphism test. We
validated the GJI and WL-align on data from the Human Connectome Project
database, inferring a strategy for choosing a suitable parcellation for
structural connectivity studies. Code and data are publicly available.
An important part of our current understanding of the structure of the human
brain relies on the concept of brain network, which is obtained by looking at
how different brain regions are connected with each other. In this paper we
present a strategy for choosing a suitable parcellation of the brain for
structural connectivity studies by making use of the concepts of network
alignment and similarity. To do so, we design a novel similarity measure between
weighted networks called graph Jaccard index, and a new network alignment
technique called WL-align. By assessing the possibility to retrieve graph
matchings that provide highly similar graphs, we show that morphology- and
structure-based atlases define brain networks that are more topologically robust
across a wide range of resolutions.