A touch of hierarchy: population receptive fields reveal fingertip integration in Brodmann areas in human primary somatosensory cortex

Several neuroimaging studies have shown the somatotopy of body part representations in primary somatosensory cortex (S1), but the functional hierarchy of distinct subregions in human S1 has not been adequately addressed. The current study investigates the functional hierarchy of cyto-architectonically distinct regions, Brodmann areas BA3, BA1, and BA2, in human S1. During functional MRI experiments, we presented participants with vibrotactile stimulation of the fingertips at three different vibration frequencies. Using population Receptive Field (pRF) modeling of the fMRI BOLD activity, we identified the hand region in S1 and the somatotopy of the fingertips. For each voxel, the pRF center indicates the finger that most effectively drives the BOLD signal, and the pRF size measures the spatial somatic pooling of fingertips. We find a systematic relationship of pRF sizes from lower-order areas to higher-order areas. Specifically, we found that pRF sizes are smallest in BA3, increase slightly towards BA1, and are largest in BA2, paralleling the increase in visual receptive field size as one ascends the visual hierarchy. Additionally, we find that the time-to-peak of the hemodynamic response in BA3 is roughly 0.5 s earlier compared to BA1 and BA2, further supporting the notion of a functional hierarchy of subregions in S1. These results were obtained during stimulation of different mechanoreceptors, suggesting that different afferent fibers leading up to S1 feed into the same cortical hierarchy

A Novel 2D Standard Cartesian Representation for the Human Sensorimotor Cortex

For some experimental approaches in brain imaging, the existing normalization techniques are not always sufficient. This may be the case if the anatomical shape of the region of interest varies substantially across subjects, or if one needs to compare the left and right hemisphere in the same subject. Here we propose a new standard representation, building upon existing normalization methods: Cgrid (Cartesian geometric representation with isometric dimensions). Cgrid is based on imposing a Cartesian grid over a cortical region of interest that is bounded by anatomical (atlas-based) landmarks. We applied this new representation to the sensorimotor cortex and we evaluated its performance by studying the similarity of activation patterns for hand, foot and tongue movements between subjects, and similarity between hemispheres within subjects. The Cgrid similarities were benchmarked against the similarities of activation patterns when transformed into standard MNI space using SPM, and to similarities from FreeSurfer's surface-based normalization. For both between-subject and between-hemisphere comparisons, similarity scores in Cgrid were high, similar to those from FreeSurfer normalization and higher than similarity scores from SPM's MNI normalization. This indicates that Cgrid allows for a straightforward way of representing and comparing sensorimotor activity patterns across subjects and between hemispheres of the same subjects

A fully synthetic three-dimensional human cerebrovascular model based on histological characteristics to investigate the hemodynamic fingerprint of the layer BOLD fMRI signal formation

Recent advances in functional magnetic resonance imaging (fMRI) at ultra-high field (≥7 tesla), novel hardware, and data analysis methods have enabled detailed research on neurovascular function, such as cortical layer-specific activity, in both human and nonhuman species. A widely used fMRI technique relies on the blood oxygen level-dependent (BOLD) signal. BOLD fMRI offers insights into brain function by measuring local changes in cerebral blood volume, cerebral blood flow, and oxygen metabolism induced by increased neuronal activity. Despite its potential, interpreting BOLD fMRI data is challenging as it is only an indirect measurement of neuronal activity. Computational modeling can help interpret BOLD data by simulating the BOLD signal formation. Current developments have focused on realistic 3D vascular models based on rodent data to understand the spatial and temporal BOLD characteristics. While such rodent-based vascular models highlight the impact of the angioarchitecture on the BOLD signal amplitude, anatomical differences between the rodent and human vasculature necessitate the development of human-specific models. Therefore, a computational framework integrating human cortical vasculature, hemodynamic changes, and biophysical properties is essential. Here, we present a novel computational approach: a three-dimensional VAscular MOdel based on Statistics (3D VAMOS), enabling the investigation of the hemodynamic fingerprint of the BOLD signal within a model encompassing a fully synthetic human 3D cortical vasculature and hemodynamics. Our algorithm generates microvascular and macrovascular architectures based on morphological and topological features from the literature on human cortical vasculature. By simulating specific oxygen saturation states and biophysical interactions, our framework characterizes the intravascular and extravascular signal contributions across cortical depth and voxel-wise levels for gradient-echo and spin-echo readouts. Thereby, the 3D VAMOS computational framework demonstrates that using human characteristics significantly affects the BOLD fingerprint, making it an essential step in understanding the fundamental underpinnings of layer-specific fMRI experiments

Translatability of preclinical to early clinical tolerable and pharmacologically active dose ranges for central nervous system active drugs

The primary purpose of this study was to assess the translatability of preclinical to early clinical tolerable and pharmacologically active dose ranges for central nervous system (CNS) active drugs. As a part of this, IBs were reviewed on reporting quality. Investigator's Brochures (IBs) of studies performed at the Centre for Human Drug Research (CHDR) reporting statistically significant results of CNS activity related to the drug's mechanism of action were included. The quality of IBs was assessed based on the presence of a rationale for the chosen animal model, completeness of pharmacokinetic (PK) results in reporting and internal validity information of the preclinical evidence. The IB-derisk tool was used to generate preclinical and early clinical data overviews data. For each compound, the overlap between pharmacologically active dose ranges and well-tolerated levels was calculated for three pharmacokinetic (PK) parameters: human equivalent dose (HED), maximum plasma concentration (C max) and area under the curve (AUC). Twenty-five IBs were included. In general, the quality of reporting in IBs was assessed as poor. About a third of studies did not explore the entire concentration-effect curve (pre)clinically. Single dose tolerability ranges were most accurately predicted by C max. Human equivalent dose and AUC were the best predictors of pharmacologically active ranges. Tolerable and pharmacologically active dose ranges in healthy volunteers can be reasonably well predicted from preclinical data with the IB-derisk tool. The translatability of preclinical studies can be improved by applying a higher reporting standard in IBs including comparable PK measurements across all preclinical and clinical studies

A global water resources ensemble of hydrological models: the eartH2Observe Tier-1 dataset

The dataset presented here consists of an ensemble of 10 global hydrological and land surface models for the period 1979–2012 using a reanalysis-based meteorological forcing dataset (0.5° resolution). The current dataset serves as a state of the art in current global hydrological modelling and as a benchmark for further improvements in the coming years. A signal-to-noise ratio analysis revealed low inter-model agreement over (i) snow-dominated regions and (ii) tropical rainforest and monsoon areas. The large uncertainty of precipitation in the tropics is not reflected in the ensemble runoff. Verification of the results against benchmark datasets for evapotranspiration, snow cover, snow water equivalent, soil moisture anomaly and total water storage anomaly using the tools from The International Land Model Benchmarking Project (ILAMB) showed overall useful model performance, while the ensemble mean generally outperformed the single model estimates. The results also show that there is currently no single best model for all variables and that model performance is spatially variable. In our unconstrained model runs the ensemble mean of total runoff into the ocean was 46 268 km3 yr−1 (334 kg m−2 yr−1), while the ensemble mean of total evaporation was 537 kg m−2 yr−1. All data are made available openly through a Water Cycle Integrator portal (WCI, wci.earth2observe.eu), and via a direct http and ftp download. The portal follows the protocols of the open geospatial consortium such as OPeNDAP, WCS and WMS. The DOI for the data is https://doi.org/10.1016/10.5281/zenodo.167070

Detailed somatotopy in primary motor and somatosensory cortex revealed by Gaussian population receptive fields

The relevance of human primary motor cortex (M1) for motor actions has long been established. However, it is still unknown how motor actions are represented, and whether M1 contains an ordered somatotopy at the mesoscopic level. In the current study we show that a detailed within-limb somatotopy can be obtained in M1 during finger movements using Gaussian population Receptive Field (pRF) models. Similar organizations were also obtained for primary somatosensory cortex (S1), showing that individual finger representations are interconnected throughout sensorimotor cortex. The current study additionally estimates receptive field sizes of neuronal populations, showing differences between finger digit representations, between M1 and S1, and additionally between finger digit flexion and extension. Using the Gaussian pRF approach, the detailed somatotopic organization of M1 can be obtained including underlying characteristics, allowing for the in-depth investigation of cortical motor representation and sensorimotor integration

Detailed somatotopy in primary motor and somatosensory cortex revealed by Gaussian population receptive fields

The relevance of human primary motor cortex (M1) for motor actions has long been established. However, it is still unknown how motor actions are represented, and whether M1 contains an ordered somatotopy at the mesoscopic level. In the current study we show that a detailed within-limb somatotopy can be obtained in M1 during finger movements using Gaussian population Receptive Field (pRF) models. Similar organizations were also obtained for primary somatosensory cortex (S1), showing that individual finger representations are interconnected throughout sensorimotor cortex. The current study additionally estimates receptive field sizes of neuronal populations, showing differences between finger digit representations, between M1 and S1, and additionally between finger digit flexion and extension. Using the Gaussian pRF approach, the detailed somatotopic organization of M1 can be obtained including underlying characteristics, allowing for the in-depth investigation of cortical motor representation and sensorimotor integration

Moving in on human motor cortex. Characterizing the relationship between body parts with non-rigid population response fields.

For cortical motor activity, the relationships between different body part representations is unknown. Through reciprocal body part relationships, functionality of cortical motor areas with respect to whole body motor control can be characterized. In the current study, we investigate the relationship between body part representations within individual neuronal populations in motor cortices, following a 7 Tesla fMRI 18-body-part motor experiment in combination with our newly developed non-rigid population Response Field (pRF) model and graph theory. The non-rigid pRF metrics reveal somatotopic structures in all included motor cortices covering frontal, parietal, medial and insular cortices and that neuronal populations in primary sensorimotor cortex respond to fewer body parts than secondary motor cortices. Reciprocal body part relationships are estimated in terms of uniqueness, clique-formation, and influence. We report unique response profiles for the knee, a clique of body parts surrounding the ring finger, and a central role for the shoulder and wrist. These results reveal associations among body parts from the perspective of the central nervous system, while being in agreement with intuitive notions of body part usage

Knowing left from right : asymmetric functional connectivity during resting state

The functional organization of left and right hemispheres is different, and hemispheric asymmetries are thought to underlie variations in brain function across individuals. In this study, we assess how differences between hemispheres are reflected in Asymmetric Functional Connectivity (AFC), which provides a full description of how the brain's connectivity structure during resting state differs from that of the same brain mirrored over the longitudinal fissure. In addition, we assess how AFC varies across subjects. Data were provided by the Human Connectome Project, including 423 resting state and combined language task fMRI data sets, and the pattern of AFC was established for all subjects. While we could quantify the symmetry of brain connectivity at 95%, significant asymmetries were observed, consisting foremost of: (1) higher correlations between language areas in the left hemisphere than between their right hemisphere homologues. (2) Higher correlations between language homologue areas in the right hemisphere and left default mode network, than between language areas in the left hemisphere and the default mode network in the right hemisphere. The extent to which subjects exhibited this pattern correlated with language lateralization and handedness. Further exploration in intersubject variation in AFC revealed several additional patterns, one involving entire hemispheres, and another correlations with limbic areas. These results show that language is an important, but not only determinant of AFC. The additional patterns of AFC require further research to be linked to specific asymmetric neuronal states or events

Knowing left from right : asymmetric functional connectivity during resting state

The functional organization of left and right hemispheres is different, and hemispheric asymmetries are thought to underlie variations in brain function across individuals. In this study, we assess how differences between hemispheres are reflected in Asymmetric Functional Connectivity (AFC), which provides a full description of how the brain's connectivity structure during resting state differs from that of the same brain mirrored over the longitudinal fissure. In addition, we assess how AFC varies across subjects. Data were provided by the Human Connectome Project, including 423 resting state and combined language task fMRI data sets, and the pattern of AFC was established for all subjects. While we could quantify the symmetry of brain connectivity at 95%, significant asymmetries were observed, consisting foremost of: (1) higher correlations between language areas in the left hemisphere than between their right hemisphere homologues. (2) Higher correlations between language homologue areas in the right hemisphere and left default mode network, than between language areas in the left hemisphere and the default mode network in the right hemisphere. The extent to which subjects exhibited this pattern correlated with language lateralization and handedness. Further exploration in intersubject variation in AFC revealed several additional patterns, one involving entire hemispheres, and another correlations with limbic areas. These results show that language is an important, but not only determinant of AFC. The additional patterns of AFC require further research to be linked to specific asymmetric neuronal states or events