Sub-millimeter fMRI reveals multiple topographical digit representations that form action maps in human motor cortex

Tactile stimulation designs adapted to clinical settings result in reliable fMRI-based somatosensory digit maps

Movement constraints in stroke survivors are often accompanied by additional impairments in related somatosensory perception. A complex interplay between the primary somatosensory and motor cortices is essential for adequate and precise movements. This necessitates investigating the role of the primary somatosensory cortex in movement deficits of stroke survivors. The first step towards this goal could be a fast and reliable functional Magnetic Resonance Imaging (fMRI)-based mapping of the somatosensory cortex applicable for clinical settings. Here, we compare two 3 T fMRI-based somatosensory digit mapping techniques adapted for clinical usage in seven neurotypical volunteers and two sessions, to assess their validity and retest-reliability. Both, the traveling wave and the blocked design approach resulted in complete digit maps in both sessions of all participants, showing the expected layout. Similarly, no evidence for differences in the volume of activation, nor the activation overlap between neighboring activations could be detected, indicating the general feasibility of the clinical adaptation and their validity. Retest-reliability, indicated by the Dice coefficient, exhibited reasonable values for the spatial correspondence of single digit activations across sessions, but low values for the spatial correspondence of the area of overlap between neighboring digits across sessions. Parameters describing the location of the single digit activations exhibited very high correlations across sessions, while activation volume and overlap only exhibited medium to low correlations. The feasibility and high retest-reliabilities for the parameters describing the location of the single digit activations are promising concerning the implementation into a clinical context to supplement diagnosis and treatment stratification in upper limb stroke patients.

Artificial embodiment displaces cortical neuromagnetic somatosensory responses

Integrating artificial limbs as part of one's body involves complex neuroplastic changes resulting from various sensory inputs. While somatosensory feedback is crucial, plastic processes that enable embodiment remain unknown. We investigated this using somatosensory evoked fields (SEFs) in the primary somatosensory cortex (S1) following the Rubber Hand Illusion (RHI), known to quickly induce artificial limb embodiment. During electrical stimulation of the little finger and thumb, 19 adults underwent neuromagnetic recordings before and after the RHI. We found early SEF displacement, including an illusion-brain correlation between extent of embodiment and specific changes to the first cortical response at 20 ms in Area 3b, within S1. Furthermore, we observed a posteriorly directed displacement at 35 ms towards Area 1, known to be important for visual integration during touch perception. That this second displacement was unrelated to extent of embodiment implies a functional distinction between neuroplastic changes of these components and areas. The earlier shift in Area 3b may shape extent of limb ownership, while subsequent displacement into Area 1 may relate to early visual-tactile integration that initiates embodiment. Here we provide evidence for multiple neuroplastic processes in S1-lasting beyond the illusion-supporting integration of artificial limbs like prostheses within the body representation.

Ipsilateral stimulation shows somatotopy of thumb and shoulder auricular points on the left primary somatosensory cortex using high-density fNIRS

Significance

Auriculotherapy is a technique based on stimulation applied to specific ear points. Its mechanism of active and clinical efficacy remain to be established. This study aims to assess the role that primary somatosensory cortex may play to validate auriculotherapy mechanisms.

Aim

This study examined whether tactile stimulation at specific auricular points is correlated with distinct cortical activation in the primary somatosensory cortex.

Approach

Seventeen healthy adults participated in the study. Tactile stimuli were delivered to the thumb, shoulder, and skin master points on the ear using von Frey filaments. Functional near-infrared spectroscopy was used to measure and spatially map cortical responses.

Results

This study revealed distinct hemodynamic activity patterns in response to ear point stimulation, consistent with the classic homunculus model of somatotopic organization. Ipsilateral stimulation showed specific cortical activations for the thumb and shoulder points, while contralateral stimulation showed less significant activity. Functional near-infrared spectroscopy effectively captured localized cortical responses to ear tactile stimuli, supporting the somatotopic mapping hypothesis.

Conclusion

These findings enhance the understanding of sensory processing with auricular stimulation and supports the concepts of auricular cartography that underpins some schools of auriculotherapy practice. Future research should explore bilateral cortical mapping and the integration of other neuroimaging techniques.

Principles of cortical areas and their implications for neuroimaging

Cortical organization should constrain the study of how the brain performs behavior and cognition. A fundamental concept in cortical organization is that of arealization: that the cortex is parceled into discrete areas. In part one of this report, we review how non-human animal studies have illuminated principles of cortical arealization by revealing: (1) what defines a cortical area, (2) how cortical areas are formed, (3) how cortical areas interact with one another, and (4) what "computations" or "functions" areas perform. In part two, we discuss how these principles apply to neuroimaging research. In doing so, we highlight several examples where the commonly accepted interpretation of neuroimaging observations requires assumptions that violate the principles of arealization, including nonstationary areas that move on short time scales, large-scale gradients as organizing features, and cortical areas with singular functionality that perfectly map psychological constructs. Our belief is that principles of neurobiology should strongly guide the nature of computational explanations.

Evolution, biomechanics, and neurobiology converge to explain selective finger motor control

Humans use their fingers to perform a variety of tasks, from simple grasping to manipulating objects, to typing and playing musical instruments, a variety wider than any other species. The more sophisticated the task, the more it involves individuated finger movements, those in which one or more selected fingers perform an intended action while the motion of other digits is constrained. Here we review the neurobiology of such individuated finger movements. We consider their evolutionary origins, the extent to which finger movements are in fact individuated, and the evolved features of neuromuscular control that both enable and limit individuation. We go on to discuss other features of motor control that combine with individuation to create dexterity, the impairment of individuation by disease, and the broad extent of capabilities that individuation confers on humans. We comment on the challenges facing the development of a truly dexterous bionic hand. We conclude by identifying topics for future investigation that will advance our understanding of how neural networks interact across multiple regions of the central nervous system to create individuated movements for the skills humans use to express their cognitive activity.

Laminar CBV and BOLD response-characteristics over time and space in the human primary somatosensory cortex at 7T

Uncovering the cortical representation of the body has been at the core of human brain mapping for decades, with special attention given to the digits. In the last decade, advances in functional magnetic resonance imaging (fMRI) technologies have opened the possibility of noninvasively unraveling the 3rd dimension of digit representations in humans along cortical layers. In laminar fMRI it is common to combine the use of the highly sensitive blood oxygen level dependent (BOLD) contrast with cerebral blood volume sensitive measurements, like vascular space occupancy (VASO), that are more specific to the underlying neuronal populations. However, the spatial and temporal VASO response characteristics across cortical depth to passive stimulation of the digits are still unknown. Therefore, we characterized haemodynamic responses to vibrotactile stimulation of individual digit-tips across cortical depth at 0.75 mm in-plane spatial resolution using BOLD and VASO fMRI at 7T. We could identify digit-specific regions of interest (ROIs) in putative Brodmann area 3b, following the known anatomical organization. In the ROIs, the BOLD response increased towards the cortical surface due to the draining vein effect, while the VASO response was more shifted towards middle cortical layers, likely reflecting bottom-up input from the thalamus, as expected. Interestingly, we also found slightly negative BOLD and VASO responses for non-preferred digits in the ROIs, potentially indicating neuronal surround inhibition. Finally, we explored the temporal signal dynamics for BOLD and VASO as a function of distance from activation peaks resulting from stimulation of contralateral digits. With this analysis, we showed a triphasic response consisting of an initial peak and a subsequent negative deflection during stimulation, followed by a positive post-stimulus response in BOLD and to some extent in VASO. While similar responses were reported with invasive methods in animal models, here we demonstrate a potential neuronal excitation-inhibition mechanism in a center-surround architecture across layers in the human somatosensory cortex. Given that, unlike in animals, human experiments do not rely on anesthesia and can readily implement extensive behavioral testing, obtaining this effect in humans is an important step towards further uncovering the functional significance of the different aspects of the triphasic response.

Submillimeter fMRI Acquisition Techniques for Detection of Laminar and Columnar Level Brain Activation

Since the first demonstration in the early 1990s, functional MRI (fMRI) has emerged as one of the most powerful, noninvasive neuroimaging tools to probe brain functions. Subsequently, fMRI techniques have advanced remarkably, enabling the acquisition of functional signals with a submillimeter voxel size. This innovation has opened the possibility of investigating subcortical neural activities with respect to the cortical depths or cortical columns. For this purpose, numerous previous works have endeavored to design suitable functional contrast mechanisms and dedicated imaging techniques. Depending on the choice of the functional contrast, functional signals can be detected with high sensitivity or with improved spatial specificity to the actual activation site, and the pertaining issues have been discussed in a number of earlier works. This review paper primarily aims to provide an overview of the subcortical fMRI techniques that allow the acquisition of functional signals with a submillimeter resolution. Here, the advantages and disadvantages of the imaging techniques will be described and compared. We also summarize supplementary imaging techniques that assist in the analysis of the subcortical brain activation for more accurate mapping with reduced geometric deformation. This review suggests that there is no single universally accepted method as the gold standard for subcortical fMRI. Instead, the functional contrast and the corresponding readout imaging technique should be carefully determined depending on the purpose of the study. Due to the technical limitations of current fMRI techniques, most subcortical fMRI studies have only targeted partial brain regions. As a future prospect, the spatiotemporal resolution of fMRI will be pushed to satisfy the community's need for a deeper understanding of whole-brain functions and the underlying connectivity in order to achieve the ultimate goal of a time-resolved and layer-specific spatial scale. EVIDENCE LEVEL: 1 TECHNICAL EFFICACY: Stage 1.

Modulation of motor surround inhibition during motor tasks

Surround inhibition (SI) in the motor system is important in individuation of actions, but is sometimes difficult to demonstrate. It has also not been evaluated in real life tasks. In this study, we use real life tasks and a new method where excitability of the surround muscle is assessed with respect to its current activity level rather than when it is at rest. Motor evoked potential (MEP) amplitudes were measured in the abductor digiti minimi (ADM) muscle while participants performed several motor tasks: "writing" on paper, "holding a pen" precisely and, "holding a water bottle" against gravity. These MEPs were compared to ADM MEPs amplitudes measured during a fifth finger abduction (ADM being the center muscle). SI was also measured in the traditional way, by comparing ADM MEPs during an index finger flexion and at rest. For the "writing" and "holding a pen" tasks, but not the "holding bottle" task, the MEP amplitudes were significantly smaller when compared to MEP amplitudes when the ADM was the center muscle with the same level of activation. The ADM MEP amplitudes were not different between rest and during index finger flexion. The new method employed here shows, that motor SI can be measured during tonic movements. The findings also show motor SI during two real-life motor tasks: "writing" and "holding a pen". The lack of modulation of MEP amplitude during "holding bottle" task seems to indicate that SI is action specific rather than muscle specific.

Visual neuroscience: A brain area tuned for processing social interactions

Socialising with others is part of everyday life. A new study demonstrates that a brain area specialised for visual body perception is attuned to processing social interactions between two people. Intriguingly, this area is lateralised in the left hemisphere.

Reduced dimension stimulus decoding and column-based modeling reveal architectural differences of primary somatosensory finger maps between younger and older adults

The primary somatosensory cortex (SI) contains fine-grained tactile representations of the body, arranged in an orderly fashion. The use of ultra-high resolution fMRI data to detect group differences, for example between younger and older adults' SI maps, is challenging, because group alignment often does not preserve the high spatial detail of the data. Here, we use robust-shared response modeling (rSRM) that allows group analyses by mapping individual stimulus-driven responses to a lower dimensional shared feature space, to detect age-related differences in tactile representations between younger and older adults using 7T-fMRI data. Using this method, we show that finger representations are more precise in Brodmann-Area (BA) 3b and BA1 compared to BA2 and motor areas, and that this hierarchical processing is preserved across age groups. By combining rSRM with column-based decoding (C-SRM), we further show that the number of columns that optimally describes finger maps in SI is higher in younger compared to older adults in BA1, indicating a greater columnar size in older adults' SI. Taken together, we conclude that rSRM is suitable for finding fine-grained group differences in ultra-high resolution fMRI data, and we provide first evidence that the columnar architecture in SI changes with increasing age.

A custom MR-compatible dataglove for fMRI of the human motor cortex at 7T

We present a custom-built MR-compatible data glove to capture hand motion during concurrent fMRI experiments at 7 Tesla. Thermal and phantom tests showed our data glove can be used safely and without degradation of image quality. Subject-specific Blood Oxygen Level Dependent (BOLD) signal models, for use in fMRI analysis, were constructed based on recorded kinematic measurements. Experiments revealed the relative fMRI BOLD signal contribution of flexing, extending, and sustained isotonic extension. The ability to evaluate subject performance in real-time and create subject-specific BOLD signal models enables a wide range of experimental paradigms with improved data quality.Clinical Relevance- Using an MR compatible dataglove, subject specific Blood Oxygen Signal Level Dependent (BOLD) signal models can be constructed to study how the brain implements fine motor control.

A CMOS-based highly scalable flexible neural electrode interface

Perception, thoughts, and actions are encoded by the coordinated activity of large neuronal populations spread over large areas. However, existing electrophysiological devices are limited by their scalability in capturing this cortex-wide activity. Here, we developed an electrode connector based on an ultra-conformable thin-film electrode array that self-assembles onto silicon microelectrode arrays enabling multithousand channel counts at a millimeter scale. The interconnects are formed using microfabricated electrode pads suspended by thin support arms, termed Flex2Chip. Capillary-assisted assembly drives the pads to deform toward the chip surface, and van der Waals forces maintain this deformation, establishing Ohmic contact. Flex2Chip arrays successfully measured extracellular action potentials ex vivo and resolved micrometer scale seizure propagation trajectories in epileptic mice. We find that seizure dynamics in absence epilepsy in the Scn8a+/- model do not have constant propagation trajectories.

The 3D Structural Architecture of the Human Hand Area Is Nontopographic

The functional topography of the human primary somatosensory cortex hand area is a widely studied model system to understand sensory organization and plasticity. It is so far unclear whether the underlying 3D structural architecture also shows a topographic organization. We used 7 Tesla (7T) magnetic resonance imaging (MRI) data to quantify layer-specific myelin, iron, and mineralization in relation to population receptive field maps of individual finger representations in Brodman area 3b (BA 3b) of human S1 in female and male younger adults. This 3D description allowed us to identify a characteristic profile of layer-specific myelin and iron deposition in the BA 3b hand area, but revealed an absence of structural differences, an absence of low-myelin borders, and high similarity of 3D microstructure profiles between individual fingers. However, structural differences and borders were detected between the hand and face areas. We conclude that the 3D structural architecture of the human hand area is nontopographic, unlike in some monkey species, which suggests a high degree of flexibility for functional finger organization and a new perspective on human topographic plasticity.SIGNIFICANCE STATEMENT Using ultra-high-field MRI, we provide the first comprehensive in vivo description of the 3D structural architecture of the human BA 3b hand area in relation to functional population receptive field maps. High similarity of precise finger-specific 3D profiles, together with an absence of structural differences and an absence of low-myelin borders between individual fingers, reveals the 3D structural architecture of the human hand area to be nontopographic. This suggests reduced structural limitations to cortical plasticity and reorganization and allows for shared representational features across fingers.

A somato-cognitive action network alternates with effector regions in motor cortex

Motor cortex (M1) has been thought to form a continuous somatotopic homunculus extending down the precentral gyrus from foot to face representations1,2, despite evidence for concentric functional zones3 and maps of complex actions4. Here, using precision functional magnetic resonance imaging (fMRI) methods, we find that the classic homunculus is interrupted by regions with distinct connectivity, structure and function, alternating with effector-specific (foot, hand and mouth) areas. These inter-effector regions exhibit decreased cortical thickness and strong functional connectivity to each other, as well as to the cingulo-opercular network (CON), critical for action5 and physiological control6, arousal7, errors8 and pain9. This interdigitation of action control-linked and motor effector regions was verified in the three largest fMRI datasets. Macaque and pediatric (newborn, infant and child) precision fMRI suggested cross-species homologues and developmental precursors of the inter-effector system. A battery of motor and action fMRI tasks documented concentric effector somatotopies, separated by the CON-linked inter-effector regions. The inter-effectors lacked movement specificity and co-activated during action planning (coordination of hands and feet) and axial body movement (such as of the abdomen or eyebrows). These results, together with previous studies demonstrating stimulation-evoked complex actions4 and connectivity to internal organs10 such as the adrenal medulla, suggest that M1 is punctuated by a system for whole-body action planning, the somato-cognitive action network (SCAN). In M1, two parallel systems intertwine, forming an integrate-isolate pattern: effector-specific regions (foot, hand and mouth) for isolating fine motor control and the SCAN for integrating goals, physiology and body movement.

Evolution of cortical geometry and its link to function, behaviour and ecology

Studies in comparative neuroanatomy and of the fossil record demonstrate the influence of socio-ecological niches on the morphology of the cerebral cortex, but have led to oftentimes conflicting theories about its evolution. Here, we study the relationship between the shape of the cerebral cortex and the topography of its function. We establish a joint geometric representation of the cerebral cortices of ninety species of extant Euarchontoglires, including commonly used experimental model organisms. We show that variability in surface geometry relates to species' ecology and behaviour, independent of overall brain size. Notably, ancestral shape reconstruction of the cortical surface and its change during evolution enables us to trace the evolutionary history of localised cortical expansions, modal segregation of brain function, and their association to behaviour and cognition. We find that individual cortical regions follow different sequences of area increase during evolutionary adaptations to dynamic socio-ecological niches. Anatomical correlates of this sequence of events are still observable in extant species, and relate to their current behaviour and ecology. We decompose the deep evolutionary history of the shape of the human cortical surface into spatially and temporally conscribed components with highly interpretable functional associations, highlighting the importance of considering the evolutionary history of cortical regions when studying their anatomy and function.

Development of top-down cortical propagations in youth

Hierarchical processing requires activity propagating between higher- and lower-order cortical areas. However, functional neuroimaging studies have chiefly quantified fluctuations within regions over time rather than propagations occurring over space. Here, we leverage advances in neuroimaging and computer vision to track cortical activity propagations in a large sample of youth (n = 388). We delineate cortical propagations that systematically ascend and descend a cortical hierarchy in all individuals in our developmental cohort, as well as in an independent dataset of densely sampled adults. Further, we demonstrate that top-down, descending hierarchical propagations become more prevalent with greater demands for cognitive control as well as with development in youth. These findings emphasize that hierarchical processing is reflected in the directionality of propagating cortical activity and suggest top-down propagations as a potential mechanism of neurocognitive maturation in youth.

Evaluating increases in sensitivity from NORDIC for diverse fMRI acquisition strategies

As the neuroimaging field moves towards detecting smaller effects at higher spatial resolutions, and faster sampling rates, there is increased attention given to the deleterious contribution of unstructured, thermal noise. Here, we critically evaluate the performance of a recently developed reconstruction method, termed NORDIC, for suppressing thermal noise using datasets acquired with various field strengths, voxel sizes, sampling rates, and task designs. Following minimal preprocessing, statistical activation (t-values) of NORDIC processed data was compared to the results obtained with alternative denoising methods. Additionally, we examined the consistency of the estimates of task responses at the single-voxel, single run level, using a finite impulse response (FIR) model. To examine the potential impact on effective image resolution, the overall smoothness of the data processed with different methods was estimated. Finally, to determine if NORDIC alters or removes temporal information important for modeling responses, we employed an exhaustive leave-p-out cross validation approach, using FIR task responses to predict held out timeseries, quantified using R². After NORDIC, the t-values are increased, an improvement comparable to what could be achieved by 1.5 voxels smoothing, and task events are clearly visible and have less cross-run error. These advantages are achieved with smoothness estimates increasing by less than 4%, while 1.5 voxel smoothing is associated with increases of over 140%. Cross-validated R²s based on the FIR models show that NORDIC is not measurably distorting the temporal structure of the data under this approach and is the best predictor of non-denoised time courses. The results demonstrate that analyzing 1 run of data after NORDIC produces results equivalent to using 2 to 3 original runs and that NORDIC performs equally well across a diverse array of functional imaging protocols. Significance Statement: For functional neuroimaging, the increasing availability of higher field strengths and ever higher spatiotemporal resolutions has led to concomitant increase in concerns about the deleterious effects of thermal noise. Historically this noise source was suppressed using methods that reduce spatial precision such as image blurring or averaging over a large number of trials or sessions, which necessitates large data collection efforts. Here, we critically evaluate the performance of a recently developed reconstruction method, termed NORDIC, which suppresses thermal noise. Across datasets varying in field strength, voxel sizes, sampling rates, and task designs, NORDIC produces substantial gains in data quality. Both conventional t-statistics derived from general linear models and coefficients of determination for predicting unseen data are improved. These gains match or even exceed those associated with 1 voxel Full Width Half Max image smoothing, however, even such small amounts of smoothing are associated with a 52% reduction in estimates of spatial precision, whereas the measurable difference in spatial precision is less than 4% following NORDIC.

Mechanisms of speed encoding in the human middle temporal cortex measured by 7T fMRI

Perception of dynamic scenes in our environment results from the evaluation of visual features such as the fundamental spatial and temporal frequency components of a moving object. The ratio between these two components represents the object's speed of motion. The human middle temporal cortex hMT+ has a crucial biological role in the direct encoding of object speed. However, the link between hMT+ speed encoding and the spatiotemporal frequency components of a moving object is still under explored. Here, we recorded high resolution 7T blood oxygen level-dependent BOLD responses to different visual motion stimuli as a function of their fundamental spatial and temporal frequency components. We fitted each hMT+ BOLD response with a 2D Gaussian model allowing for two different speed encoding mechanisms: (1) distinct and independent selectivity for the spatial and temporal frequencies of the visual motion stimuli; (2) pure tuning for the speed of motion. We show that both mechanisms occur but in different neuronal groups within hMT+, with the largest subregion of the complex showing separable tuning for the spatial and temporal frequency of the visual stimuli. Both mechanisms were highly reproducible within participants, reconciling single cell recordings from MT in animals that have showed both encoding mechanisms. Our findings confirm that a more complex process is involved in the perception of speed than initially thought and suggest that hMT+ plays a primary role in the evaluation of the spatial features of the moving visual input.

Targeted Optical Neural Stimulation: A New Era for Personalized Medicine

Targeted optical neural stimulation comprises infrared neural stimulation and optogenetics, which affect the nervous system through induced thermal transients and activation of light-sensitive proteins, respectively. The main advantage of this pair of optical tools is high functional selectivity, which conventional electrical stimulation lacks. Over the past 15 years, the mechanism, safety, and feasibility of optical stimulation techniques have undergone continuous investigation and development. When combined with other methods like optical imaging and high-field functional magnetic resonance imaging, the translation of optical stimulation to clinical practice adds high value. We review the theoretical foundations and current state of optical stimulation, with a particular focus on infrared neural stimulation as a potential bridge linking optical stimulation to personalized medicine.

Generalization indicates asymmetric and interactive control networks for multi-finger dexterous movements

Finger dexterity is manifested by coordinated patterns of muscle activity and generalization of learning across contexts. Some fingers flex, others extend, and some are immobile. Whether or not the neural control processes of these direction-specific actions are independent remains unclear. We characterized behavioral principles underlying learning and generalization of dexterous flexion and extension movements, within and across hands, using an isometric dexterity task that precisely measured finger individuation, force accuracy, and temporal synchronization. Two cohorts of participants trained for 3 days in either the flexion or extension direction. All dexterity measures in both groups showed post-training improvement, although finger extension exhibited inferior dexterity. Surprisingly, learning of finger extension generalized to the untrained flexion direction, but not vice versa. This flexion bias was also evident in the untrained hand. Our study indicates direction-specific control circuits for learning of finger flexion and extension that interact by partially, but asymmetrically, transferring between directions.

Morphology, Connectivity, and Encoding Features of Tactile and Motor Representations of the Fingers in the Human Precentral and Postcentral Gyrus

Despite the tight coupling between sensory and motor processing for fine manipulation in humans, it is not yet totally clear which specific properties of the fingers are mapped in the precentral and postcentral gyrus. We used fMRI to compare the morphology, connectivity, and encoding of the motor and tactile finger representations (FRs) in the precentral and postcentral gyrus of 25 5-fingered participants (8 females). Multivoxel pattern and structural and functional connectivity analyses demonstrated the existence of distinct motor and tactile FRs within both the precentral and postcentral gyrus, integrating finger-specific motor and tactile information. Using representational similarity analysis, we found that the motor and tactile FRs in the sensorimotor cortex were described by the perceived structure of the hand better than by the actual hand anatomy or other functional models (finger kinematics, muscles synergies). We then studied a polydactyly individual (i.e., with a congenital 6-fingered hand) showing superior manipulation abilities and divergent anatomic-functional hand properties. The perceived hand model was still the best model for tactile representations in the precentral and postcentral gyrus, while finger kinematics better described motor representations in the precentral gyrus. We suggest that, under normal conditions (i.e., in subjects with a standard hand anatomy), the sensorimotor representations of the 5 fingers in humans converge toward a model of perceived hand anatomy, deviating from the real hand structure, as the best synthesis between functional and structural features of the hand.SIGNIFICANCE STATEMENT Distinct motor and tactile finger representations exist in both the precentral and postcentral gyrus, supported by a finger-specific pattern of anatomic and functional connectivity across modalities. At the representational level, finger representations reflect the perceived structure of the hand, which might result from an adapting process harmonizing (i.e., uniformizing) the encoding of hand function and structure in the precentral and postcentral gyrus. The same analyses performed in an extremely rare polydactyly subject showed that the emergence of such representational geometry is also found in neuromechanical variants with different hand anatomy and function. However, the harmonization process across the precentral and postcentral gyrus might not be possible because of divergent functional-structural properties of the hand and associated superior manipulation abilities.

Brain activity associated with taste stimulation: A mechanism for neuroplastic change?

PURPOSE

Neuroplasticity may be enhanced by increasing brain activation and bloodflow in neural regions relevant to the target behavior. We administered precisely formulated and dosed taste stimuli to determine whether the associated brain activity patterns included areas that underlie swallowing control.

METHODS

Five taste stimuli (unflavored, sour, sweet-sour, lemon, and orange suspensions) were administered in timing-regulated and temperature-controlled 3 mL doses via a customized pump/tubing system to 21 healthy adults during functional magnetic resonance imaging (fMRI). Whole-brain analyses of fMRI data assessed main effects of taste stimulation as well as differential effects of taste profile.

RESULTS

Differences in brain activity associated with taste stimulation overall as well as specific stimulus type were observed in key taste and swallowing regions including the orbitofrontal cortex, insula, cingulate, and pre- and postcentral gyri. Overall, taste stimulation elicited increased activation in swallowing-related brain regions compared to unflavored trials. Different patterns of blood oxygen level-dependent (BOLD) signal were noted by taste profile. For most areas, sweet-sour and sour trials elicited increases in BOLD compared to unflavored trials within that region, whereas lemon and orange trials yielded reductions in BOLD. This was despite identical concentrations of citric acid and sweetener in the lemon, orange, and sweet-sour solutions.

CONCLUSIONS

These results suggest that neural activity in swallowing-relevant regions can be amplified with taste stimuli and may be differentially affected by specific properties within very similar taste profiles. These findings provide critical foundational information for interpreting disparities in previous studies of taste effects on brain activity and swallowing function, defining optimal stimuli to increase brain activity in swallowing-relevant regions, and harnessing taste to enhance neuroplasticity and recovery for persons with swallowing disorders.

Laminar VASO fMRI in focal hand dystonia patients

Focal Hand Dystonia (FHD) is a disabling movement disorder characterized by involuntary movements, cramps and spasms. It is associated with pathological neural microcircuits in the cortical somatosensory system. While invasive preclinical modalities allow researchers to probe specific neural microcircuits of cortical layers and columns, conventional functional magnetic resonance imaging (fMRI) cannot resolve such small neural computational units. In this study, we take advantage of recent developments in ultra-high-field MRI hardware and MR-sequences to capture altered digit representations and laminar processing in FHD patients. We aim to characterize the capability and challenges of layer-specific imaging and analysis tools in resolving laminar and columnar structures in clinical research setups. We scanned N = 4 affected and N = 5 unaffected hemispheres at 7T and found consistent results of altered neural microcircuitry in FHD patients: 1) In affected hemispheres of FHD patients, we found a breakdown of ordered finger representation in the primary somatosensory cortex, as suggested from previous low-resolution fMRI. 2) In affected primary motor cortices of FHD patients, we furthermore found increased fMRI activity in superficial cortico-cortical neural input layers (II/III), compared to relatively weaker activity in the cortico-spinal output layers (Vb/VI). Overall, we show that layer-fMRI acquisition and analysis tools have the potential to address clinically-driven neuroscience research questions about altered computational mechanisms at the spatial scales that were previously only accessible in animal models. We believe that this study paves the way for easier translation of preclinical work into clinical research in focal hand dystonia and beyond.

Towards functional spin-echo BOLD line-scanning in humans at 7T

OBJECTIVE

Neurons cluster into sub-millimeter spatial structures and neural activity occurs at millisecond resolutions; hence, ultimately, high spatial and high temporal resolutions are required for functional MRI. In this work, we implemented a spin-echo line-scanning (SELINE) sequence to use in high spatial and temporal resolution fMRI.

MATERIALS AND METHODS

A line is formed by simply rotating the spin-echo refocusing gradient to a plane perpendicular to the excited slice and by removing the phase-encoding gradient. This technique promises a combination of high spatial and temporal resolution (250 μm, 500 ms) and microvascular specificity of functional responses. We compared SELINE data to a corresponding gradient-echo version (GELINE).

RESULTS

We demonstrate that SELINE showed much-improved line selection (i.e. a sharper line profile) compared to GELINE, albeit at the cost of a significant drop in functional sensitivity.

DISCUSSION

This low functional sensitivity needs to be addressed before SELINE can be applied for neuroscientific purposes.

Improved Selectivity in 7 T Digit Mapping Using VASO-CBV

Functional magnetic resonance imaging (fMRI) at Ultra-high field (UHF, ≥ 7 T) benefits from significant gains in the BOLD contrast-to-noise ratio (CNR) and temporal signal-to-noise ratio (tSNR) compared to conventional field strengths (3 T). Although these improvements enabled researchers to study the human brain to unprecedented spatial resolution, the blood pooling effect reduces the spatial specificity of the widely-used gradient-echo BOLD acquisitions. In this context, vascular space occupancy (VASO-CBV) imaging may be advantageous since it is proposed to have a higher spatial specificity than BOLD. We hypothesized that the assumed higher specificity of VASO-CBV imaging would translate to reduced overlap in fine-scale digit representation maps compared to BOLD-based digit maps. We used sub-millimeter resolution VASO fMRI at 7 T to map VASO-CBV and BOLD responses simultaneously in the motor and somatosensory cortices during individual finger movement tasks. We assessed the cortical overlap in different ways, first by calculating similarity coefficient metrics (DICE and Jaccard) and second by calculating selectivity measures. In addition, we demonstrate a consistent topographical organization of the targeted digit representations (thumb-index-little finger) in the motor areas. We show that the VASO-CBV responses yielded less overlap between the digit clusters than BOLD, and other selectivity measures were higher for VASO-CBV too. In summary, these results were consistent across metrics and participants, confirming the higher spatial specificity of VASO-CBV compared to BOLD.

Thalamocortical and corticothalamic connections of multiple functional domains in posterior parietal cortex of galagos

In prosimian galagos, the posterior parietal cortex (PPC) is subdivided into a number of functional domains where long-train intracortical microstimulation evoked different types of complex movements. Here, we placed anatomical tracers in multiple locations of PPC to reveal the origins and targets of thalamic connections of four PPC domains for different types of hindlimb, forelimb, or face movements. Thalamic connections of all four domains included nuclei of the motor thalamus, ventral anterior and ventral lateral nuclei, as well as parts of the sensory thalamus, the anterior pulvinar, posterior and ventral posterior superior nuclei, consistent with the sensorimotor functions of PPC domains. PPC domains also projected to the thalamic reticular nucleus in a somatotopic pattern. Quantitative differences in the distributions of labeled neurons in thalamic nuclei suggested that connectional patterns of these domains differed from each other.

Improved laminar specificity and sensitivity by combining SE and GE BOLD signals

The most widely used gradient-echo (GE) blood oxygenation level-dependent (BOLD) contrast has high sensitivity, but low specificity due to draining vein contributions, while spin-echo (SE) BOLD approach at ultra-high magnetic fields is highly specific to neural active sites but has lower sensitivity. To obtain high specificity and sensitivity, we propose to utilize a vessel-size-sensitive filter to the GE-BOLD signal, which suppresses macrovascular contributions and to combine selectively retained microvascular GE-BOLD signals with the SE-BOLD signals. To investigate our proposed idea, fMRI with 0.8 mm isotropic resolution was performed on the primary motor and sensory cortices in humans at 7 T by implementing spin- and gradient-echo (SAGE) echo planar imaging (EPI) acquisition. Microvascular-passed sigmoidal filters were designed based upon the vessel-size-sensitive ΔR₂*/ΔR₂ value for retaining GE-BOLD signals originating from venous vessels with ≤ 45 μm and ≤ 65 μm diameter. Unlike GE-BOLD fMRI, the laminar profile of SAGE-BOLD fMRI with the vessel-size-sensitive filter peaked at ∼ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity over GE-BOLD fMRI. Also, the functional sensitivity of SAGE BOLD at middle layers (0.75-1.5 mm) was improved by ∼ 80% to ∼100% when compared with SE BOLD. In summary, we showed that combined GE- and SE-BOLD fMRI with the vessel-size-sensitive filter indeed yielded improved laminar specificity and sensitivity and is therefore an excellent tool for high spatial resolution ultra-high filed (UHF)-fMRI studies for resolving mesoscopic functional units.

Magnetization transfer weighted laminar fMRI with multi-echo FLASH

Laminar functional magnetic resonance imaging (fMRI) using the gradient echo (GRE) blood oxygenation level dependent (BOLD) contrast is prone to signal changes arising from large unspecific venous vessels. Alternatives based on changes of cerebral blood volume (CBV) become more popular since it is expected that this hemodynamic response is dominant in microvasculature. One approach to sensitize the signal toward changes in CBV, and to simultaneously reduce unwanted extravascular (EV) BOLD blurring, is to selectively reduce gray matter (GM) signal via magnetization transfer (MT). In this work, we use off-resonant MT-pulses with a 3D FLASH readout to perform MT-prepared (MT-prep) laminar fMRI of the primary visual cortex (V1) at multiple echo times at 7 T. With a GRE-BOLD contrast without additional MT-weighting as reference, we investigated the influence of the MT-preparation on the shape and the echo time dependency of laminar profiles. Through numerical simulations, we optimized the sequence parameters to increase the sensitivity toward signal changes induced by changes in arterial CBV and to delineate the contributions of different compartments to the signal. We show that at 7 T, GM signals can be reduced by 30 %. Our laminar fMRI responses exhibit an increased signal change in the parenchyma at very short TE compared to a BOLD-only reference as a result of reduced EV signal intensity. By varying echo times, we could show that MT-prep results in less sensitivity toward unwanted signal changes based on changes in T₂*. We conclude that when accounting for nuclear overhauser enhancement effects in blood, off-resonant MT-prep combined with efficient short TE readouts can become a promising method to reduce unwanted EV venous contributions in GRE-BOLD and/or to allow scanning at much shorter echo times without incurring a sensitivity penalty in laminar fMRI.

Individual finger movement decoding using a novel ultra-high-density electroencephalography-based brain-computer interface system

Brain-Computer Interface (BCI) technology enables users to operate external devices without physical movement. Electroencephalography (EEG) based BCI systems are being actively studied due to their high temporal resolution, convenient usage, and portability. However, fewer studies have been conducted to investigate the impact of high spatial resolution of EEG on decoding precise body motions, such as finger movements, which are essential in activities of daily living. Low spatial sensor resolution, as found in common EEG systems, can be improved by omitting the conventional standard of EEG electrode distribution (the international 10-20 system) and ordinary mounting structures (e.g., flexible caps). In this study, we used newly proposed flexible electrode grids attached directly to the scalp, which provided ultra-high-density EEG (uHD EEG). We explored the performance of the novel system by decoding individual finger movements using a total of 256 channels distributed over the contralateral sensorimotor cortex. Dense distribution and small-sized electrodes result in an inter-electrode distance of 8.6 mm (uHD EEG), while that of conventional EEG is 60 to 65 mm on average. Five healthy subjects participated in the experiment, performed single finger extensions according to a visual cue, and received avatar feedback. This study exploits mu (8-12 Hz) and beta (13-25 Hz) band power features for classification and topography plots. 3D ERD/S activation plots for each frequency band were generated using the MNI-152 template head. A linear support vector machine (SVM) was used for pairwise finger classification. The topography plots showed regular and focal post-cue activation, especially in subjects with optimal signal quality. The average classification accuracy over subjects was 64.8 (6.3)%, with the middle versus ring finger resulting in the highest average accuracy of 70.6 (9.4)%. Further studies are required using the uHD EEG system with real-time feedback and motor imagery tasks to enhance classification performance and establish the basis for BCI finger movement control of external devices.

No robust online effects of transcranial direct current stimulation on corticospinal excitability

Transcranial direct current stimulation (tDCS) has been used for over twenty years to modulate cortical (particularly motor corticospinal) excitability both during (online) and outlasting (offline) the stimulation, with the former effects associated to the latter. However, tDCS effects are highly variable, partially because stimulation intensity is commonly not adjusted individually (in contrast to transcranial magnetic stimulation, TMS). In Experiment 1, we therefore explored an empirical approach of personalizing tDCS intensity for the primary motor cortex (M1) based on dose-response curves (DRCs), individually relating tDCS Intensity (in steps from 0.3 to 2.0 mA) and Polarity (anodal, cathodal) to the online modulation of concurrent TMS motor evoked potentials (MEP), assessing DRC reliability across two separate days. No robust DRCs could be observed, neither at the individual nor at the group level, with the only robust effect being a (paradoxical) MEP facilitation during cathodal tDCS at 2.0 mA, but no modulation at traditional intensities of or near 1 mA. In Experiment 2, we therefore attempted to replicate the classical bidirectional online MEP modulation during 1 mA tDCS that had been reported by several of the early seminal tDCS papers. We either closely recreated stimulation parameters and temporal protocol of these original studies (Experiment 2A) or slightly modernized them according to current standards (Experiment 2B). In neither experiment did we observed any significant online MEP modulation. We conclude that an empirical titration of individually effective tDCS intensities may not be feasible as online tDCS effects do not appear to be sufficiently robust.

An fMRI dataset for whole-body somatotopic mapping in humans

The somatotopic representation of the body is a well-established organizational principle in the human brain. Classic invasive direct electrical stimulation for somatotopic mapping cannot be used to map the whole-body topographical representation of healthy individuals. Functional magnetic resonance imaging (fMRI) has become an indispensable tool for the noninvasive investigation of somatotopic organization of the human brain using voluntary movement tasks. Unfortunately, body movements during fMRI scanning often cause large head motion artifacts. Consequently, there remains a lack of publicly accessible fMRI datasets for whole-body somatotopic mapping. Here, we present public high-resolution fMRI data to map the somatotopic organization based on motor movements in a large cohort of healthy adults (N = 62). In contrast to previous studies that were mostly designed to distinguish few body representations, most body parts are considered, including toe, ankle, leg, finger, wrist, forearm, upper arm, jaw, lip, tongue, and eyes. Moreover, the fMRI data are denoised by combining spatial independent component analysis with manual identification to clean artifacts from head motion associated with body movements.

Tactile perception of right middle fingertip suppresses excitability of motor cortex supplying right first dorsal interosseous muscle

The present study examined whether tactile perception of the fingertip modulates excitability of the motor cortex supplying the intrinsic hand muscle and whether this modulation is specific to the fingertip stimulated and the muscle and hand tested. Tactile stimulation was given to one of the five fingertips in the left or right hand, and transcranial magnetic stimulation eliciting motor evoked potential in the first dorsal interosseous muscle (FDI) or abductor digiti minimi was given 200 ms after the onset of tactile stimulation. The corticospinal excitability of the FDI at rest was suppressed by the tactile stimulation of the right middle fingertip, but such suppression was absent for the other fingers stimulated and for the other muscle or hand tested. The persistence and amplitude of the F-wave was not significantly influenced by tactile stimulation of the fingertip in the right hand. These findings indicate that tactile perception of the right middle fingertip suppresses excitability of the motor cortex supplying the right FDI at rest. The suppression of corticospinal excitability was absent during tonic contraction of the right FDI, indicating that the motor execution process interrupts the tactile perception-induced suppression of motor cortical excitability supplying the right FDI. These findings are in line with a view that the tactile perception of the right middle finger induces surround inhibition of the motor cortex supplying the prime mover of the finger neighboring the stimulated finger.

Pre-processing of Sub-millimeter GE-BOLD fMRI Data for Laminar Applications

Over the past 30 years, brain function has primarily been evaluated non-invasively using functional magnetic resonance imaging (fMRI) with gradient-echo (GE) sequences to measure blood-oxygen-level-dependent (BOLD) signals. Despite the multiple advantages of GE sequences, e.g., higher signal-to-noise ratio, faster acquisitions, etc., their relatively inferior spatial localization compromises the routine use of GE-BOLD in laminar applications. Here, in an attempt to rescue the benefits of GE sequences, we evaluated the effect of existing pre-processing methods on the spatial localization of signals obtained with EPIK, a GE sequence that affords voxel volumes of 0.25 mm3 with near whole-brain coverage. The methods assessed here apply to both task and resting-state fMRI data assuming the availability of reconstructed magnitude and phase images.

Mapping of whole-cerebrum resting-state networks using ultra-high resolution acquisition protocols

Resting‐state functional magnetic resonance imaging (fMRI) has been used in numerous studies to map networks in the brain that employ spatially disparate regions. However, attempts to map networks with high spatial resolution have been hampered by conflicting technical demands and associated problems. Results from recent fMRI studies have shown that spatial resolution remains around 0.7 × 0.7 × 0.7 mm³, with only partial brain coverage. Therefore, this work aims to present a novel fMRI technique that was developed based on echo‐planar‐imaging with keyhole (EPIK) combined with repetition‐time‐external (TR‐external) EPI phase correction. Each technique has been previously shown to be effective in enhancing the spatial resolution of fMRI, and in this work, the combination of the two techniques into TR‐external EPIK provided a nominal spatial resolution of 0.51 × 0.51 × 1.00 mm³ (0.26 mm³ voxel) with whole‐cerebrum coverage. Here, the feasibility of using half‐millimetre in‐plane TR‐external EPIK for resting‐state fMRI was validated using 13 healthy subjects and the corresponding reproducible mapping of resting‐state networks was demonstrated. Furthermore, TR‐external EPIK enabled the identification of various resting‐state networks distributed throughout the brain from a single fMRI session, with mapping fidelity onto the grey matter at 7T. The high‐resolution functional image further revealed mesoscale anatomical structures, such as small cerebral vessels and the internal granular layer of the cortex within the postcentral gyrus.

Beyond body maps: Information content of specific body parts is distributed across the somatosensory homunculus

The homunculus in primary somatosensory cortex (S1) is famous for its body part selectivity, but this dominant feature may eclipse other representational features, e.g., information content, also relevant for S1 organization. Using multivariate fMRI analysis, we ask whether body part information content can be identified in S1 beyond its primary region. Throughout S1, we identify significant representational dissimilarities between body parts but also subparts in distant non-primary regions (e.g., between the hand and the lips in the foot region and between different face parts in the foot region). Two movements performed by one body part (e.g., the hand) could also be dissociated well beyond its primary region (e.g., in the foot and face regions), even within Brodmann area 3b. Our results demonstrate that information content is more distributed across S1 than selectivity maps suggest. This finding reveals underlying information contents in S1 that could be harnessed for rehabilitation and brain-machine interfaces.

Analysis of Intracerebral Activity during Reflex Locomotion Stimulation According to Vojta’s Principle

Vojta’s therapy is a widely used approach in both the prevention and therapy of musculoskeletal disorders. Changes in the musculoskeletal system have been described repeatedly, but the principles of the approach have not yet been clarified. The objective of our study was to evaluate changes of intracerebral activity using electromagnetic tomography (sLORETA) that arise during reflex locomotion stimulation of the breast trigger zone according to Vojta’s therapy. Seventeen healthy women took part in the experiment (aged 20–30 years old). EEG activity was recorded 5 min prior to the reflex locomotion stimulation, during stimulation, and 5 min after the stimulation. The obtained data were subsequently processed in the sLORETA program and statistically evaluated at the significance level p ≤ 0.05. The analysis found statistically significant differences in the frequency bands alpha-2, beta-1, and beta-2 between the condition prior to stimulation and the actual stimulation in BAs 6, 7, 23, 24, and 31 and between the resting condition prior to stimulation, and the condition after the stimulation was terminated in the frequency bands alpha-1, alpha-2, beta-1, and beta-2 in BAs 3, 4, 6, and 24. The results showed that reflex locomotion stimulation according to Vojta’s therapy modulates electrical activity in the brain areas responsible for movement planning and regulating and performing the movement.

Representational Similarity Scores of Digits in the Sensorimotor Cortex Are Associated with Behavioral Performance

Previous studies aimed to unravel a digit-specific somatotopy in the primary sensorimotor (SM1) cortex. However, it remains unknown whether digit somatotopy is associated with motor preparation and/or motor execution during different types of tasks. We adopted multivariate representational similarity analysis to explore digit activation patterns in response to a finger tapping task (FTT). Sixteen healthy young adults underwent magnetic resonance imaging, and additionally performed an out-of-scanner choice reaction time task (CRTT) to assess digit selection performance. During both the FTT and CRTT, force data of all digits were acquired using force transducers. This allowed us to assess execution-related interference (i.e., digit enslavement; obtained from FTT & CRTT), as well as planning-related interference (i.e., digit selection deficit; obtained from CRTT) and determine their correlation with digit representational similarity scores of SM1. Findings revealed that digit enslavement during FTT was associated with contralateral SM1 representational similarity scores. During the CRTT, digit enslavement of both hands was also associated with representational similarity scores of the contralateral SM1. In addition, right hand digit selection performance was associated with representational similarity scores of left S1. In conclusion, we demonstrate a cortical origin of digit enslavement, and uniquely reveal that digit selection is associated with digit representations in primary somatosensory cortex (S1). Significance statement In current systems neuroscience, it is of critical importance to understand the relationship between brain function and behavioral outcome. With the present work, we contribute significantly to this understanding by uniquely assessing how digit representations in the sensorimotor cortex are associated with planning- and execution-related digit interference during a continuous finger tapping and a choice reaction time task. We observe that digit enslavement (i.e., execution-related interference) finds its origin in contralateral digit representations of SM1, and that deficits in digit selection (i.e., planning-related interference) in the right hand during a choice reaction time task are associated with more overlapping digit representations in left S1. This knowledge sheds new light on the functional contribution of the sensorimotor cortex to everyday motor skills.

Response Systems, Antagonistic Responses, and the Behavioral Repertoire

While response systems are often mentioned in the behavioral and physiological literature, an explicit discussion of what response systems are is lacking. Here we argue that response systems can be understood as an interaction between anatomically constrained behavioral topographies occasioned by currently present stimuli and a history of reinforcement. “New” response systems can develop during the lifetime as the organism gains instrumental control of new fine-grained topographies. Within this framework, antagonistic responses compete within each response system based on environmental stimulation, and competition is resolved at the striatum-thalamo-cortical loops level. While response systems can be by definition independent from one another, separate systems are often recruited at the same time to engage in complex responses, which themselves may be selected by reinforcement as functional units.

Comparing BOLD and VASO-CBV population receptive field estimates in human visual cortex

Vascular Space Occupancy (VASO) is an alternative fMRI approach based on changes in Cerebral Blood Volume (CBV). VASO-CBV fMRI can provide higher spatial specificity than the blood oxygenation level-dependent (BOLD) method because the CBV response is thought to be limited to smaller vessels. To investigate how this technique compares to BOLD fMRI for cognitive neuroscience applications, we compared population receptive field (pRF) mapping estimates between BOLD and VASO-CBV. We hypothesized that VASO-CBV would elicit distinct pRF properties compared to BOLD. Specifically, since pRF size estimates also depend on vascular sources, we hypothesized that reduced vascular blurring might yield narrower pRFs for VASO-CBV measurements. We used a VASO sequence with a double readout 3D EPI sequence at 7T to simultaneously measure VASO-CBV and BOLD responses in the visual cortex while participants viewed conventional pRF mapping stimuli. Both VASO-CBV and BOLD images show similar eccentricity and polar angle maps across all participants. Compared to BOLD-based measurements, VASO-CBV yielded lower tSNR and variance explained. The pRF size changed with eccentricity similarly for VASO-CBV and BOLD, and the pRF size estimates were similar for VASO-CBV and BOLD, even when we equate variance explained between VASO-CBV and BOLD. This result suggests that the vascular component of the pRF size is not dominating in either VASO-CBV or BOLD.

Cortical depth-dependent modeling of visual hemodynamic responses

A physiologically based three-dimensional (3D) hemodynamic model is developed to predict the experimentally observed blood oxygen level dependent (BOLD) responses versus the cortical depth induced by visual stimuli. Prior 2D approximations are relaxed in order to analyze 3D blood flow dynamics as a function of cortical depth. Comparison of the predictions with experimental data for evoked stimuli demonstrates that the full 3D model performs at least as well as previous approaches while remaining parsimonious. In particular, the 3D model requires significantly fewer assumptions and model parameters than previous models such that there is no longer need to define depth-specific parameter values for spatial spreading, peak amplitude, and hemodynamic velocity.

A dual larynx motor networks hypothesis

Humans are vocal modulators par excellence. This ability is supported in part by the dual representation of the laryngeal muscles in the motor cortex. Movement, however, is not the product of motor cortex alone but of a broader motor network. This network consists of brain regions that contain somatotopic maps that parallel the organization in motor cortex. We therefore present a novel hypothesis that the dual laryngeal representation is repeated throughout the broader motor network. In support of the hypothesis, we review existing literature that demonstrates the existence of network-wide somatotopy and present initial evidence for the hypothesis' plausibility. Understanding how this uniquely human phenotype in motor cortex interacts with broader brain networks is an important step toward understanding how humans evolved the ability to speak. We further suggest that this system may provide a means to study how individual components of the nervous system evolved within the context of neuronal networks. This article is part of the theme issue 'Voice modulation: from origin and mechanism to social impact (Part I)'.

Efficient high-resolution TMS mapping of the human motor cortex by nonlinear regression

Transcranial magnetic stimulation (TMS) is a powerful tool to investigate causal structure-function relationships in the human brain. However, a precise delineation of the effectively stimulated neuronal populations is notoriously impeded by the widespread and complex distribution of the induced electric field. Here, we propose a method that allows rapid and feasible cortical localization at the individual subject level. The functional relationship between electric field and behavioral effect is quantified by combining experimental data with numerically modeled fields to identify the cortical origin of the modulated effect. Motor evoked potentials (MEPs) from three finger muscles were recorded for a set of random stimulations around the primary motor area. All induced electric fields were nonlinearly regressed against the elicited MEPs to identify their cortical origin. We could distinguish cortical muscle representation with high spatial resolution and localized them primarily on the crowns and rims of the precentral gyrus. A post-hoc analysis revealed exponential convergence of the method with the number of stimulations, yielding a minimum of about 180 random stimulations to obtain stable results. Establishing a functional link between the modulated effect and the underlying mode of action, the induced electric field, is a fundamental step to fully exploit the potential of TMS. In contrast to previous approaches, the presented protocol is particularly easy to implement, fast to apply, and very robust due to the random coil positioning and therefore is suitable for practical and clinical applications.

Four methods of brain pattern analyses of fMRI signals associated with wrist extension versus wrist flexion studied for potential use in future motor learning BCI

OBJECTIVE

In stroke survivors, a treatment-resistant problem is inability to volitionally differentiate upper limb wrist extension versus flexion. When one intends to extend the wrist, the opposite occurs, wrist flexion, rendering the limb non-functional. Conventional therapeutic approaches have had limited success in achieving functional recovery of patients with chronic and severe upper extremity impairments. Functional magnetic resonance imaging (fMRI) neurofeedback is an emerging strategy that has shown potential for stroke rehabilitation. There is a lack of information regarding unique blood-oxygenation-level dependent (BOLD) cortical activations uniquely controlling execution of wrist extension versus uniquely controlling wrist flexion. Therefore, a first step in providing accurate neural feedback and training to the stroke survivor is to determine the feasibility of classifying (or differentiating) brain activity uniquely associated with wrist extension from that of wrist flexion, first in healthy adults.

APPROACH

We studied brain signal of 10 healthy adults, who performed wrist extension and wrist flexion during fMRI data acquisition. We selected four types of analyses to study the feasibility of differentiating brain signal driving wrist extension versus wrist flexion, as follows: 1) general linear model (GLM) analysis; 2) support vector machine (SVM) classification; 3) 'Winner Take All'; and 4) Relative Dominance.

RESULTS

With these four methods and our data, we found that few voxels were uniquely active during either wrist extension or wrist flexion. SVM resulted in only minimal classification accuracies. There was no significant difference in activation magnitude between wrist extension versus flexion; however, clusters of voxels showed extension signal > flexion signal and other clusters vice versa. Spatial patterns of activation differed among subjects.

SIGNIFICANCE

We encountered a number of obstacles to obtaining clear group results in healthy adults. These obstacles included the following: high variability across healthy adults in all measures studied; close proximity of uniquely active voxels to voxels that were common to both the extension and flexion movements; in general, higher magnitude of signal for the voxels common to both movements versus the magnitude of any given uniquely active voxel for one type of movement. Our results indicate that greater precision in imaging will be required to develop a truly effective method for differentiating wrist extension versus wrist flexion from fMRI data.

LayNii: A software suite for layer-fMRI

High-resolution fMRI in the sub-millimeter regime allows researchers to resolve brain activity across cortical layers and columns non-invasively. While these high-resolution data make it possible to address novel questions of directional information flow within and across brain circuits, the corresponding data analyses are challenged by MRI artifacts, including image blurring, image distortions, low SNR, and restricted coverage. These challenges often result in insufficient spatial accuracy of conventional analysis pipelines. Here we introduce a new software suite that is specifically designed for layer-specific functional MRI: LayNii. This toolbox is a collection of command-line executable programs written in C/C++ and is distributed opensource and as pre-compiled binaries for Linux, Windows, and macOS. LayNii is designed for layer-fMRI data that suffer from SNR and coverage constraints and thus cannot be straightforwardly analyzed in alternative software packages. Some of the most popular programs of LayNii contain 'layerification' and columnarization in the native voxel space of functional data as well as many other layer-fMRI specific analysis tasks: layer-specific smoothing, model-based vein mitigation of GE-BOLD data, quality assessment of artifact dominated sub-millimeter fMRI, as well as analyses of VASO data.

Evoking highly focal percepts in the fingertips through targeted stimulation of sulcal regions of the brain for sensory restoration

Background:

Paralysis and neuropathy, affecting millions of people worldwide, can be accompanied by significant loss of somatosensation. With tactile sensation being central to achieving dexterous movement, brain-computer interface (BCI) researchers have used intracortical and cortical surface electrical stimulation to restore somatotopically-relevant sensation to the hand. However, these approaches are restricted to stimulating the gyral areas of the brain. Since representation of distal regions of the hand extends into the sulcal regions of human primary somatosensory cortex (S1), it has been challenging to evoke sensory percepts localized to the fingertips.

Objective/hypothesis:

Targeted stimulation of sulcal regions of S1, using stereoelectroencephalography (SEEG) depth electrodes, can evoke focal sensory percepts in the fingertips.

Methods:

Two participants with intractable epilepsy received cortical stimulation both at the gyri via high-density electrocorticography (HD-ECoG) grids and in the sulci via SEEG depth electrode leads. We characterized the evoked sensory percepts localized to the hand.

Results:

We show that highly focal percepts can be evoked in the fingertips of the hand through sulcal stimulation. fMRI, myelin content, and cortical thickness maps from the Human Connectome Project elucidated specific cortical areas and sub-regions within S1 that evoked these focal percepts. Within-participant comparisons showed that percepts evoked by sulcal stimulation via SEEG electrodes were significantly more focal (80% less area; p = 0.02) and localized to the fingertips more often, than by gyral stimulation via HD-ECoG electrodes. Finally, sulcal locations with consistent modulation of high-frequency neural activity during mechanical tactile stimulation of the fingertips showed the same somatotopic correspondence as cortical stimulation.

Conclusions:

Our findings indicate minimally invasive sulcal stimulation via SEEG electrodes could be a clinically viable approach to restoring sensation.

Improvement of sensitivity and specificity for laminar BOLD fMRI with double spin-echo EPI in humans at 7 T

Mapping mesoscopic cortical functional units such as columns or laminae is increasingly pursued by ultra-high field (UHF) functional magnetic resonance imaging (fMRI). The most popular approach for high-resolution fMRI is currently gradient-echo (GE) blood oxygenation level-dependent (BOLD) fMRI. However, its spatial accuracy is reduced due to its sensitivity to draining vessels, including pial veins, whereas spin-echo (SE) BOLD signal is expected to have higher spatial accuracy, albeit with lower sensitivity than the GE-BOLD signal. Here, we introduce a new double spin-echo (dSE) echo-planar imaging (EPI) method to improve the sensitivity of SE-BOLD contrast by averaging two spin-echoes using three radiofrequency pulses. Human fMRI experiments were performed with slices perpendicular to the central sulcus between motor and sensory cortices at 7 T during fist-clenching with touching. First, we evaluated the feasibility of single-shot dSE-EPI for BOLD fMRI with 1.5 mm isotropic resolution and found that dSE-BOLD fMRI has higher signal-to-noise ratio (SNR), temporal SNR (tSNR), and higher functional sensitivity than conventional SE-BOLD fMRI. Second, to investigate the laminar specificity of dSE-BOLD fMRI, we implemented a multi-shot approach to achieve 0.8-mm isotropic resolution with sliding-window reconstruction. Unlike GE-BOLD fMRI, the cortical profile of dSE-BOLD fMRI peaked at ~ 1.0 mm from the surface of the primary motor and sensory cortices, demonstrating an improvement of laminar specificity in humans over GE-BOLD fMRI. The proposed multi-shot dSE-EPI method is viable for high spatial resolution UHF-fMRI studies in the pursuit of resolving mesoscopic functional units.

Human Somatosensory Processing and Artificial Somatosensation

In the past few years, we have gained a better understanding of the information processing mechanism in the human brain, which has led to advances in artificial intelligence and humanoid robots. However, among the various sensory systems, studying the somatosensory system presents the greatest challenge. Here, we provide a comprehensive review of the human somatosensory system and its corresponding applications in artificial systems. Due to the uniqueness of the human hand in integrating receptor and actuator functions, we focused on the role of the somatosensory system in object recognition and action guidance. First, the low-threshold mechanoreceptors in the human skin and somatotopic organization principles along the ascending pathway, which are fundamental to artificial skin, were summarized. Second, we discuss high-level brain areas, which interacted with each other in the haptic object recognition. Based on this close-loop route, we used prosthetic upper limbs as an example to highlight the importance of somatosensory information. Finally, we present prospective research directions for human haptic perception, which could guide the development of artificial somatosensory systems.

Typical somatomotor physiology of the hand is preserved in a patient with an amputated arm: An ECoG case study

Electrophysiological signals in the human motor system may change in different ways after deafferentation, with some studies emphasizing reorganization while others propose retained physiology. Understanding whether motor electrophysiology is retained over longer periods of time can be invaluable for patients with paralysis (e.g. ALS or brainstem stroke) when signals from sensorimotor areas may be used for communication or control over neural prosthetic devices. In addition, a maintained electrophysiology can potentially benefit the treatment of phantom limb pains through prolonged use of these signals in a brain-machine interface (BCI). Here, we were presented with the unique opportunity to investigate the physiology of the sensorimotor cortex in a patient with an amputated arm using electrocorticographic (ECoG) measurements. While implanted with an ECoG grid for clinical evaluation of electrical stimulation for phantom limb pain, the patient performed attempted finger movements with the contralateral (lost) hand and executed finger movements with the ipsilateral (healthy) hand. The electrophysiology of the sensorimotor cortex contralateral to the amputated hand remained very similar to that of hand movement in healthy people, with a spatially focused increase of high-frequency band (65-175 Hz; HFB) power over the hand region and a distributed decrease in low-frequency band (15-28 Hz; LFB) power. The representation of the three different fingers (thumb, index and little) remained intact and HFB patterns could be decoded using support vector learning at single-trial classification accuracies of >90%, based on the first 1-3 s of the HFB response. These results indicate that hand representations are largely retained in the motor cortex. The intact physiological response of the amputated hand, the high distinguishability of the fingers and fast temporal peak are encouraging for neural prosthetic devices that target the sensorimotor cortex.

Validating layer-specific VASO across species

Cerebral blood volume (CBV) has been shown to be a robust and important physiological parameter for quantitative interpretation of functional (f)MRI, capable of delivering highly localized mapping of neural activity. Indeed, with recent advances in ultra-high-field (≥7T) MRI hardware and associated sequence libraries, it has become possible to capture non-invasive CBV weighted fMRI signals across cortical layers. One of the most widely used approaches to achieve this (in humans) is through vascular-space-occupancy (VASO) fMRI. Unfortunately, the exact contrast mechanisms of layer-dependent VASO fMRI have not been validated for human fMRI and thus interpretation of such data is confounded. Here we validate the signal source of layer-dependent SS-SI VASO fMRI using multi-modal imaging in a rat model in response to neuronal activation (somatosensory cortex) and respiratory challenge (hypercapnia). In particular VASO derived CBV measures are directly compared to concurrent measures of total haemoglobin changes from high resolution intrinsic optical imaging spectroscopy (OIS). Quantified cortical layer profiling is demonstrated to be in agreement between VASO and contrast enhanced fMRI (using monocrystalline iron oxide nanoparticles, MION). Responses show high spatial localisation to layers of cortical processing independent of confounding large draining veins which can hamper BOLD fMRI studies, (depending on slice positioning). Thus, a cross species comparison is enabled using VASO as a common measure. We find increased VASO based CBV reactivity (3.1 ± 1.2 fold increase) in humans compared to rats. Together, our findings confirm that the VASO contrast is indeed a reliable estimate of layer-specific CBV changes. This validation study increases the neuronal interpretability of human layer-dependent VASO fMRI as an appropriate method in neuroscience application studies, in which the presence of large draining intracortical and pial veins limits neuroscientific inference with BOLD fMRI.