Mapping the human subcortical auditory system using histology, postmortem MRI and in vivo MRI at 7T

Functional Localization of the Human Auditory and Visual Thalamus Using a Thalamic Localizer Functional Magnetic Resonance Imaging Task

Functional magnetic resonance imaging (fMRI) of the auditory and visual sensory systems of the human brain is an active area of investigation in the study of human health and disease. The medial geniculate nucleus (MGN) and lateral geniculate nucleus (LGN) are key thalamic nuclei involved in the processing and relay of auditory and visual information, respectively, and are the subject of blood-oxygen-level-dependent (BOLD) fMRI studies of neural activation and functional connectivity in human participants. However, localization of BOLD fMRI signal originating from neural activity in MGN and LGN remains a technical challenge, due in part to the poor definition of boundaries of these thalamic nuclei in standard T1-weighted and T2-weighted magnetic resonance imaging sequences. Here, we report the development and evaluation of an auditory and visual sensory thalamic localizer (TL) fMRI task that produces participant-specific functionally-defined regions of interest (fROIs) of both MGN and LGN, using 3 Tesla multiband fMRI and a clustered-sparse temporal acquisition sequence, in less than 16 minutes of scan time. We demonstrate the use of MGN and LGN fROIs obtained from the TL fMRI task in standard resting-state functional connectivity (RSFC) fMRI analyses in the same participants. In RSFC analyses, we validated the specificity of MGN and LGN fROIs for signals obtained from primary auditory and visual cortex, respectively, and benchmark their performance against alternative atlas- and segmentation-based localization methods. The TL fMRI task and analysis code (written in Presentation and MATLAB, respectively) have been made freely available to the wider research community.

Role of autonomic, nociceptive, and limbic brainstem nuclei in core autism features

Although multiple theories have speculated about the brainstem reticular formation's involvement in autistic behaviors, the in vivo imaging of brainstem nuclei needed to test these theories has proven technologically challenging. Using methods to improve brainstem imaging in children, this study set out to elucidate the role of the autonomic, nociceptive, and limbic brainstem nuclei in the autism features of 145 children (74 autistic children, 6.0-10.9 years). Participants completed an assessment of core autism features and diffusion- and T1-weighted imaging optimized to improve brainstem images. After data reduction via principal component analysis, correlational analyses examined associations among autism features and the microstructural properties of brainstem clusters. Independent replication was performed in 43 adolescents (24 autistic, 13.0-17.9 years). We found specific nuclei, most robustly the parvicellular reticular formation-alpha (PCRtA) and to a lesser degree the lateral parabrachial nucleus (LPB) and ventral tegmental parabrachial pigmented complex (VTA-PBP), to be associated with autism features. The PCRtA and some of the LPB associations were independently found in the replication sample, but the VTA-PBP associations were not. Consistent with theoretical perspectives, the findings suggest that individual differences in pontine reticular formation nuclei contribute to the prominence of autistic features. Specifically, the PCRtA, a nucleus involved in mastication, digestion, and cardio-respiration in animal models, was associated with social communication in children, while the LPB, a pain-network nucleus, was associated with repetitive behaviors. These findings highlight the contributions of key autonomic brainstem nuclei to the expression of core autism features.

Multiple Concurrent Predictions Inform Prediction Error in the Human Auditory Pathway

The key assumption of the predictive coding framework is that internal representations are used to generate predictions on how the sensory input will look like in the immediate future. These predictions are tested against the actual input by the so-called prediction error units, which encode the residuals of the predictions. What happens to prediction errors, however, if predictions drawn by different stages of the sensory hierarchy contradict each other? To answer this question, we conducted two fMRI experiments while female and male human participants listened to sequences of sounds: pure tones in the first experiment and frequency-modulated sweeps in the second experiment. In both experiments, we used repetition to induce predictions based on stimulus statistics (stats-informed predictions) and abstract rules disclosed in the task instructions to induce an orthogonal set of (task-informed) predictions. We tested three alternative scenarios: neural responses in the auditory sensory pathway encode prediction error with respect to (1) the stats-informed predictions, (2) the task-informed predictions, or (3) a combination of both. Results showed that neural populations in all recorded regions (bilateral inferior colliculus, medial geniculate body, and primary and secondary auditory cortices) encode prediction error with respect to a combination of the two orthogonal sets of predictions. The findings suggest that predictive coding exploits the non-linear architecture of the auditory pathway for the transmission of predictions. Such non-linear transmission of predictions might be crucial for the predictive coding of complex auditory signals like speech.Significance Statement Sensory systems exploit our subjective expectations to make sense of an overwhelming influx of sensory signals. It is still unclear how expectations at each stage of the processing pipeline are used to predict the representations at the other stages. The current view is that this transmission is hierarchical and linear. Here we measured fMRI responses in auditory cortex, sensory thalamus, and midbrain while we induced two sets of mutually inconsistent expectations on the sensory input, each putatively encoded at a different stage. We show that responses at all stages are concurrently shaped by both sets of expectations. The results challenge the hypothesis that expectations are transmitted linearly and provide for a normative explanation of the non-linear physiology of the corticofugal sensory system.

Segmentation of supragranular and infragranular layers in ultra-high resolution 7T ex vivo MRI of the human cerebral cortex

Accurate labeling of specific layers in the human cerebral cortex is crucial for advancing our understanding of neurodevelopmental and neurodegenerative disorders. Leveraging recent advancements in ultra-high resolution ex vivo MRI, we present a novel semi-supervised segmentation model capable of identifying supragranular and infragranular layers in ex vivo MRI with unprecedented precision. On a dataset consisting of 17 whole-hemisphere ex vivo scans at 120 μm, we propose a multi-resolution U-Nets framework (MUS) that integrates global and local structural information, achieving reliable segmentation maps of the entire hemisphere, with Dice scores over 0.8 for supra- and infragranular layers. This enables surface modeling, atlas construction, anomaly detection in disease states, and cross-modality validation, while also paving the way for finer layer segmentation. Our approach offers a powerful tool for comprehensive neuroanatomical investigations and holds promise for advancing our mechanistic understanding of progression of neurodegenerative diseases.

Impact of inner ear malformation and cochlear nerve deficiency on the development of auditory-language network in children with profound sensorineural hearing loss

Profound congenital sensorineural hearing loss (SNHL) prevents children from developing spoken language. Cochlear implantation and auditory brainstem implantation can provide partial hearing sensation, but language development outcomes can vary, particularly for patients with inner ear malformations and/or cochlear nerve deficiency (IEM&CND). Currently, the peripheral auditory structure is evaluated through visual inspection of clinical imaging, but this method is insufficient for surgical planning and prognosis. The central auditory pathway is also challenging to examine in vivo due to its delicate subcortical structures. Previous attempts to locate subcortical auditory nuclei using fMRI responses to sounds are not applicable to patients with profound hearing loss as no auditory brainstem responses can be detected in these individuals, making it impossible to capture corresponding blood oxygen signals in fMRI. In this study, we developed a new pipeline for mapping the auditory pathway using structural and diffusional MRI. We used a fixel-based approach to investigate the structural development of the auditory-language network for profound SNHL children with normal peripheral structure and those with IEM&CND under 6 years old. Our findings indicate that the language pathway is more sensitive to peripheral auditory condition than the central auditory pathway, highlighting the importance of early intervention for profound SNHL children to provide timely speech inputs. We also propose a comprehensive pre-surgical evaluation extending from the cochlea to the auditory-language network, showing significant correlations between age, gender, Cn.VIII median contrast value, and the language network with post-implant qualitative outcomes.

Mixed methodology in human brain research: integrating MRI and histology

Postmortem magnetic resonance imaging (MRI) can provide a bridge between histological observations and the in vivo anatomy of the human brain. Approaches aimed at the co-registration of data derived from the two techniques are gaining interest. Optimal integration of the two research fields requires detailed knowledge of the tissue property requirements for individual research techniques, as well as a detailed understanding of the consequences of tissue fixation steps on the imaging quality outcomes for both MRI and histology. Here, we provide an overview of existing studies that bridge between state-of-the-art imaging modalities, and discuss the background knowledge incorporated into the design, execution and interpretation of postmortem studies. A subset of the discussed challenges transfer to animal studies as well. This insight can contribute to furthering our understanding of the normal and diseased human brain, and to facilitate discussions between researchers from the individual disciplines.

The effects of exposure to road traffic noise at school on central auditory pathway functional connectivity

As the world becomes more urbanized, more people become exposed to traffic and the risks associated with a higher exposure to road traffic noise increase. Excessive exposure to environmental noise could potentially interfere with functional maturation of the auditory brain in developing individuals. The aim of the present study was to assess the association between exposure to annual average road traffic noise (LAeq) in schools and functional connectivity of key elements of the central auditory pathway in schoolchildren. A total of 229 children from 34 representative schools in the city of Barcelona with ages between 8 and 12 years (49.2% girls) were evaluated. LAeq was obtained as the mean of 2-consecutive day measurements inside classrooms before lessons started following standard procedures to obtain an indicator of long-term road traffic noise levels. A region-of-interest functional connectivity Magnetic Resonance Imaging (MRI) approach was adopted. Functional connectivity maps were generated for the inferior colliculus, medial geniculate body of the thalamus and primary auditory cortex as key levels of the central auditory pathway. Road traffic noise in schools was significantly associated with stronger connectivity between the inferior colliculus and a bilateral thalamic region adjacent to the medial geniculate body, and with stronger connectivity between the medial geniculate body and a bilateral brainstem region adjacent to the inferior colliculus. Such a functional connectivity strengthening effect did not extend to the cerebral cortex. The anatomy of the association implicating subcortical relays suggests that prolonged road traffic noise exposure in developing individuals may accelerate maturation in the basic elements of the auditory pathway. Future research is warranted to establish whether such a faster maturation in early pathway levels may ultimately reduce the developing potential in the whole auditory system.

Interaction between spatial perception and temporal perception enables preservation of cause-effect relationship: Visual psychophysics and neuronal dynamics

INTRODUCTION

Visual perception of moving objects is integral to our day-to-day life, integrating visual spatial and temporal perception. Most research studies have focused on finding the brain regions activated during motion perception. However, an empirically validated general mathematical model is required to understand the modulation of the motion perception. Here, we develop a mathematical formulation of the modulation of the perception of a moving object due to a change in speed, under the formulation of the invariance of causality.

METHODS

We formulated the perception of a moving object as the coordinate transformation from a retinotopic space onto perceptual space and derived a quantitative relationship between spatiotemporal coordinates. To validate our model, we undertook the analysis of two experiments: (i) the perceived length of the moving arc, and (ii) the perceived time while observing moving stimuli. We performed a magnetic resonance imaging (MRI) tractography investigation of subjects to demarcate the anatomical correlation of the modulation of the perception of moving objects.

RESULTS

Our theoretical model shows that the interaction between visual-spatial and temporal perception, during the perception of moving object is described by coupled linear equations; and experimental observations validate our model. We observed that cerebral area V5 may be an anatomical correlate for this interaction. The physiological basis of interaction is shown by a Lotka-Volterra system delineating interplay between acetylcholine and dopamine neurotransmitters, whose concentrations vary periodically with the orthogonal phase shift between them, occurring at the axodendritic synapse of complex cells at area V5.

CONCLUSION

Under the invariance of causality in the representation of events in retinotopic space and perceptual space, the speed modulates the perception of a moving object. This modulation may be due to variations of the tuning properties of complex cells at area V5 due to the dynamic interaction between acetylcholine and dopamine. Our analysis is the first significant study, to our knowledge, that establishes a mathematical linkage between motion perception and causality invariance.

BrainStat: A toolbox for brain-wide statistics and multimodal feature associations

Analysis and interpretation of neuroimaging datasets has become a multidisciplinary endeavor, relying not only on statistical methods, but increasingly on associations with respect to other brain-derived features such as gene expression, histological data, and functional as well as cognitive architectures. Here, we introduce BrainStat - a toolbox for (i) univariate and multivariate linear models in volumetric and surface-based brain imaging datasets, and (ii) multidomain feature association of results with respect to spatial maps of post-mortem gene expression and histology, task-based fMRI meta-analysis, as well as resting-state fMRI motifs across several common surface templates. The combination of statistics and feature associations into a turnkey toolbox streamlines analytical processes and accelerates cross-modal research. The toolbox is implemented in both Python and MATLAB, two widely used programming languages in the neuroimaging and neuroinformatics communities. BrainStat is openly available and complemented by an expandable documentation.

Effects of acute ischemic stroke on binaural perception

Stroke-induced lesions at different locations in the brain can affect various aspects of binaural hearing, including spatial perception. Previous studies found impairments in binaural hearing, especially in patients with temporal lobe tumors or lesions, but also resulting from lesions all along the auditory pathway from brainstem nuclei up to the auditory cortex. Currently, structural magnetic resonance imaging (MRI) is used in the clinical treatment routine of stroke patients. In combination with structural imaging, an analysis of binaural hearing enables a better understanding of hearing-related signaling pathways and of clinical disorders of binaural processing after a stroke. However, little data are currently available on binaural hearing in stroke patients, particularly for the acute phase of stroke. Here, we sought to address this gap in an exploratory study of patients in the acute phase of ischemic stroke. We conducted psychoacoustic measurements using two tasks of binaural hearing: binaural tone-in-noise detection, and lateralization of stimuli with interaural time- or level differences. The location of the stroke lesion was established by previously acquired MRI data. An additional general assessment included three-frequency audiometry, cognitive assessments, and depression screening. Fifty-five patients participated in the experiments, on average 5 days after their stroke onset. Patients whose lesions were in different locations were tested, including lesions in brainstem areas, basal ganglia, thalamus, temporal lobe, and other cortical and subcortical areas. Lateralization impairments were found in most patients with lesions within the auditory pathway. Lesioned areas at brainstem levels led to distortions of lateralization in both hemifields, thalamus lesions were correlated with a shift of the whole auditory space, whereas some cortical lesions predominantly affected the lateralization of stimuli contralateral to the lesion and resulted in more variable responses. Lateralization performance was also found to be affected by lesions of the right, but not the left, basal ganglia, as well as by lesions in non-auditory cortical areas. In general, altered lateralization was common in the stroke group. In contrast, deficits in tone-in-noise detection were relatively scarce in our sample of lesion patients, although a significant number of patients with multiple lesion sites were not able to complete the task.

Multimodal Assessment of Bottlenose Dolphin Auditory Nuclei Using 7-Tesla MRI, Immunohistochemistry and Stereology

The importance of assessing neurochemical processes in the cetacean brain as a tool for monitoring their cognitive health and to indirectly model human neurodegenerative conditions is increasingly evident, although available data are largely semiquantitative. High-resolution MRI for post-mortem brains and stereology allow for quantitative assessments of the cetacean brain. In this study, we scanned two brains of bottlenose dolphins in a 7-Tesla (7T) MR scanner and assessed the connectivity of the inferior colliculi and ventral cochlear nuclei using diffusion tensor imaging (DTI). Serial thick sections were investigated stereologically in one of the dolphins to generate rigorous quantitative estimates of identifiable cell types according to their morphology and expression of molecular markers, yielding reliable cell counts with most coefficients of error <10%. Fibronectin immunoreactivity in the dolphin resembled the pattern in a human chronic traumatic encephalopathy brain, suggesting that neurochemical compensation for insults such as hypoxia may constitute a noxious response in humans, while being physiological in dolphins. These data contribute to a growing body of knowledge on the morphological and neurochemical properties of the dolphin brain and highlight a stereological and neuroimaging workflow that may enable quantitative and translational assessment of pathological processes in the dolphin brain in the future.

Going down to go up: understanding human auditory cognition by investigating cortical-subcortical interactions

Human auditory cognition spans everything from detecting a creaking door in the night to enjoying beautiful music. Neurofunctional models of these processes tend to focus on cortical networks, but how do subcortical circuits contribute to auditory cognition? Answering this question will lead to a richer understanding of how we process the complex auditory world.

Selective attention sharpens population receptive fields in human auditory cortex

Selective attention enables the preferential processing of relevant stimulus aspects. Invasive animal studies have shown that attending a sound feature rapidly modifies neuronal tuning throughout the auditory cortex. Human neuroimaging studies have reported enhanced auditory cortical responses with selective attention. To date, it remains unclear how the results obtained with functional magnetic resonance imaging (fMRI) in humans relate to the electrophysiological findings in animal models. Here we aim to narrow the gap between animal and human research by combining a selective attention task similar in design to those used in animal electrophysiology with high spatial resolution ultra-high field fMRI at 7 Tesla. Specifically, human participants perform a detection task, whereas the probability of target occurrence varies with sound frequency. Contrary to previous fMRI studies, we show that selective attention resulted in population receptive field sharpening, and consequently reduced responses, at the attended sound frequencies. The difference between our results to those of previous fMRI studies supports the notion that the influence of selective attention on auditory cortex is diverse and may depend on context, stimulus, and task.

Micapipe: A pipeline for multimodal neuroimaging and connectome analysis

Multimodal magnetic resonance imaging (MRI) has accelerated human neuroscience by fostering the analysis of brain microstructure, geometry, function, and connectivity across multiple scales and in living brains. The richness and complexity of multimodal neuroimaging, however, demands processing methods to integrate information across modalities and to consolidate findings across different spatial scales. Here, we present micapipe, an open processing pipeline for multimodal MRI datasets. Based on BIDS-conform input data, micapipe can generate i) structural connectomes derived from diffusion tractography, ii) functional connectomes derived from resting-state signal correlations, iii) geodesic distance matrices that quantify cortico-cortical proximity, and iv) microstructural profile covariance matrices that assess inter-regional similarity in cortical myelin proxies. The above matrices can be automatically generated across established 18 cortical parcellations (100-1000 parcels), in addition to subcortical and cerebellar parcellations, allowing researchers to replicate findings easily across different spatial scales. Results are represented on three different surface spaces (native, conte69, fsaverage5), and outputs are BIDS-conform. Processed outputs can be quality controlled at the individual and group level. micapipe was tested on several datasets and is available at https://github.com/MICA-MNI/micapipe, documented at https://micapipe.readthedocs.io/, and containerized as a BIDS App http://bids-apps.neuroimaging.io/apps/. We hope that micapipe will foster robust and integrative studies of human brain microstructure, morphology, function, cand connectivity.

A specialized channel for encoding auditory transients in the magnocellular division of the human medial geniculate nucleus

We test the hypothesis that there exists a generalized magnocellular system in the brain optimized for temporal processing. In the visual system, it is well known that the magnocellular layers in the lateral geniculate nucleus (LGN) are strongly activated by transients and quickly habituate. However, little is known about the perhaps analogous magnocellular division of the medial geniculate nucleus (MGN), the auditory relay in the thalamus. We measured the functional responses of the MGN in 11 subjects who passively listened to sustained and transient nonlinguistic sounds, using functional MRI. We observed that voxels in the ventromedial portion of the MGN, corresponding to the magnocellular division, exhibited a robust preference to transient sounds, consistently across subjects, whereas the remainder of the MGN did not discriminate between sustained and transient sounds. We conclude that the magnocellular neurons in the MGN parallel the magnocellular neurons in its visual counterpart, LGN, and constitute an information stream specialized for encoding auditory dynamics.

A unified 3D map of microscopic architecture and MRI of the human brain

We present the first three-dimensional (3D) concordance maps of cyto- and fiber architecture of the human brain, combining histology, immunohistochemistry, and 7-T quantitative magnetic resonance imaging (MRI), in two individual specimens. These 3D maps each integrate data from approximately 800 microscopy sections per brain, showing neuronal and glial cell bodies, nerve fibers, and interneuronal populations, as well as ultrahigh-field quantitative MRI, all coaligned at the 200-μm scale to the stacked blockface images obtained during sectioning. These unprecedented 3D multimodal datasets are shared without any restrictions and provide a unique resource for the joint study of cell and fiber architecture of the brain, detailed anatomical atlasing, or modeling of the microscopic underpinnings of MRI contrasts.

Structural Connectivity of Human Inferior Colliculus Subdivisions Using in vivo and post mortem Diffusion MRI Tractography

Inferior colliculus (IC) is an obligatory station along the ascending auditory pathway that also has a high degree of top-down convergence via efferent pathways, making it a major computational hub. Animal models have attributed critical roles for the IC in in mediating auditory plasticity, egocentric selection, and noise exclusion. IC contains multiple functionally distinct subdivisions. These include a central nucleus that predominantly receives ascending inputs and external and dorsal nuclei that receive more heterogeneous inputs, including descending and multisensory connections. Subdivisions of human IC have been challenging to identify and quantify using standard brain imaging techniques such as MRI, and the connectivity of each of these subnuclei has not been identified in the human brain. In this study, we estimated the connectivity of human IC subdivisions with diffusion MRI (dMRI) tractography, using both anatomical-based seed analysis as well as unsupervised k-means clustering. We demonstrate sensitivity of tractography to overall IC connections in both high resolution post mortem and in vivo datasets. k-Means clustering of the IC streamlines in both the post mortem and in vivo datasets generally segregated streamlines based on their terminus beyond IC, such as brainstem, thalamus, or contralateral IC. Using fine-grained anatomical segmentations of the major IC subdivisions, the post mortem dataset exhibited unique connectivity patterns from each subdivision, including commissural connections through dorsal IC and lateral lemniscal connections to central and external IC. The subdivisions were less distinct in the context of in vivo connectivity, although lateral lemniscal connections were again highest to central and external IC. Overall, the unsupervised and anatomically driven methods provide converging evidence for distinct connectivity profiles for each of the IC subdivisions in both post mortem and in vivo datasets, suggesting that dMRI tractography with high quality data is sensitive to neural pathways involved in auditory processing as well as top-down control of incoming auditory information.

Cytoarchitectonic Maps of the Human Metathalamus in 3D Space

The human metathalamus plays an important role in processing visual and auditory information. Understanding its layers and subdivisions is important to gain insights in its function as a subcortical relay station and involvement in various pathologies. Yet, detailed histological references of the microanatomy in 3D space are still missing. We therefore aim at providing cytoarchitectonic maps of the medial geniculate body (MGB) and its subdivisions in the BigBrain – a high-resolution 3D-reconstructed histological model of the human brain, as well as probabilistic cytoarchitectonic maps of the MGB and lateral geniculate body (LGB). Therefore, histological sections of ten postmortem brains were studied. Three MGB subdivisions (MGBv, MGBd, MGBm) were identified on every 5th BigBrain section, and a deep-learning based tool was applied to map them on every remaining section. The maps were 3D-reconstructed to show the shape and extent of the MGB and its subdivisions with cellular precision. The LGB and MGB were additionally identified in nine other postmortem brains. Probabilistic cytoarchitectonic maps in the MNI “Colin27” and MNI ICBM152 reference spaces were computed which reveal an overall low interindividual variability in topography and extent. The probabilistic maps were included into the Julich-Brain atlas, and are freely available. They can be linked to other 3D data of human brain organization and serve as an anatomical reference for diagnostic, prognostic and therapeutic neuroimaging studies of healthy brains and patients. Furthermore, the high-resolution MGB BigBrain maps provide a basis for data integration, brain modeling and simulation to bridge the larger scale involvement of thalamocortical and local subcortical circuits.

Improving Imaging of the Brainstem and Cerebellum in Autistic Children: Transformation-Based High-Resolution Diffusion MRI (TiDi-Fused) in the Human Brainstem

Diffusion-weighted magnetic resonance imaging (dMRI) of the brainstem is technically challenging, especially in young autistic children as nearby tissue-air interfaces and motion (voluntary and physiological) can lead to artifacts. This limits the availability of high-resolution images, which are desirable for improving the ability to study brainstem structures. Furthermore, inherently low signal-to-noise ratios, geometric distortions, and sensitivity to motion not related to molecular diffusion have resulted in limited techniques for high-resolution data acquisition compared to other modalities such as T1-weighted imaging. Here, we implement a method for achieving increased apparent spatial resolution in pediatric dMRI that hinges on accurate geometric distortion correction and on high fidelity within subject image registration between dMRI and magnetization prepared rapid acquisition gradient echo (MPnRAGE) images. We call this post-processing pipeline T1 weighted-diffusion fused, or "TiDi-Fused". Data used in this work consists of dMRI data (2.4 mm resolution, corrected using FSL's Topup) and T1-weighted (T1w) MPnRAGE anatomical data (1 mm resolution) acquired from 128 autistic and non-autistic children (ages 6-10 years old). Accurate correction of geometric distortion permitted for a further increase in apparent resolution of the dMRI scan via boundary-based registration to the MPnRAGE T1w. Estimation of fiber orientation distributions and further analyses were carried out in the T1w space. Data processed with the TiDi-Fused method were qualitatively and quantitatively compared to data processed with conventional dMRI processing methods. Results show the advantages of the TiDi-Fused pipeline including sharper brainstem gray-white matter tissue contrast, improved inter-subject spatial alignment for group analyses of dMRI based measures, accurate spatial alignment with histology-based imaging of the brainstem, reduced variability in brainstem-cerebellar white matter tracts, and more robust biologically plausible relationships between age and brainstem-cerebellar white matter tracts. Overall, this work identifies a promising pipeline for achieving high-resolution imaging of brainstem structures in pediatric and clinical populations who may not be able to endure long scan times. This pipeline may serve as a gateway for feasibly elucidating brainstem contributions to autism and other conditions.

A methodological scoping review of the integration of fMRI to guide dMRI tractography. What has been done and what can be improved: A 20-year perspective

Combining MRI modalities is a growing trend in neurosciences. It provides opportunities to investigate the brain architecture supporting cognitive functions. Integrating fMRI activation to guide dMRI tractography offers potential advantages over standard tractography methods. A quick glimpse of the literature on this topic reveals that this technique is challenging, and no consensus or "best practices" currently exist, at least not within a single document. We present the first attempt to systematically analyze and summarize the literature of 80 studies that integrated task-based fMRI results to guide tractography, over the last two decades. We report 19 findings that cover challenges related to sample size, microstructure modelling, seeding methods, multimodal space registration, false negatives/positives, specificity/validity, gray/white matter interface and more. These findings will help the scientific community (1) understand the strengths and limitations of the approaches, (2) design studies using this integrative framework, and (3) motivate researchers to fill the gaps identified. We provide references toward best practices, in order to improve the overall result's replicability, sensitivity, specificity, and validity.

Altered processing of communication signals in the subcortical auditory sensory pathway in autism

Autism spectrum disorder (ASD) is characterised by social communication difficulties. These difficulties have been mainly explained by cognitive, motivational, and emotional alterations in ASD. The communication difficulties could, however, also be associated with altered sensory processing of communication signals. Here, we assessed the functional integrity of auditory sensory pathway nuclei in ASD in three independent functional magnetic resonance imaging experiments. We focused on two aspects of auditory communication that are impaired in ASD: voice identity perception, and recognising speech‐in‐noise. We found reduced processing in adults with ASD as compared to typically developed control groups (pairwise matched on sex, age, and full‐scale IQ) in the central midbrain structure of the auditory pathway (inferior colliculus [IC]). The right IC responded less in the ASD as compared to the control group for voice identity, in contrast to speech recognition. The right IC also responded less in the ASD as compared to the control group when passively listening to vocal in contrast to non‐vocal sounds. Within the control group, the left and right IC responded more when recognising speech‐in‐noise as compared to when recognising speech without additional noise. In the ASD group, this was only the case in the left, but not the right IC. The results show that communication signal processing in ASD is associated with reduced subcortical sensory functioning in the midbrain. The results highlight the importance of considering sensory processing alterations in explaining communication difficulties, which are at the core of ASD.

Using high spatial resolution fMRI to understand representation in the auditory network

Following rapid methodological advances, ultra-high field (UHF) functional and anatomical magnetic resonance imaging (MRI) has been repeatedly and successfully used for the investigation of the human auditory system in recent years. Here, we review this work and argue that UHF MRI is uniquely suited to shed light on how sounds are represented throughout the network of auditory brain regions. That is, the provided gain in spatial resolution at UHF can be used to study the functional role of the small subcortical auditory processing stages and details of cortical processing. Further, by combining high spatial resolution with the versatility of MRI contrasts, UHF MRI has the potential to localize the primary auditory cortex in individual hemispheres. This is a prerequisite to study how sound representation in higher-level auditory cortex evolves from that in early (primary) auditory cortex. Finally, the access to independent signals across auditory cortical depths, as afforded by UHF, may reveal the computations that underlie the emergence of an abstract, categorical sound representation based on low-level acoustic feature processing. Efforts on these research topics are underway. Here we discuss promises as well as challenges that come with studying these research questions using UHF MRI, and provide a future outlook.

Mapping cortico-subcortical sensitivity to 4 Hz amplitude modulation depth in human auditory system with functional MRI

Temporal modulations in the envelope of acoustic waveforms at rates around 4 Hz constitute a strong acoustic cue in speech and other natural sounds. It is often assumed that the ascending auditory pathway is increasingly sensitive to slow amplitude modulation (AM), but sensitivity to AM is typically considered separately for individual stages of the auditory system. Here, we used blood oxygen level dependent (BOLD) fMRI in twenty human subjects (10 male) to measure sensitivity of regional neural activity in the auditory system to 4 Hz temporal modulations. Participants were exposed to AM noise stimuli varying parametrically in modulation depth to characterize modulation-depth effects on BOLD responses. A Bayesian hierarchical modeling approach was used to model potentially nonlinear relations between AM depth and group-level BOLD responses in auditory regions of interest (ROIs). Sound stimulation activated the auditory brainstem and cortex structures in single subjects. BOLD responses to noise exposure in core and belt auditory cortices scaled positively with modulation depth. This finding was corroborated by whole-brain cluster-level inference. Sensitivity to AM depth variations was particularly pronounced in the Heschl's gyrus but also found in higher-order auditory cortical regions. None of the sound-responsive subcortical auditory structures showed a BOLD response profile that reflected the parametric variation in AM depth. The results are compatible with the notion that early auditory cortical regions play a key role in processing low-rate modulation content of sounds in the human auditory system.

Advances in spiral fMRI: A high-resolution study with single-shot acquisition

Spiral fMRI has been put forward as a viable alternative to rectilinear echo-planar imaging, in particular due to its enhanced average k-space speed and thus high acquisition efficiency. This renders spirals attractive for contemporary fMRI applications that require high spatiotemporal resolution, such as laminar or columnar fMRI. However, in practice, spiral fMRI is typically hampered by its reduced robustness and ensuing blurring artifacts, which arise from imperfections in both static and dynamic magnetic fields. Recently, these limitations have been overcome by the concerted application of an expanded signal model that accounts for such field imperfections, and its inversion by iterative image reconstruction. In the challenging ultra-high field environment of 7 Tesla, where field inhomogeneity effects are aggravated, both multi-shot and single-shot 2D spiral imaging at sub-millimeter resolution was demonstrated with high depiction quality and anatomical congruency. In this work, we further these advances towards a time series application of spiral readouts, namely, single-shot spiral BOLD fMRI at 0.8 mm in-plane resolution. We demonstrate that high-resolution spiral fMRI at 7 T is not only feasible, but delivers both excellent image quality, BOLD sensitivity, and spatial specificity of the activation maps, with little artifactual blurring. Furthermore, we show the versatility of the approach with a combined in/out spiral readout at a more typical resolution (1.5 mm), where the high acquisition efficiency allows to acquire two images per shot for improved sensitivity by echo combination.

Evaluation of the whole auditory pathway using high-resolution and functional MRI at 7T parallel-transmit

Purpose

The aim of the present study is to show a MR procedure for the evaluation of simultaneous left and right auditory functions with functional MRI, and high-resolution acquisition of anatomical auditory pathway using parallel-transmit (pTx) methods at 7T.

Methods

The time-efficient MR acquisition included two steps: RF weights were optimized for the regions-of-interest and high-resolution MR images of the inner-ear were acquired for the first 30 min (400 μm-iso resolution) followed by functional MRI acquisitions along the whole auditory pathway during the next 20 minutes. Data was processed with a linear cross-correlation analysis to define frequency preferences for each voxel in the auditory relays.

Results

Tonotopic maps revealed ordered bilateral frequency gradients in the auditory relays whereas at the level of the cochlear nuclei and superior olivary complexes the frequency gradients were less evident. A 21% increase in transmit-field efficiency was achieved over the left/right inner-ear regions and thus its main structures were clearly discernible using the pTx methods, compared to a single transmit RF coil.

Conclusion

Using 7T pTx allows a fast (less than 60 min in total) and qualitative evaluation of the simultaneous left and right auditory response along the entire auditory pathway, together with high-resolution anatomical images of the inner-ear. This could be further used for patient examination at 7T.

The BigBrainWarp toolbox for integration of BigBrain 3D histology with multimodal neuroimaging

Neuroimaging stands to benefit from emerging ultrahigh-resolution 3D histological atlases of the human brain; the first of which is 'BigBrain'. Here, we review recent methodological advances for the integration of BigBrain with multi-modal neuroimaging and introduce a toolbox, 'BigBrainWarp', that combines these developments. The aim of BigBrainWarp is to simplify workflows and support the adoption of best practices. This is accomplished with a simple wrapper function that allows users to easily map data between BigBrain and standard MRI spaces. The function automatically pulls specialised transformation procedures, based on ongoing research from a wide collaborative network of researchers. Additionally, the toolbox improves accessibility of histological information through dissemination of ready-to-use cytoarchitectural features. Finally, we demonstrate the utility of BigBrainWarp with three tutorials and discuss the potential of the toolbox to support multi-scale investigations of brain organisation.

A high-resolution interactive atlas of the human brainstem using magnetic resonance imaging

Conventional atlases of the human brainstem are limited by the inflexible, sparsely-sampled, two-dimensional nature of histology, or the low spatial resolution of conventional magnetic resonance imaging (MRI). Postmortem high-resolution MRI circumvents the challenges associated with both modalities. A single human brainstem specimen extending from the rostral diencephalon through the caudal medulla was prepared for imaging after the brain was removed from a 65-year-old male within 24 h of death. The specimen was formalin-fixed for two weeks, then rehydrated and placed in a custom-made MRI compatible tube and immersed in liquid fluorocarbon. MRI was performed in a 7-Tesla scanner with 120 unique diffusion directions. Acquisition time for anatomic and diffusion images were 14 h and 208 h, respectively. Segmentation was performed manually. Deterministic fiber tractography was done using strategically chosen regions of interest and avoidance, with manual editing using expert knowledge of human neuroanatomy. Anatomic and diffusion images were rendered with isotropic resolutions of 50 μm and 200 μm, respectively. Ninety different structures were segmented and labeled, and 11 different fiber bundles were rendered with tractography. The complete atlas is available online for interactive use at https://www.civmvoxport.vm.duke.edu/voxbase/login.php?return_url=%2Fvoxbase%2F. This atlas presents multiple contrasting datasets and selected tract reconstruction with unprecedented resolution for MR imaging of the human brainstem. There are immediate applications in neuroanatomical education, with the potential to serve future applications for neuroanatomical research and enhanced neurosurgical planning through "safe" zones of entry into the human brainstem.

In Vivo Super-Resolution Track-Density Imaging for Thalamic Nuclei Identification

The development of novel techniques for the in vivo, non-invasive visualization and identification of thalamic nuclei has represented a major challenge for human neuroimaging research in the last decades. Thalamic nuclei have important implications in various key aspects of brain physiology and many of them show selective alterations in various neurologic and psychiatric disorders. In addition, both surgical stimulation and ablation of specific thalamic nuclei have been proven to be useful for the treatment of different neuropsychiatric diseases. The present work aimed at describing a novel protocol for histologically guided delineation of thalamic nuclei based on short-tracks track-density imaging (stTDI), which is an advanced imaging technique exploiting high angular resolution diffusion tractography to obtain super-resolved white matter maps. We demonstrated that this approach can identify up to 13 distinct thalamic nuclei bilaterally with very high inter-subject (ICC: 0.996, 95% CI: 0.993-0.998) and inter-rater (ICC:0.981; 95% CI:0.963-0.989) reliability, and that both subject-based and group-level thalamic parcellation show a fair share of similarity to a recent standard-space histological thalamic atlas. Finally, we showed that stTDI-derived thalamic maps can be successfully employed to study structural and functional connectivity of the thalamus and may have potential implications both for basic and translational research, as well as for presurgical planning purposes.

Ex Vivo MR Histology and Cytometric Feature Mapping Connect Three-dimensional in Vivo MR Images to Two-dimensional Histopathologic Images of Murine Sarcomas

Purpose To establish a platform for quantitative tissue-based interpretation of cytoarchitecture features from tumor MRI measurements. Materials and Methods In a pilot preclinical study, multicontrast in vivo MRI of murine soft-tissue sarcomas in 10 mice, followed by ex vivo MRI of fixed tissues (termed MR histology), was performed. Paraffin-embedded limb cross-sections were stained with hematoxylin-eosin, digitized, and registered with MRI. Registration was assessed by using binarized tumor maps and Dice similarity coefficients (DSCs). Quantitative cytometric feature maps from histologic slides were derived by using nuclear segmentation and compared with registered MRI, including apparent diffusion coefficients and transverse relaxation times as affected by magnetic field heterogeneity (T2* maps). Cytometric features were compared with each MR image individually by using simple linear regression analysis to identify the features of interest, and the goodness of fit was assessed on the basis of R² values. Results Registration of MR images to histopathologic slide images resulted in mean DSCs of 0.912 for ex vivo MR histology and 0.881 for in vivo MRI. Triplicate repeats showed high registration repeatability (mean DSC, >0.9). Whole-slide nuclear segmentations were automated to detect nuclei on histopathologic slides (DSC = 0.8), and feature maps were generated for correlative analysis with MR images. Notable trends were observed between cell density and in vivo apparent diffusion coefficients (best line fit: R² = 0.96, P < .001). Multiple cytoarchitectural features exhibited linear relationships with in vivo T2* maps, including nuclear circularity (best line fit: R² = 0.99, P < .001) and variance in nuclear circularity (best line fit: R² = 0.98, P < .001). Conclusion An infrastructure for registering and quantitatively comparing in vivo tumor MRI with traditional histologic analysis was successfully implemented in a preclinical pilot study of soft-tissue sarcomas. Keywords: MRI, Pathology, Animal Studies, Tissue Characterization Supplemental material is available for this article. © RSNA, 2021.

Multi-contrast anatomical subcortical structures parcellation

The human subcortex is comprised of more than 450 individual nuclei which lie deep in the brain. Due to their small size and close proximity, up until now only 7% have been depicted in standard MRI atlases. Thus, the human subcortex can largely be considered as terra incognita. Here, we present a new open-source parcellation algorithm to automatically map the subcortex. The new algorithm has been tested on 17 prominent subcortical structures based on a large quantitative MRI dataset at 7 Tesla. It has been carefully validated against expert human raters and previous methods, and can easily be extended to other subcortical structures and applied to any quantitative MRI dataset. In sum, we hope this novel parcellation algorithm will facilitate functional and structural neuroimaging research into small subcortical nuclei and help to chart terra incognita.

Depth relationships and measures of tissue thickness in dorsal midbrain

Dorsal human midbrain contains two nuclei with clear laminar organization, the superior and inferior colliculi. These nuclei extend in depth between the superficial dorsal surface of midbrain and a deep midbrain nucleus, the periaqueductal gray matter (PAG). The PAG, in turn, surrounds the cerebral aqueduct (CA). This study examined the use of two depth metrics to characterize depth and thickness relationships within dorsal midbrain using the superficial surface of midbrain and CA as references. The first utilized nearest-neighbor Euclidean distance from one reference surface, while the second used a level-set approach that combines signed distances from both reference surfaces. Both depth methods provided similar functional depth profiles generated by saccadic eye movements in a functional MRI task, confirming their efficacy for delineating depth for superficial functional activity. Next, the boundaries of the PAG were estimated using Euclidean distance together with elliptical fitting, indicating that the PAG can be readily characterized by a smooth surface surrounding PAG. Finally, we used the level-set approach to measure tissue depth between the superficial surface and the PAG, thus characterizing the variable thickness of the colliculi. Overall, this study demonstrates depth-mapping schemes for human midbrain that enables accurate segmentation of the PAG and consistent depth and thickness estimates of the superior and inferior colliculi.

Abstract rules drive adaptation in the subcortical sensory pathway

The subcortical sensory pathways are the fundamental channels for mapping the outside world to our minds. Sensory pathways efficiently transmit information by adapting neural responses to the local statistics of the sensory input. The long-standing mechanistic explanation for this adaptive behaviour is that neural activity decreases with increasing regularities in the local statistics of the stimuli. An alternative account is that neural coding is directly driven by expectations of the sensory input. Here, we used abstract rules to manipulate expectations independently of local stimulus statistics. The ultra-high-field functional-MRI data show that abstract expectations can drive the response amplitude to tones in the human auditory pathway. These results provide first unambiguous evidence of abstract processing in a subcortical sensory pathway. They indicate that the neural representation of the outside world is altered by our prior beliefs even at initial points of the processing hierarchy.

What are the gray and white matter volumes of the human spinal cord?

The gray matter of the spinal cord is the seat of somata of various types of neurons devoted to the sensory and motor activities of the limbs and trunk as well as a part of the autonomic nervous system. The volume of the spinal gray matter is an indicator of the local neuronal processing, and this can decrease due to atrophy associated with degenerative diseases and injury. Nevertheless, the absolute volume of the human spinal cord has rarely been reported, if ever. Here, we use high-resolution magnetic resonance imaging, with a cross-sectional resolution of 50 × 50 μm and a voxel size of 0.0005 mm³ to estimate the total gray and white matter volume of a post mortem human female spinal cord. Segregation of gray and white matter was accomplished using deep learning image segmentation. Furthermore, we include data from a male spinal cord of a previously published study. The gray and white matter volumes were found to be 2.87 and 11.33 mL, respectively, for the female and 3.55 and 19.33 mL, respectively, for the male. The gray and white matter profiles along the vertebral axis were found to be strikingly similar, and the volumes of the cervical, thoracic, and lumbosacral sections were almost equal.NEW & NOTEWORTHY Here, we combine high-field MRI (9.4 T) and deep learning for a post mortem reconstruction of the gray and white matter in human spinal cords. We report a minuscule total gray matter volume of 2.87 mL for a female and 3.55 mL for a male. For comparison, these volumes correspond approximately to the distal digit of the little finger.

Exercise alters brain activation in Gulf War Illness and Myalgic Encephalomyelitis/Chronic Fatigue Syndrome

Gulf War Illness affects 25-30% of American veterans deployed to the 1990-91 Persian Gulf War and is characterized by cognitive post-exertional malaise following physical effort. Gulf War Illness remains controversial since cognitive post-exertional malaise is also present in the more common Myalgic Encephalomyelitis/Chronic Fatigue Syndrome. An objective dissociation between neural substrates for cognitive post-exertional malaise in Gulf War Illness and Myalgic Encephalomyelitis/Chronic Fatigue Syndrome would represent a biological basis for diagnostically distinguishing these two illnesses. Here, we used functional magnetic resonance imaging to measure neural activity in healthy controls and patients with Gulf War Illness and Myalgic Encephalomyelitis/Chronic Fatigue Syndrome during an N-back working memory task both before and after exercise. Whole brain activation during working memory (2-Back > 0-Back) was equal between groups prior to exercise. Exercise had no effect on neural activity in healthy controls yet caused deactivation within dorsal midbrain and cerebellar vermis in Gulf War Illness relative to Myalgic Encephalomyelitis/Chronic Fatigue Syndrome patients. Further, exercise caused increased activation among Myalgic Encephalomyelitis/Chronic Fatigue Syndrome patients within the dorsal midbrain, left operculo-insular cortex (Rolandic operculum) and right middle insula. These regions-of-interest underlie threat assessment, pain, interoception, negative emotion and vigilant attention. As they only emerge post-exercise, these regional differences likely represent neural substrates of cognitive post-exertional malaise useful for developing distinct diagnostic criteria for Gulf War Illness and Myalgic Encephalomyelitis/Chronic Fatigue Syndrome.

Improving a probabilistic cytoarchitectonic atlas of auditory cortex using a novel method for inter-individual alignment

The human superior temporal plane, the site of the auditory cortex, displays high inter-individual macro-anatomical variation. This questions the validity of curvature-based alignment (CBA) methods for in vivo imaging data. Here, we have addressed this issue by developing CBA+, which is a cortical surface registration method that uses prior macro-anatomical knowledge. We validate this method by using cytoarchitectonic areas on 10 individual brains (which we make publicly available). Compared to volumetric and standard surface registration, CBA+ results in a more accurate cytoarchitectonic auditory atlas. The improved correspondence of micro-anatomy following the improved alignment of macro-anatomy validates the superiority of CBA+ compared to CBA. In addition, we use CBA+ to align in vivo and postmortem data. This allows projection of functional and anatomical information collected in vivo onto the cytoarchitectonic areas, which has the potential to contribute to the ongoing debate on the parcellation of the human auditory cortex.

Investigating the Reliability of Population Receptive Field Size Estimates Using fMRI

In functional MRI (fMRI), population receptive field (pRF) models allow a quantitative description of the response as a function of the features of the stimuli that are relevant for each voxel. The most popular pRF model used in fMRI assumes a Gaussian shape in the features space (e.g., the visual field) reducing the description of the voxel’s pRF to the Gaussian mean (the pRF preferred feature) and standard deviation (the pRF size). The estimation of the pRF mean has been proven to be highly reliable. However, the estimate of the pRF size has been shown not to be consistent within and between subjects. While this issue has been noted experimentally, here we use an optimization theory perspective to describe how the inconsistency in estimating the pRF size is linked to an inherent property of the Gaussian pRF model. When fitting such models, the goodness of fit is less sensitive to variations in the pRF size than to variations in the pRF mean. We also show how the same issue can be considered from a bias-variance perspective. We compare different estimation procedures in terms of the reliability of their estimates using simulated and real fMRI data in the visual (using the Human Connectome Project database) and auditory domain. We show that, the reliability of the estimate of the pRF size can be improved considering a linear combination of those pRF models with similar goodness of fit or a permutation based approach. This increase in reliability of the pRF size estimate does not affect the reliability of the estimate of the pRF mean and the prediction accuracy.

In vivo Probabilistic Structural Atlas of the Inferior and Superior Colliculi, Medial and Lateral Geniculate Nuclei and Superior Olivary Complex in Humans Based on 7 Tesla MRI

Despite extensive neuroimaging research of primary sensory cortices involved in auditory and visual functions, subcortical structures within these domains, such as the inferior and superior colliculi, the medial and lateral geniculate nuclei and the superior olivary complex, are currently understudied with magnetic resonance imaging (MRI) in living humans. This is because a precise localization of these nuclei is hampered by the limited contrast and sensitivity of conventional neuroimaging methods for deep brain nuclei. In this work, we used 7 Tesla multi-modal (T₂-weighted and diffusion fractional anisotropy) 1.1 mm isotropic resolution MRI to achieve high sensitivity and contrast for single-subject brainstem and thalamic nuclei delineation. After precise coregistration to stereotactic space, we generated an in vivo human probabilistic atlas of auditory (medial geniculate nucleus, inferior colliculus, and superior olivary complex) and visual (lateral geniculate nucleus and superior colliculus) subcortical nuclei. We foresee the use of this atlas as a tool to precisely identify the location and shape of auditory/visual deep nuclei in research as well as clinical human studies.