A global optimisation method for robust affine registration of brain images

Motor or non-motor speech interference? A multimodal fMRI and direct cortical stimulation mapping study

We retrospectively analyzed data from 15 patients, with a normal pre-operative cognitive performance, undergoing awake surgery for left fronto-temporal low-grade glioma. We combined a pre-surgical measure (fMRI maps of motor- and language-related centers) with intra-surgical measures (MNI-registered cortical sites data obtained during intra-operative direct electrical stimulation, DES, while they performed the two most common language tasks: number counting and picture naming). Selective DES effects along the precentral gyrus/inferior frontal gyrus (and/or the connected speech articulation network) were obtained. DES of the precentral gyrus evoked the motor speech arrest, i.e., anarthria (with apparent mentalis muscle movements). We calculated the number of shared voxels between the lip-tongue and overt counting related- and silent naming-related fMRI maps and the Volumes of Interest (VOIs) obtained by merging together the MNI sites at which a given speech disturbance was observed, normalized on their mean the values (i.e., Z score). Both tongue- and lips-related movements fMRI maps maximally overlapped (Z = 1.05 and Z = 0.94 for lips and tongue vs. 0.16 and -1.003 for counting and naming) with the motor speech arrest seed. DES of the inferior frontal gyrus, pars opercularis and the rolandic operculum induced speech arrest proper (without apparent mentalis muscle movements). This area maximally overlapped with overt counting-related fMRI map (Z = -0.11 and Z = 0.09 for lips and tongue vs. 0.9 and 0.0006 for counting and naming). Interestingly, our fMRI maps indicated reduced Broca's area activity during silent speech compared to overt speech. Lastly, DES of the inferior frontal gyrus, pars opercularis and triangularis evoked variations of the output, i.e., dysarthria, a motor speech disorder occurring when patients cannot control the muscles used to produce articulated sounds (phonemes). Silent object naming-related fMRI map maximally overlapped (Z = -0.93 and Z = -1.04 for lips and tongue vs. -1.07 and 0.99 for counting and naming) with this seed. Speech disturbances evoked by DES may be thought of as selective interferences with specific recruitment of left inferior frontal gyrus and precentral cortex which are differentiable in terms of the specific interference induced.

Behavioral Bias for Exploration Is Associated with Enhanced Signaling in the Lateral and Medial Frontopolar Cortex

Should we keep doing what we know works for us, or should we risk trying something new as it could work even better? The exploration-exploitation dilemma is ubiquitous in daily life decision-making, and balancing between the two is crucial for adaptive behavior. Yet, we only have started to unravel the neurocognitive mechanisms that help us to find this balance in practice. Analyzing BOLD signals of healthy young adults during virtual foraging, we could show that a behavioral tendency for prolonged exploitation was associated with weakened signaling during exploration in central node points of the frontoparietal attention network, plus the frontopolar cortex. These results provide an important link between behavioral heuristics that we use to balance between exploitation and exploration and the brain function that supports shifts from one tendency to the other. Importantly, they stress that interindividual differences in behavioral strategies are reflected in differences in brain activity during exploration and should thus be more in the focus of basic research that aims at delineating general laws governing visual attention.

Neural Tracking of Perceived Parent, but Not Peer, Norms Is Associated with Longitudinal Changes in Adolescent Attitudes about Externalizing Behaviors

Adolescents' perceptions of parent and peer norms about externalizing behaviors influence the extent to which they adopt similar attitudes, yet little is known about how the trajectories of perceived parent and peer norms are related to trajectories of personal attitudes across adolescence. Neural development of midline regions implicated in self-other processing may underlie developmental changes in parent and peer influence. Here, we examined whether neural processing of perceived parent and peer norms in midline regions during self-evaluations would be associated with trajectories of personal attitudes about externalizing behaviors. Trajectories of adolescents' perceived parent and peer norms were examined longitudinally with functional neuroimaging (n = 165; ages 11-16 years across three waves; 86 girls, 79 boys; 29.7% White, 21.8% Black, 35.8% Latinx, 12.7% other/multiracial). Behavioral results showed perceived parent norms were less permissive than adolescents' own attitudes about externalizing behaviors, whereas perceived peer norms were more permissive than adolescents' own attitudes, effects that increased from early to middle adolescence. Although younger adolescents reported less permissive attitudes when they spontaneously tracked perceived parent norms in the ventromedial and medial pFCs during self-evaluations, this effect weakened as they aged. No brain-behavior effects were found when tracking perceived peer norms. These findings elucidate how perceived parent and peer norms change in parallel with personal attitudes about externalizing behaviors from early to middle adolescence and underscore the importance of spontaneous neural tracking of perceived parent norms during self-evaluations for buffering permissive personal attitudes, particularly in early adolescence.

Joint MAPLE: Accelerated joint T1 and T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ mapping with scan-specific self-supervised networks

PURPOSE

Quantitative MRI finds important applications in clinical and research studies. However, it is encoding intensive and may suffer from prohibitively long scan times. Accelerated MR parameter mapping techniques have been developed to help address these challenges. Here, an accelerated joint T1,T 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ , frequency and proton density mapping technique with scan-specific self-supervised network reconstruction is proposed to synergistically combine parallel imaging, model-based, and deep learning approaches to speed up parameter mapping.

METHODS

Proposed framework, Joint MAPLE, includes parallel imaging, signal modeling, and data consistency blocks which are optimized jointly in a combined loss function. A scan-specific self-supervised reconstruction is embedded into the framework, which takes advantage of multi-contrast data from a multi-echo, multi-flip angle, gradient echo acquisition.

RESULTS

In comparison with parallel reconstruction techniques powered by low-rank methods, emerging scan specific networks, and model-basedT 2 * $$ {{\mathrm{T}}_2}^{\ast } $$ estimation approaches, the proposed framework reduces the reconstruction error in parameter maps by approximately two-fold on average at acceleration rates as high as R = 16 with uniform sampling. It can outperform evaluated parallel reconstruction techniques up to four-fold on average in the presence of challenging sub-sampling masks. It is observed that Joint MAPLE performs well at extreme acceleration rates of R = 25 and R = 36 with error values less than 20%.

CONCLUSION

Joint MAPLE enables higher fidelity parameter estimation at high acceleration rates by synergistically combining parallel imaging and model-based parameter mapping and exploiting multi-echo, multi-flip angle datasets. Utilizing a scan-specific self-supervised reconstruction obviates the need for large data sets for training while improving the parameter estimation ability.

Characteristic BOLD signals are detectable in white matter of the spinal cord at rest and after a stimulus

We report the reliable detection of reproducible patterns of blood-oxygenation-level-dependent (BOLD) MRI signals within the white matter (WM) of the spinal cord during a task and in a resting state. Previous functional MRI studies have shown that BOLD signals are robustly detectable not only in gray matter (GM) in the brain but also in cerebral WM as well as the GM within the spinal cord, but similar signals in WM of the spinal cord have been overlooked. In this study, we detected BOLD signals in the WM of the spinal cord in squirrel monkeys and studied their relationships with the locations and functions of ascending and descending WM tracts. Tactile sensory stimulus -evoked BOLD signal changes were detected in the ascending tracts of the spinal cord using a general-linear model. Power spectral analysis confirmed that the amplitude at the fundamental frequency of the response to a periodic stimulus was significantly higher in the ascending tracts than the descending ones. Independent component analysis of resting-state signals identified coherent fluctuations from eight WM hubs which correspond closely to the known anatomical locations of the major WM tracts. Resting-state analyses showed that the WM hubs exhibited correlated signal fluctuations across spinal cord segments in reproducible patterns that correspond well with the known neurobiological functions of WM tracts in the spinal cord. Overall, these findings provide evidence of a functional organization of intraspinal WM tracts and confirm that they produce hemodynamic responses similar to GM both at baseline and under stimulus conditions.

Life-course neighbourhood deprivation and brain structure in older adults: the Lothian Birth Cohort 1936

Neighbourhood disadvantage may be associated with brain health but the importance of exposure at different stages of the life course is poorly understood. Utilising the Lothian Birth Cohort 1936, we explored the relationship between residential neighbourhood deprivation from birth to late adulthood, and global and local neuroimaging measures at age 73. A total of 689 participants had at least one valid brain measures (53% male); to maximise the sample size structural equation models with full information maximum likelihood were conducted. Residing in disadvantaged neighbourhoods in mid- to late adulthood was associated with smaller total brain (β = -0.06; SE = 0.02; sample size[N] = 658; number of pairwise complete observations[n]=390), grey matter (β = -0.11; SE = 0.03; N = 658; n = 390), and normal-appearing white matter volumes (β = -0.07; SE = 0.03; N = 658; n = 390), thinner cortex (β = -0.14; SE = 0.06; N = 636; n = 379), and lower general white matter fractional anisotropy (β = -0.19; SE = 0.06; N = 665; n = 388). We also found some evidence on the accumulating impact of neighbourhood deprivation from birth to late adulthood on age 73 total brain (β = -0.06; SE = 0.02; N = 658; n = 276) and grey matter volumes (β = -0.10; SE = 0.04; N = 658; n = 276). Local analysis identified affected focal cortical areas and specific white matter tracts. Among individuals belonging to lower social classes, the brain-neighbourhood associations were particularly strong, with the impact of neighbourhood deprivation on total brain and grey matter volumes, and general white matter fractional anisotropy accumulating across the life course. Our findings suggest that living in deprived neighbourhoods across the life course, but especially in mid- to late adulthood, is associated with adverse brain morphologies, with lower social class amplifying the vulnerability.

Microstructure predicts non-motor outcomes following deep brain stimulation in Parkinson's disease

Deep brain stimulation of the subthalamic nucleus (STN-DBS) effectively treats motor and non-motor symptoms in advanced Parkinson's disease (PD). As considerable interindividual variability of outcomes exists, neuroimaging-based biomarkers, including microstructural metrics, have been proposed to anticipate treatment response. In this prospective open-label study, we sought to detect microstructural properties of brain areas associated with short-term non-motor outcomes following STN-DBS. Thirty-seven PD patients underwent diffusion MRI and clinical assessments at preoperative baseline and 6-month follow-up. Whole brain voxel-wise analysis assessed associations between microstructural metrics and non-motor outcomes. Intact microstructure within specific areas, including the right insular cortex, right putamen, right cingulum, and bilateral corticospinal tract were associated with greater postoperative improvement of non-motor symptom burden. Furthermore, microstructural properties of distinct brain regions were associated with postoperative changes in sleep, attention/memory, urinary symptoms, and apathy. In conclusion, diffusion MRI could support preoperative patient counselling by identifying patients with above- or below-average non-motor responses.

Distinctive and Complementary Roles of Default Mode Network Subsystems in Semantic Cognition

The default mode network (DMN) typically deactivates to external tasks, yet supports semantic cognition. It comprises medial temporal (MT), core, and frontotemporal (FT) subsystems, but its functional organization is unclear: the requirement for perceptual coupling versus decoupling, input modality (visual/verbal), type of information (social/spatial), and control demands all potentially affect its recruitment. We examined the effect of these factors on activation and deactivation of DMN subsystems during semantic cognition, across four task-based human functional magnetic resonance imaging (fMRI) datasets, and localized these responses in whole-brain state space defined by gradients of intrinsic connectivity. FT showed activation consistent with a central role across domains, tasks, and modalities, although it was most responsive to abstract, verbal tasks; this subsystem uniquely showed more "tuned" states characterized by increases in both activation and deactivation when semantic retrieval demands were higher. MT also activated to both perceptually coupled (scenes) and decoupled (autobiographical memory) tasks and showed stronger responses to picture associations, consistent with a role in scene construction. Core DMN consistently showed deactivation, especially to externally oriented tasks. These diverse contributions of DMN subsystems to semantic cognition were related to their location on intrinsic connectivity gradients: activation was closer to the sensory-motor cortex than deactivation, particularly for FT and MT, while activation for core DMN was distant from both visual cortex and cognitive control. These results reveal distinctive yet complementary DMN responses: MT and FT support different memory-based representations that are accessed externally and internally, while deactivation in core DMN is associated with demanding, external semantic tasks.

Premotor cortical beta synchronization and the network neuromodulation of externally paced finger tapping in Parkinson's disease

Parkinson's disease (PD) is characterized by the disruption of repetitive, concurrent and sequential motor actions due to compromised timing-functions principally located in cortex-basal ganglia (BG) circuits. Increasing evidence suggests that motor impairments in untreated PD patients are linked to an excessive synchronization of cortex-BG activity at beta frequencies (13-30 Hz). Levodopa and subthalamic nucleus deep brain stimulation (STN-DBS) suppress pathological beta-band reverberation and improve the motor symptoms in PD. Yet a dynamic tuning of beta oscillations in BG-cortical loops is fundamental for movement-timing and synchronization, and the impact of PD therapies on sensorimotor functions relying on neural transmission in the beta frequency-range remains controversial. Here, we set out to determine the differential effects of network neuromodulation through dopaminergic medication (ON and OFF levodopa) and STN-DBS (ON-DBS, OFF-DBS) on tapping synchronization and accompanying cortical activities. To this end, we conducted a rhythmic finger-tapping study with high-density EEG-recordings in 12 PD patients before and after surgery for STN-DBS and in 12 healthy controls. STN-DBS significantly ameliorated tapping parameters as frequency, amplitude and synchrony to the given auditory rhythms. Aberrant neurophysiologic signatures of sensorimotor feedback in the beta-range were found in PD patients: their neural modulation was weaker, temporally sluggish and less distributed over the right cortex in comparison to controls. Levodopa and STN-DBS boosted the dynamics of beta-band modulation over the right hemisphere, hinting to an improved timing of movements relying on tactile feedback. The strength of the post-event beta rebound over the supplementary motor area correlated significantly with the tapping asynchrony in patients, thus indexing the sensorimotor match between the external auditory pacing signals and the performed taps. PD patients showed an excessive interhemispheric coherence in the beta-frequency range during the finger-tapping task, while under DBS-ON the cortico-cortical connectivity in the beta-band was normalized. Ultimately, therapeutic DBS significantly ameliorated the auditory-motor coupling of PD patients, enhancing the electrophysiological processing of sensorimotor feedback-information related to beta-band activity, and thus allowing a more precise cued-tapping performance.

Mechanical Properties of White Matter Tracts in Aging Assessed via Anisotropic MR Elastography

Magnetic resonance elastography (MRE) is a promising neuroimaging technique to probe tissue microstructure, which has revealed widespread softening with loss of structural integrity in the aging brain. Traditional MRE approaches assume mechanical isotropy. However, white matter is known to be anisotropic from aligned, myelinated axonal bundles, which can lead to uncertainty in mechanical property estimates in these areas when using isotropic MRE. Recent advances in anisotropic MRE now allow for estimation of shear and tensile anisotropy, along with substrate shear modulus, in white matter tracts. The objective of this study was to investigate age-related differences in anisotropic mechanical properties in human brain white matter tracts for the first time. Anisotropic mechanical properties in all tracts were found to be significantly lower in older adults compared to young adults, with average property differences ranging between 0.028-0.107 for shear anisotropy and between 0.139-0.347 for tensile anisotropy. Stiffness perpendicular to the axonal fiber direction was also significantly lower in older age, but only in certain tracts. When compared with fractional anisotropy measures from diffusion tensor imaging, we found that anisotropic MRE measures provided additional, complementary information in describing differences between the white matter integrity of young and older populations. Anisotropic MRE provides a new tool for studying white matter structural integrity in aging and neurodegeneration.

Counterfactual thinking induces different neural patterns of memory modification in anxious individuals

Episodic counterfactual thinking (eCFT) is the process of mentally simulating alternate versions of experiences, which confers new phenomenological properties to the original memory and may be a useful therapeutic target for trait anxiety. However, it remains unclear how the neural representations of a memory change during eCFT. We hypothesized that eCFT-induced memory modification is associated with changes to the neural pattern of a memory primarily within the default mode network, moderated by dispositional anxiety levels. We tested this proposal by examining the representational dynamics of eCFT for 39 participants varying in trait anxiety. During eCFT, lateral parietal regions showed progressively more distinct activity patterns, whereas medial frontal neural activity patterns became more similar to those of the original memory. Neural pattern similarity in many default mode network regions was moderated by trait anxiety, where highly anxious individuals exhibited more generalized representations for upward eCFT (better counterfactual outcomes), but more distinct representations for downward eCFT (worse counterfactual outcomes). Our findings illustrate the efficacy of examining eCFT-based memory modification via neural pattern similarity, as well as the intricate interplay between trait anxiety and eCFT generation.

Constructing personalized characterizations of structural brain aberrations in patients with dementia using explainable artificial intelligence

Deep learning approaches for clinical predictions based on magnetic resonance imaging data have shown great promise as a translational technology for diagnosis and prognosis in neurological disorders, but its clinical impact has been limited. This is partially attributed to the opaqueness of deep learning models, causing insufficient understanding of what underlies their decisions. To overcome this, we trained convolutional neural networks on structural brain scans to differentiate dementia patients from healthy controls, and applied layerwise relevance propagation to procure individual-level explanations of the model predictions. Through extensive validations we demonstrate that deviations recognized by the model corroborate existing knowledge of structural brain aberrations in dementia. By employing the explainable dementia classifier in a longitudinal dataset of patients with mild cognitive impairment, we show that the spatially rich explanations complement the model prediction when forecasting transition to dementia and help characterize the biological manifestation of disease in the individual brain. Overall, our work exemplifies the clinical potential of explainable artificial intelligence in precision medicine.

Quantifying 3D MR fingerprinting (3D-MRF) reproducibility across subjects, sessions, and scanners automatically using MNI atlases

PURPOSE

Quantitative MRI techniques such as MR fingerprinting (MRF) promise more objective and comparable measurements of tissue properties at the point-of-care than weighted imaging. However, few direct cross-modal comparisons of MRF's repeatability and reproducibility versus weighted acquisitions have been performed. This work proposes a novel fully automated pipeline for quantitatively comparing cross-modal imaging performance in vivo via atlas-based sampling.

METHODS

We acquire whole-brain 3D-MRF, turbo spin echo, and MPRAGE sequences three times each on two scanners across 10 subjects, for a total of 60 multimodal datasets. The proposed automated registration and analysis pipeline uses linear and nonlinear registration to align all qualitative and quantitative DICOM stacks to Montreal Neurological Institute (MNI) 152 space, then samples each dataset's native space through transformation inversion to compare performance within atlas regions across subjects, scanners, and repetitions.

RESULTS

Voxel values within MRF-derived maps were found to be more repeatable (σT1  = 1.90, σT2  = 3.20) across sessions than vendor-reconstructed MPRAGE (σT1w  = 6.04) or turbo spin echo (σT2w  = 5.66) images. Additionally, MRF was found to be more reproducible across scanners (σT1  = 2.21, σT2  = 3.89) than either qualitative modality (σT1w  = 7.84, σT2w  = 7.76). Notably, differences between repeatability and reproducibility of in vivo MRF were insignificant, unlike the weighted images.

CONCLUSION

MRF data from many sessions and scanners can potentially be treated as a single dataset for harmonized analysis or longitudinal comparisons without the additional regularization steps needed for qualitative modalities.

Default mode network shows distinct emotional and contextual responses yet common effects of retrieval demands across tasks

The default mode network (DMN) lies towards the heteromodal end of the principal gradient of intrinsic connectivity, maximally separated from the sensory-motor cortex. It supports memory-based cognition, including the capacity to retrieve conceptual and evaluative information from sensory inputs, and to generate meaningful states internally; however, the functional organisation of DMN that can support these distinct modes of retrieval remains unclear. We used fMRI to examine whether activation within subsystems of DMN differed as a function of retrieval demands, or the type of association to be retrieved, or both. In a picture association task, participants retrieved semantic associations that were either contextual or emotional in nature. Participants were asked to avoid generating episodic associations. In the generate phase, these associations were retrieved from a novel picture, while in the switch phase, participants retrieved a new association for the same image. Semantic context and emotion trials were associated with dissociable DMN subnetworks, indicating that a key dimension of DMN organisation relates to the type of association being accessed. The frontotemporal and medial temporal DMN showed a preference for emotional and semantic contextual associations, respectively. Relative to the generate phase, the switch phase recruited clusters closer to the heteromodal apex of the principal gradient-a cortical hierarchy separating unimodal and heteromodal regions. There were no differences in this effect between association types. Instead, memory switching was associated with a distinct subnetwork associated with controlled internal cognition. These findings delineate distinct patterns of DMN recruitment for different kinds of associations yet common responses across tasks that reflect retrieval demands.

Activity in the pontine reticular nuclei scales with handgrip force in humans

The neural pathways that contribute to force production in humans are currently poorly understood, as the relative roles of the corticospinal tract and brainstem pathways, such as the reticulospinal tract (RST), vary substantially across species. Using functional magnetic resonance imaging (fMRI), we aimed to measure activation in the pontine reticular nuclei (PRN) during different submaximal handgrip contractions to determine the potential role of the PRN in force modulation. Thirteen neurologically intact participants (age: 28 ± 6 yr) performed unilateral handgrip contractions at 25%, 50%, 75% of maximum voluntary contraction during brain scans. We quantified the magnitude of PRN activation from the contralateral and ipsilateral sides during each of the three contraction intensities. A repeated-measures ANOVA demonstrated a significant main effect of force (P = 0.012, [Formula: see text] = 0.307) for PRN activation, independent of side (i.e., activation increased with force for both contralateral and ipsilateral nuclei). Further analyses of these data involved calculating the linear slope between the magnitude of activation and handgrip force for each region of interest (ROI) at the individual-level. One-sample t tests on the slopes revealed significant group-level scaling for the PRN bilaterally, but only the ipsilateral PRN remained significant after correcting for multiple comparisons. We show evidence of task-dependent activation in the PRN that was positively related to handgrip force. These data build on a growing body of literature that highlights the RST as a functionally relevant motor pathway for force modulation in humans.NEW & NOTEWORTHY In this study, we used a task-based functional magnetic resonance imaging (fMRI) paradigm to show that activity in the pontine reticular nuclei scales linearly with increasing force during a handgrip task. These findings directly support recently proposed hypotheses that the reticulospinal tract may play an important role in modulating force production in humans.

Extensive T1-weighted MRI preprocessing improves generalizability of deep brain age prediction models

Brain age is an estimate of chronological age obtained from T1-weighted magnetic resonance images (T1w MRI), representing a straightforward diagnostic biomarker of brain aging and associated diseases. While the current best accuracy of brain age predictions on T1w MRIs of healthy subjects ranges from two to three years, comparing results across studies is challenging due to differences in the datasets, T1w preprocessing pipelines, and evaluation protocols used. This paper investigates the impact of T1w image preprocessing on the performance of four deep learning brain age models from recent literature. Four preprocessing pipelines, which differed in terms of registration transform, grayscale correction, and software implementation, were evaluated. The results showed that the choice of software or preprocessing steps could significantly affect the prediction error, with a maximum increase of 0.75 years in mean absolute error (MAE) for the same model and dataset. While grayscale correction had no significant impact on MAE, using affine rather than rigid registration to brain atlas statistically significantly improved MAE. Models trained on 3D images with isotropic 1mm3 resolution exhibited less sensitivity to the T1w preprocessing variations compared to 2D models or those trained on downsampled 3D images. Our findings indicate that extensive T1w preprocessing improves MAE, especially when predicting on a new dataset. This runs counter to prevailing research literature, which suggests that models trained on minimally preprocessed T1w scans are better suited for age predictions on MRIs from unseen scanners. We demonstrate that, irrespective of the model or T1w preprocessing used during training, applying some form of offset correction is essential to enable the model's performance to generalize effectively on datasets from unseen sites, regardless of whether they have undergone the same or different T1w preprocessing as the training set.

Multi-Template Meta-Information Regularized Network for Alzheimer's Disease Diagnosis Using Structural MRI

Structural magnetic resonance imaging (sMRI) has been widely applied in computer-aided Alzheimer's disease (AD) diagnosis, owing to its capabilities in providing detailed brain morphometric patterns and anatomical features in vivo. Although previous works have validated the effectiveness of incorporating metadata (e.g., age, gender, and educational years) for sMRI-based AD diagnosis, existing methods solely paid attention to metadata-associated correlation to AD (e.g., gender bias in AD prevalence) or confounding effects (e.g., the issue of normal aging and metadata-related heterogeneity). Hence, it is difficult to fully excavate the influence of metadata on AD diagnosis. To address these issues, we constructed a novel Multi-template Meta-information Regularized Network (MMRN) for AD diagnosis. Specifically, considering diagnostic variation resulting from different spatial transformations onto different brain templates, we first regarded different transformations as data augmentation for self-supervised learning after template selection. Since the confounding effects may arise from excessive attention to meta-information owing to its correlation with AD, we then designed the modules of weakly supervised meta-information learning and mutual information minimization to learn and disentangle meta-information from learned class-related representations, which accounts for meta-information regularization for disease diagnosis. We have evaluated our proposed MMRN on two public multi-center cohorts, including the Alzheimer's Disease Neuroimaging Initiative (ADNI) with 1,950 subjects and the National Alzheimer's Coordinating Center (NACC) with 1,163 subjects. The experimental results have shown that our proposed method outperformed the state-of-the-art approaches in both tasks of AD diagnosis, mild cognitive impairment (MCI) conversion prediction, and normal control (NC) vs. MCI vs. AD classification.

A biophysically constrained brain connectivity model based on stimulation-evoked potentials

BACKGROUND

Single-pulse electrical stimulation (SPES) is an established technique used to map functional effective connectivity networks in treatment-refractory epilepsy patients undergoing intracranial-electroencephalography monitoring. While the connectivity path between stimulation and recording sites has been explored through the integration of structural connectivity, there are substantial gaps, such that new modeling approaches may advance our understanding of connectivity derived from SPES studies.

NEW METHOD

Using intracranial electrophysiology data recorded from a single patient undergoing stereo-electroencephalography (sEEG) evaluation, we employ an automated detection method to identify early response components, C1, from pulse-evoked potentials (PEPs) induced by SPES. C1 components were utilized for a novel topology optimization method, modeling 3D electrical conductivity to infer neural pathways from stimulation sites. Additionally, PEP features were compared with tractography metrics, and model results were analyzed with respect to anatomical features.

RESULTS

The proposed optimization model resolved conductivity paths with low error. Specific electrode contacts displaying high error correlated with anatomical complexities. The C1 component strongly correlated with additional PEP features and displayed stable, weak correlations with tractography measures.

COMPARISON WITH EXISTING METHOD

Existing methods for estimating neural signal pathways are imaging-based and thus rely on anatomical inferences.

CONCLUSIONS

These results demonstrate that informing topology optimization methods with human intracranial SPES data is a feasible method for generating 3D conductivity maps linking electrical pathways with functional neural ensembles. PEP-estimated effective connectivity is correlated with but distinguished from structural connectivity. Modeled conductivity resolves connectivity pathways in the absence of anatomical priors.

Three-dimensional architecture and moment arms of human rotator cuff muscles in vivo: Interindividual, intermuscular, and intramuscular variations

The human rotator cuff consists of four muscles, each with a complex, multipennate architecture. Despite the functional and clinical importance, the architecture of the human rotator cuff has yet to be clearly described in humans in vivo. The purpose of this study was to investigate the intramuscular, intermuscular, and interindividual variations in architecture and moment arms of the human rotator cuff. Muscle volumes, fascicle lengths, physiological cross-sectional areas (PCSAs), pennation angles, and moment arms of all four rotator cuff muscles were measured from mDixon and diffusion tensor imaging (DTI) scans of the right shoulders of 20 young adults. In accordance with the most detailed dissections available to date, we found substantial intramuscular variation in fascicle length (coefficients of variation (CVs) ranged from 26% to 40%) and pennation angles (CVs ranged from 56% to 62%) in all rotator cuff muscles. We also found substantial intermuscular and interindividual variations in muscle volumes, but relatively consistent mean fascicle lengths, pennation angles, and moment arms (CVs for all ≤17%). Moreover, when expressed as a proportion of total rotator cuff muscle volume, the volumes of individual rotator cuff muscles were highly consistent between individuals and sexes (CVs ≤16%), suggesting that rotator cuff muscle volumes scale uniformly, at least in a younger population without musculoskeletal problems. Together, these data indicate limited interindividual and intermuscular variability in architecture, which may simplify scaling routines for musculoskeletal models. However, the substantial intramuscular variation in architecture questions the validity of previously reported mean architectural parameters to adequately describe rotator cuff function.

brainlife.io: a decentralized and open-source cloud platform to support neuroscience research

Neuroscience is advancing standardization and tool development to support rigor and transparency. Consequently, data pipeline complexity has increased, hindering FAIR (findable, accessible, interoperable and reusable) access. brainlife.io was developed to democratize neuroimaging research. The platform provides data standardization, management, visualization and processing and automatically tracks the provenance history of thousands of data objects. Here, brainlife.io is described and evaluated for validity, reliability, reproducibility, replicability and scientific utility using four data modalities and 3,200 participants.

One model to unite them all: Personalized federated learning of multi-contrast MRI synthesis

Curation of large, diverse MRI datasets via multi-institutional collaborations can help improve learning of generalizable synthesis models that reliably translate source- onto target-contrast images. To facilitate collaborations, federated learning (FL) adopts decentralized model training while mitigating privacy concerns by avoiding sharing of imaging data. However, conventional FL methods can be impaired by the inherent heterogeneity in the data distribution, with domain shifts evident within and across imaging sites. Here we introduce the first personalized FL method for MRI Synthesis (pFLSynth) that improves reliability against data heterogeneity via model specialization to individual sites and synthesis tasks (i.e., source-target contrasts). To do this, pFLSynth leverages an adversarial model equipped with novel personalization blocks that control the statistics of generated feature maps across the spatial/channel dimensions, given latent variables specific to sites and tasks. To further promote communication efficiency and site specialization, partial network aggregation is employed over later generator stages while earlier generator stages and the discriminator are trained locally. As such, pFLSynth enables multi-task training of multi-site synthesis models with high generalization performance across sites and tasks. Comprehensive experiments demonstrate the superior performance and reliability of pFLSynth in MRI synthesis against prior federated methods.

In primary visual cortex fMRI responses to chromatic and achromatic stimuli are interdependent and predict contrast detection thresholds

Chromatic and achromatic signals in primary visual cortex have historically been considered independent of each other but have since shown evidence of interdependence. Here, we investigated the combination of two components of a stimulus; an achromatic dynamically changing check background and a chromatic (L-M or S cone) target grating. We found that combinations of chromatic and achromatic signals in primary visual cortex were interdependent, with the dynamic range of responses to chromatic contrast decreasing as achromatic contrast increased. A contrast detection threshold study also revealed interdependence of background and target, with increasing chromatic contrast detection thresholds as achromatic background contrast increased. A model that incorporated a normalising effect of achromatic contrast on chromatic responses, but not vice versa, best predicted our V1 data as well as behavioural thresholds. Further along the visual hierarchy, the dynamic range of chromatic responses was maintained when compared to achromatic responses, which became increasingly compressive.

Re-awakening the brain: Forcing transitions in disorders of consciousness by external in silico perturbation

A fundamental challenge in neuroscience is accurately defining brain states and predicting how and where to perturb the brain to force a transition. Here, we investigated resting-state fMRI data of patients suffering from disorders of consciousness (DoC) after coma (minimally conscious and unresponsive wakefulness states) and healthy controls. We applied model-free and model-based approaches to help elucidate the underlying brain mechanisms of patients with DoC. The model-free approach allowed us to characterize brain states in DoC and healthy controls as a probabilistic metastable substate (PMS) space. The PMS of each group was defined by a repertoire of unique patterns (i.e., metastable substates) with different probabilities of occurrence. In the model-based approach, we adjusted the PMS of each DoC group to a causal whole-brain model. This allowed us to explore optimal strategies for promoting transitions by applying off-line in silico probing. Furthermore, this approach enabled us to evaluate the impact of local perturbations in terms of their global effects and sensitivity to stimulation, which is a model-based biomarker providing a deeper understanding of the mechanisms underlying DoC. Our results show that transitions were obtained in a synchronous protocol, in which the somatomotor network, thalamus, precuneus and insula were the most sensitive areas to perturbation. This motivates further work to continue understanding brain function and treatments of disorders of consciousness.

A roadmap for implanting microelectrode arrays to evoke tactile sensations through intracortical microstimulation

Intracortical microstimulation (ICMS) is a method for restoring sensation to people with paralysis as part of a bidirectional brain-computer interface to restore upper limb function. Evoking tactile sensations of the hand through ICMS requires precise targeting of implanted electrodes. Here we describe the presurgical imaging procedures used to generate functional maps of the hand area of the somatosensory cortex and subsequent planning that guided the implantation of intracortical microelectrode arrays. In five participants with cervical spinal cord injury, across two study locations, this procedure successfully enabled ICMS-evoked sensations localized to at least the first four digits of the hand. The imaging and planning procedures developed through this clinical trial provide a roadmap for other brain-computer interface studies to ensure successful placement of stimulation electrodes.

Accelerated epigenetic age is associated with whole-brain functional connectivity and impaired cognitive performance in older adults

While chronological age is a strong predictor for health-related risk factors, it is an incomplete metric that fails to fully characterize the unique aging process of individuals with different genetic makeup, neurodevelopment, and environmental experiences. Recent advances in epigenomic array technologies have made it possible to generate DNA methylation-based biomarkers of biological aging, which may be useful in predicting a myriad of cognitive abilities and functional brain network organization across older individuals. It is currently unclear which cognitive domains are negatively correlated with epigenetic age above and beyond chronological age, and it is unknown if functional brain organization is an important mechanism for explaining these associations. In this study, individuals with accelerated epigenetic age (i.e. AgeAccelGrim) performed worse on tasks that spanned a wide variety of cognitive faculties including both fluid and crystallized intelligence (N = 103, average age = 68.98 years, 73 females, 30 males). Additionally, fMRI connectome-based predictive models suggested a mediating mechanism of functional connectivity on epigenetic age acceleration-cognition associations primarily in medial temporal lobe and limbic structures. This research highlights the important role of epigenetic aging processes on the development and maintenance of healthy cognitive capacities and function of the aging brain.

Naturalistic drug cue reactivity in heroin use disorder: orbitofrontal synchronization as a marker of craving and recovery

Movies captivate groups of individuals (the audience), especially if they contain themes of common motivational interest to the group. In drug addiction, a key mechanism is maladaptive motivational salience attribution whereby drug cues outcompete other reinforcers within the same environment or context. We predicted that while watching a drug-themed movie, where cues for drugs and other stimuli share a continuous narrative context, fMRI responses in individuals with heroin use disorder (iHUD) will preferentially synchronize during drug scenes. Results revealed such drug-biased synchronization in the orbitofrontal cortex (OFC), ventromedial and ventrolateral prefrontal cortex, and insula. After 15 weeks of inpatient treatment, there was a significant reduction in this drug-biased shared response in the OFC, which correlated with a concomitant reduction in dynamically-measured craving, suggesting synchronized OFC responses to a drug-themed movie as a neural marker of craving and recovery in iHUD.

Greater Pattern Similarity between Mother Tongue and Second Language in the Right ATL Facilitates Understanding of Written Language

Previous research has mapped out the brain regions that respond to semantic stimuli presented visually and auditorily, but there is debate about whether semantic representation is modality-specific (only written or only spoken) or modality-invariant (both written and spoken). The mechanism of semantic representation underlying native (L1) and second language (L2) comprehension in different modalities as well as how this mechanism is influenced by L2 proficiency, remains unclear. We used functional magnetic resonance imaging (fMRI) data from the OpenNEURO database to calculate neural pattern similarity across native and second languages (Spanish and English) for different input modalities (written and spoken) and learning sessions (before and after training). The correlations between behavioral performance and cross-language pattern similarity for L1 and L2 were also calculated. Spanish-English bilingual adolescents (N = 24; ages 16-17; 19 girls) participated in a 3-month English immersion after-school program. As L2 proficiency increased, greater cross-language pattern similarity between L1 and L2 spoken words was observed in the left pars triangularis. Cross-language pattern similarity between L1 and L2 written words was observed in the right anterior temporal lobe. Brain-behavior correlations indicated that increased cross-language pattern similarity between L1 and L2 written words in the right anterior temporal lobe was associated with L2 written word comprehension. This study identified an effective neurofunctional predictor related to L2 written word comprehension.

Longitudinal microstructural changes in 18 amygdala nuclei resonate with cortical circuits and phenomics

The amygdala nuclei modulate distributed neural circuits that most likely evolved to respond to environmental threats and opportunities. So far, the specific role of unique amygdala nuclei in the context processing of salient environmental cues lacks adequate characterization across neural systems and over time. Here, we present amygdala nuclei morphometry and behavioral findings from longitudinal population data (>1400 subjects, age range 40-69 years, sampled 2-3 years apart): the UK Biobank offers exceptionally rich phenotyping along with brain morphology scans. This allows us to quantify how 18 microanatomical amygdala subregions undergo plastic changes in tandem with coupled neural systems and delineating their associated phenome-wide profiles. In the context of population change, the basal, lateral, accessory basal, and paralaminar nuclei change in lockstep with the prefrontal cortex, a region that subserves planning and decision-making. The central, medial and cortical nuclei are structurally coupled with the insular and anterior-cingulate nodes of the salience network, in addition to the MT/V5, basal ganglia, and putamen, areas proposed to represent internal bodily states and mediate attention to environmental cues. The central nucleus and anterior amygdaloid area are longitudinally tied with the inferior parietal lobule, known for a role in bodily awareness and social attention. These population-level amygdala-brain plasticity regimes in turn are linked with unique collections of phenotypes, ranging from social status and employment to sleep habits and risk taking. The obtained structural plasticity findings motivate hypotheses about the specific functions of distinct amygdala nuclei in humans.

Goal commitment is supported by vmPFC through selective attention

When striking a balance between commitment to a goal and flexibility in the face of better options, people often demonstrate strong goal perseveration. Here, using functional MRI (n = 30) and lesion patient (n = 26) studies, we argue that the ventromedial prefrontal cortex (vmPFC) drives goal commitment linked to changes in goal-directed selective attention. Participants performed an incremental goal pursuit task involving sequential decisions between persisting with a goal versus abandoning progress for better alternative options. Individuals with stronger goal perseveration showed higher goal-directed attention in an interleaved attention task. Increasing goal-directed attention also affected abandonment decisions: while pursuing a goal, people lost their sensitivity to valuable alternative goals while remaining more sensitive to changes in the current goal. In a healthy population, individual differences in both commitment biases and goal-oriented attention were predicted by baseline goal-related activity in the vmPFC. Among lesion patients, vmPFC damage reduced goal commitment, leading to a performance benefit.

Non-invasive quantification of 18F-florbetaben with total-body EXPLORER PET

BACKGROUND

Kinetic modeling of 18F-florbetaben provides important quantification of brain amyloid deposition in research and clinical settings but its use is limited by the requirement of arterial blood data for quantitative PET. The total-body EXPLORER PET scanner supports the dynamic acquisition of a full human body simultaneously and permits noninvasive image-derived input functions (IDIFs) as an alternative to arterial blood sampling. This study quantified brain amyloid burden with kinetic modeling, leveraging dynamic 18F-florbetaben PET in aorta IDIFs and the brain in an elderly cohort.

METHODS

18F-florbetaben dynamic PET imaging was performed on the EXPLORER system with tracer injection (300 MBq) in 3 individuals with Alzheimer's disease (AD), 3 with mild cognitive impairment, and 9 healthy controls. Image-derived input functions were extracted from the descending aorta with manual regions of interest based on the first 30 s after injection. Dynamic time-activity curves (TACs) for 110 min were fitted to the two-tissue compartment model (2TCM) using population-based metabolite corrected IDIFs to calculate total and specific distribution volumes (VT, Vs) in key brain regions with early amyloid accumulation. Non-displaceable binding potential ([Formula: see text] was also calculated from the multi-reference tissue model (MRTM).

RESULTS

Amyloid-positive (AD) patients showed the highest VT and VS in anterior cingulate, posterior cingulate, and precuneus, consistent with [Formula: see text] analysis. [Formula: see text]and VT from kinetic models were correlated (r² = 0.46, P < 2[Formula: see text] with a stronger positive correlation observed in amyloid-positive participants, indicating reliable model fits with the IDIFs. VT from 2TCM was highly correlated ([Formula: see text]= 0.65, P < 2[Formula: see text]) with Logan graphical VT estimation.

CONCLUSION

Non-invasive quantification of amyloid binding from total-body 18F-florbetaben PET data is feasible using aorta IDIFs with high agreement between kinetic distribution volume parameters compared to [Formula: see text]in amyloid-positive and amyloid-negative older individuals.

Comprehensive analysis of synthetic learning applied to neonatal brain MRI segmentation

Brain segmentation from neonatal MRI images is a very challenging task due to large changes in the shape of cerebral structures and variations in signal intensities reflecting the gestational process. In this context, there is a clear need for segmentation techniques that are robust to variations in image contrast and to the spatial configuration of anatomical structures. In this work, we evaluate the potential of synthetic learning, a contrast-independent model trained using synthetic images generated from the ground truth labels of very few subjects. We base our experiments on the dataset released by the developmental Human Connectome Project, for which high-quality images are available for more than 700 babies aged between 26 and 45 weeks postconception. First, we confirm the impressive performance of a standard UNet trained on a few volumes, but also confirm that such models learn intensity-related features specific to the training domain. We then confirm the robustness of the synthetic learning approach to variations in image contrast. However, we observe a clear influence of the age of the baby on the predictions. We improve the performance of this model by enriching the synthetic training set with realistic motion artifacts and over-segmentation of the white matter. Based on extensive visual assessment, we argue that the better performance of the model trained on real T2w data may be due to systematic errors in the ground truth. We propose an original experiment allowing us to show that learning from real data will reproduce any systematic bias affecting the training set, while synthetic models can avoid this limitation. Overall, our experiments confirm that synthetic learning is an effective solution for segmenting neonatal brain MRI. Our adapted synthetic learning approach combines key features that will be instrumental for large multisite studies and clinical applications.

Resting heart rate variability is associated with neural adaptation when repeatedly exposed to emotional stimuli

Higher heart rate variability (HRV) at rest is associated with better emotion regulation ability. While the neurovisceral integration model explains this by postulating that HRV can index how the brain adaptively modulates responses to emotional stimuli, neuroimaging studies directly supporting this idea are scarce. We examined the neural correlates of regulating negative and positive emotion in relation to resting HRV based on the neuroimaging and heart rate data of one hundred young adults. The results showed that those with higher HRV better recruit the medial prefrontal cortex while intensifying positive compared to negative emotion. We also examined how individual differences in resting HRV are associated with adjusting brain activity to repeated emotional stimuli. During repeated viewing of emotional images, subjects with higher resting HRV better reduced activity in the medial prefrontal cortex, posterior cingulate gyrus, and angular gyrus, most of which overlapped with the default mode network. This HRV-DMN association was observed during passively viewing emotional images rather than during actively regulating emotion. While the regulating trials can better detect task-induced changes, the viewing trials might approximate resting state, better revealing individual differences. These findings suggest two possibilities: people with higher resting HRV might have a tendency to spontaneously engage with emotion regulation or possess a trait helping emotional arousal fade away.

Role of corpus callosum in unconscious vision

The existence of unconscious visually triggered behavior in patients with cortical blindness (e.g., homonymous hemianopia) has been amply demonstrated and the neural bases of this phenomenon have been thoroughly studied. However, a crosstalk between the two hemispheres as a possible mechanism of unconscious or partially conscious vision has not been so far considered. Thus, the aim of this study was to assess the relationship between structural and functional properties of the corpus callosum (CC), as shown by probabilistic tractography (PT), behavioral detection/discrimination performance and level of perceptual awareness in the blind field of patients with hemianopia. Twelve patients were tested in two tasks with black-and-white visual square-wave gratings, one task of movement and the other of orientation. The stimuli were lateralized to one hemifield either intact or blind. A PT analysis was carried out on MRI data to extract fiber properties along the CC (genu, body, and splenium). Compared with a control group of participants without brain damage, patients showed lower FA values in all three CC sections studied. For the intact hemifield we found a significant correlation between PT values and visual detection/discrimination accuracy. For the blind hemifield the level of perceptual awareness correlated with PT values for all three CC sections in the movement task. Importantly, significant differences in all three CC sections were found also between patients with above-vs. chance detection/discrimination performance while differences in the genu were found between patients with and without perceptual awareness. Overall, our study provides evidence that the properties of CC fibers are related to the presence of unconscious stimulus detection/discrimination and to hints of perceptual awareness for stimulus presentation to the blind hemifield. These results underline the importance of information exchange between the damaged and the healthy hemisphere for possible partial or full recovery from hemianopia.

Differing patterns of cortical grey matter pathology identified by multifractal analysis in UMN-predominant ALS patients with and without corticospinal tract hyperintensity

The pathological hallmarks of amyotrophic lateral sclerosis (ALS) are degeneration of the primary motor cortex grey matter (GM) and corticospinal tract (CST) resulting in upper motor neuron (UMN) dysfunction. Conventional brain magnetic resonance imaging (MRI) shows abnormal CST hyperintensity in some UMN-predominant ALS patients (ALS-CST+) but not in others (ALS-CST-). In addition to the CST differences, we aimed to determine whether GM degeneration differs between ALS-CST+ and ALS-CST- patients by cortical thickness (CT), voxel-based morphometry (VBM) and fractal dimension analyses. We hypothesized that MRI multifractal (MF) measures could differentiate between neurologic controls (n = 14) and UMN-predominant ALS patients as well as between patient subgroups (ALS-CST+, n = 21 vs ALS-CST-, n = 27). No significant differences were observed in CT or GM VBM in any brain regions between patients and controls or between ALS subgroups. MF analyses were performed separately on GM of the whole brain, of frontal, parietal, occipital, and temporal lobes as well as of cerebellum. Estimating MF measures D (Q = 0), D (Q = 1), D (Q = 2), Δf, Δα of frontal lobe GM classified neurologic controls, ALS-CST+ and ALS-CST- groups with 98% accuracy and > 95% in F1, recall, precision and specificity scores. Classification accuracy was only 74% when using whole brain MF measures and < 70% for other brain lobes. We demonstrate that MF analysis can distinguish UMN-predominant ALS subgroups based on GM changes, which the more commonly used quantitative approaches of CT and VBM cannot.

Water Diffusion in the Live Human Brain is Gaussian at the Mesoscale

Imaging the live human brain at the mesoscopic scale is a desideratum in basic and clinical neurosciences. Despite the promise of diffusion MRI, the lack of an accurate model relating the measured signal and the associated microstructure has hampered its success. The widely used diffusion tensor MRI (DTI) model assumes an anisotropic Gaussian diffusion process in each voxel, but lacks the ability to capture intravoxel heterogeneity. This study explores the extension of the DTI model to mesoscopic length scales by use of the diffusion tensor distribution (DTD) model, which assumes a Gaussian diffusion process in each subvoxel. DTD MRI has shown promise in addressing some limitations of DTI, particularly in distinguishing among different types of brain cancers and elucidating multiple fiber populations within a voxel. However, its validity in live brain tissue has never been established. Here, multiple diffusion-encoded (MDE) data were acquired in the living human brain using a 3 Tesla MRI scanner with large diffusion weighting factors. Two different diffusion times (Δ = 37, 74 ms) were employed, with other scanning parameters fixed to assess signal decay differences. In vivo diffusion-weighted signals in gray and white matter were nearly identical at the two diffusion times. Fitting the signals to the DTD model yielded indistinguishable results, except in the cerebrospinal fluid (CSF)-filled voxels likely due to pulsatile flow. Overall, the study supports the time invariance of water diffusion at the mesoscopic scale in live brain parenchyma, extending the validity of the anisotropic Gaussian diffusion model in clinical brain imaging.

Thalamic feedback shapes brain responses evoked by cortical stimulation in mice and humans

Cortical stimulation with single pulses is a common technique in clinical practice and research. However, we still do not understand the extent to which it engages subcortical circuits which contribute to the associated evoked potentials (EPs). Here we find that cortical stimulation generates remarkably similar EPs in humans and mice, with a late component similarly modulated by the subject's behavioral state. We optogenetically dissect the underlying circuit in mice, demonstrating that the late component of these EPs is caused by a thalamic hyperpolarization and rebound. The magnitude of this late component correlates with the bursting frequency and synchronicity of thalamic neurons, modulated by the subject's behavioral state. A simulation of the thalamo-cortical circuit highlights that both intrinsic thalamic currents as well as cortical and thalamic GABAergic neurons contribute to this response profile. We conclude that the cortical stimulation engages cortico-thalamo-cortical circuits highly preserved across different species and stimulation modalities.

Early life adversity is associated with greater similarity in neural representations of ambiguous and threatening stimuli

Exposure to early life adversity (ELA) is hypothesized to sensitize threat-responsive neural circuitry. This may lead individuals to overestimate threat in the face of ambiguity, a cognitive-behavioral phenotype linked to poor mental health. The tendency to process ambiguity as threatening may stem from difficulty distinguishing between ambiguous and threatening stimuli. However, it is unknown how exposure to ELA relates to neural representations of ambiguous and threatening stimuli, or how processing of ambiguity following ELA relates to psychosocial functioning. The current fMRI study examined multivariate representations of threatening and ambiguous social cues in 41 emerging adults (aged 18 to 19 years). Using representational similarity analysis, we assessed neural representations of ambiguous and threatening images within affective neural circuitry and tested whether similarity in these representations varied by ELA exposure. Greater exposure to ELA was associated with greater similarity in neural representations of ambiguous and threatening images. Moreover, individual differences in processing ambiguity related to global functioning, an association that varied as a function of ELA. By evidencing reduced neural differentiation between ambiguous and threatening cues in ELA-exposed emerging adults and linking behavioral responses to ambiguity to psychosocial wellbeing, these findings have important implications for future intervention work in at-risk, ELA-exposed populations.

Can the neural representation of physical pain predict empathy for pain in others?

The question of whether physical pain and vicarious pain have some shared neural substrates is unresolved. Recent research has argued that physical and vicarious pain are represented by dissociable multivariate brain patterns by creating biomarkers for physical pain (Neurologic Pain Signature, NPS) and vicarious pain (Vicarious Pain Signature, VPS), respectively. In the current research, the NPS and two versions of the VPS were applied to three fMRI datasets (one new, two published) relating to vicarious pain which focused on between-subject differences in vicarious pain (Datasets 1 and 3) and within-subject manipulations of perspective taking (Dataset 2). Results show that (i) NPS can distinguish brain responses to images of pain vs no-pain and to a greater extent in vicarious pain responders who report experiencing pain when observing pain and (ii) neither version of the VPS mapped on to individual differences in vicarious pain and the two versions differed in their success in predicting vicarious pain overall. This study suggests that the NPS (created to detect physical pain) is, under some circumstances, sensitive to vicarious pain and there is significant variability in VPS measures (created to detect vicarious pain) to act as generalizable biomarkers of vicarious pain.

Commonality and variance of resting-state networks in common marmoset brains

Animal models of brain function are critical for the study of human diseases and development of effective interventions. Resting-state network (RSN) analysis is a powerful tool for evaluating brain function and performing comparisons across animal species. Several studies have reported RSNs in the common marmoset (Callithrix jacchus; marmoset), a non-human primate. However, it is necessary to identify RSNs and evaluate commonality and inter-individual variance through analyses using a larger amount of data. In this study, we present marmoset RSNs detected using > 100,000 time-course image volumes of resting-state functional magnetic resonance imaging data with careful preprocessing. In addition, we extracted brain regions involved in the composition of these RSNs to understand the differences between humans and marmosets. We detected 16 RSNs in major marmosets, three of which were novel networks that have not been previously reported in marmosets. Since these RSNs possess the potential for use in the functional evaluation of neurodegenerative diseases, the data in this study will significantly contribute to the understanding of the functional effects of neurodegenerative diseases.

A subject-specific unsupervised deep learning method for quantitative susceptibility mapping using implicit neural representation

Quantitative susceptibility mapping (QSM) is an MRI-based technique that estimates the underlying tissue magnetic susceptibility based on phase signal. Deep learning (DL)-based methods have shown promise in handling the challenging ill-posed inverse problem for QSM reconstruction. However, they require extensive paired training data that are typically unavailable and suffer from generalization problems. Recent model-incorporated DL approaches also overlook the non-local effect of the tissue phase in applying the source-to-field forward model due to patch-based training constraint, resulting in a discrepancy between the prediction and measurement and subsequently suboptimal QSM reconstruction. This study proposes an unsupervised and subject-specific DL method for QSM reconstruction based on implicit neural representation (INR), referred to as INR-QSM. INR has emerged as a powerful framework for learning a high-quality continuous representation of the signal (image) by exploiting its internal information without training labels. In INR-QSM, the desired susceptibility map is represented as a continuous function of the spatial coordinates, parameterized by a fully-connected neural network. The weights are learned by minimizing a loss function that includes a data fidelity term incorporated by the physical model and regularization terms. Additionally, a novel phase compensation strategy is proposed for the first time to account for the non-local effect of tissue phase in data consistency calculation to make the physical model more accurate. Our experiments show that INR-QSM outperforms traditional established QSM reconstruction methods and the compared unsupervised DL method both qualitatively and quantitatively, and is competitive against supervised DL methods under data perturbations.

Early Infant Prefrontal Cortical Microstructure Predicts Present and Future Emotionality

BACKGROUND

High levels of infant negative emotionality (NE) and low positive emotionality (PE) predict future emotional and behavioral problems. The prefrontal cortex (PFC) supports emotional regulation, with each PFC subregion specializing in specific emotional processes. Neurite orientation dispersion and density imaging estimates microstructural integrity and myelination via the neurite density index (NDI) and dispersion via the orientation dispersion index (ODI), with potential to more accurately evaluate microstructural alterations in the developing brain. Yet, no study has used these indices to examine associations between PFC microstructure and concurrent or developing infant emotionality.

METHODS

We modeled PFC subregional NDI and ODI at 3 months with caregiver-reported infant NE and PE at 3 months (n = 61) and at 9 months (n = 50), using multivariable and subsequent bivariate regression models.

RESULTS

The most robust statistically significant findings were positive associations among 3-month rostral anterior cingulate cortex (ACC) ODI and caudal ACC NDI and concurrent NE, a positive association between 3-month lateral orbitofrontal cortex ODI and prospective NE, and a negative association between 3-month dorsolateral PFC ODI and concurrent PE. Multivariate models also revealed that other PFC subregional microstructure measures, as well as infant and caregiver sociodemographic and clinical factors, predicted infant 3- and 9-month NE and PE.

CONCLUSIONS

Greater NDI and ODI, reflecting greater microstructural complexity, in PFC regions supporting salience perception (rostral ACC), decision making (lateral orbitofrontal cortex), action selection (caudal ACC), and attentional processes (dorsolateral PFC) might result in greater integration of these subregions with other neural networks and greater attention to salient negative external cues, thus higher NE and/or lower PE. These findings provide potential infant cortical markers of future psychopathology risk.

Quantitative Assessment of Preanalytic Variables on Clinical Evaluation of PI3/AKT/mTOR Signaling Activity in Diffuse Glioma

Biomarker-driven therapeutic clinical trials require the implementation of standardized, evidence-based practices for sample collection. In diffuse glioma, phosphatidylinositol 3 (PI3)-kinase/AKT/mTOR (PI3/AKT/mTOR) signaling is an attractive therapeutic target for which window-of-opportunity clinical trials could facilitate the identification of promising new agents. Yet, the relevant preanalytic variables and optimal tumor sampling methods necessary to measure pathway activity are unknown. To address this, we used a murine model for isocitrate dehydrogenase (IDH)-wildtype glioblastoma (GBM) and human tumor tissue, including IDH-wildtype GBM and IDH-mutant diffuse glioma. First, we determined the impact of delayed time-to-formalin fixation, or cold ischemia time (CIT), on the quantitative assessment of cellular expression of 6 phosphoproteins that are readouts of PI3K/AK/mTOR activity (phosphorylated-proline-rich Akt substrate of 40 kDa (p-PRAS40, T246), -mechanistic target of rapamycin (p-mTOR; S2448); -AKT (p-AKT, S473); -ribosomal protein S6 (p-RPS6, S240/244 and S235/236), and -eukaryotic initiation factor 4E-binding protein 1 (p-4EBP1, T37/46). With CITs ≥ 2 hours, typical of routine clinical handling, all had reduced or altered expression with p-RPS6 (S240/244) exhibiting relatively greater stability. A similar pattern was observed using patient tumor samples from the operating room with p-4EBP1 more sensitive to delayed fixation than p-RPS6 (S240/244). Many clinical trials utilize unstained slides for biomarker evaluation. Thus, we evaluated the impact of slide storage conditions on the detection of p-RPS6 (S240/244), p-4EBP1, and p-AKT. After 5 months, storage at -80°C was required to preserve the expression of p-4EBP1 and p-AKT, whereas p-RPS6 (240/244) expression was not stable regardless of storage temperature. Biomarker heterogeneity impacts optimal tumor sampling. Quantification of p-RPS6 (240/244) expression in multiple regionally distinct human tumor samples from 8 patients revealed significant intratumoral heterogeneity. Thus, the accurate assessment of PI3K/AKT/mTOR signaling in diffuse glioma must overcome intratumoral heterogeneity and multiple preanalytic factors, including time-to-formalin fixation, slide storage conditions, and phosphoprotein of interest.

Conversion map from quantitative parameter mapping to myelin water fraction: comparison with R1·R2* and myelin water fraction in white matter

OBJECTIVE

To clarify the relationship between myelin water fraction (MWF) and R1⋅R2* and to develop a method to calculate MWF directly from parameters derived from QPM, i.e., MWF converted from QPM (MWFQPM).

MATERIALS AND METHODS

Subjects were 12 healthy volunteers. On a 3 T MR scanner, dataset was acquired using spoiled gradient-echo sequence for QPM. MWF and R1⋅R2* maps were derived from the multi-gradient-echo (mGRE) dataset. Volume-of-interest (VOI) analysis using the JHU-white matter (WM) atlas was performed. All the data in the 48 WM regions measured VOI were plotted, and quadratic polynomial approximations of each region were derived from the relationship between R1·R2* and the two-pool model-MWF. The R1·R2* map was converted to MWFQPM map. MWF atlas template was generated using converted to MWF from R1·R2* per WM region.

RESULTS

The mean MWF and R1·R2* values for the 48 WM regions were 11.96 ± 6.63%, and 19.94 ± 4.59 s-2, respectively. A non-linear relationship in 48 regions of the WM between MWF and R1·R2* values was observed by quadratic polynomial approximation (R2 ≥ 0.963, P < 0.0001).

DISCUSSION

MWFQPM map improved image quality compared to the mGRE-MWF map. Myelin water atlas template derived from MWFQPM may be generated with combined multiple WM regions.

Whole tumor analysis reveals early origin of the TERT promoter mutation and intercellular heterogeneity in TERT expression

BACKGROUND

The TERT promoter mutation (TPM) is acquired in most IDH-wildtype glioblastomas (GBM) and IDH-mutant oligodendrogliomas (OD) enabling tumor cell immortality. Previous studies on TPM clonality show conflicting results. This study was performed to determine whether TPM is clonal on a tumor-wide scale.

METHODS

We investigated TPM clonality in relation to presumed early events in 19 IDH-wildtype GBM and 10 IDH-mutant OD using 3-dimensional comprehensive tumor sampling. We performed Sanger sequencing on 264 tumor samples and deep amplicon sequencing on 187 tumor samples. We obtained tumor purity and copy number estimates from whole exome sequencing. TERT expression was assessed by RNA-seq and RNAscope.

RESULTS

We detected TPM in 100% of tumor samples with quantifiable tumor purity (219 samples). Variant allele frequencies (VAF) of TPM correlate positively with chromosome 10 loss in GBM (R = 0.85), IDH1 mutation in OD (R = 0.87), and with tumor purity (R = 0.91 for GBM; R = 0.90 for OD). In comparison, oncogene amplification was tumor-wide for MDM4- and most EGFR-amplified cases but heterogeneous for MYCN and PDGFRA, and strikingly high in low-purity samples. TPM VAF was moderately correlated with TERT expression (R = 0.52 for GBM; R = 0.65 for OD). TERT expression was detected in a subset of cells, solely in TPM-positive samples, including samples equivocal for tumor.

CONCLUSIONS

On a tumor-wide scale, TPM is among the earliest events in glioma evolution. Intercellular heterogeneity of TERT expression, however, suggests dynamic regulation during tumor growth. TERT expression may be a tumor cell-specific biomarker.

[1-11C]-Butanol Positron Emission Tomography reveals an impaired brain to nasal turbinates pathway in aging amyloid positive subjects

BACKGROUND

Reduced clearance of cerebrospinal fluid (CSF) has been suggested as a pathological feature of Alzheimer's disease (AD). With extensive documentation in non-human mammals and contradictory human neuroimaging data it remains unknown whether the nasal mucosa is a CSF drainage site in humans. Here, we used dynamic PET with [1-11C]-Butanol, a highly permeable radiotracer with no appreciable brain binding, to test the hypothesis that tracer drainage from the nasal pathway reflects CSF drainage from brain. As a test of the hypothesis, we examined whether brain and nasal fluid drainage times were correlated and affected by brain amyloid.

METHODS

24 cognitively normal subjects (≥ 65 years) were dynamically PET imaged for 60 min. using [1-11C]-Butanol. Imaging with either [11C]-PiB or [18F]-FBB identified 8 amyloid PET positive (Aβ+) and 16 Aβ- subjects. MRI-determined regions of interest (ROI) included: the carotid artery, the lateral orbitofrontal (LOF) brain, the cribriform plate, and an All-turbinate region comprised of the superior, middle, and inferior turbinates. The bilateral temporalis muscle and jugular veins served as control regions. Regional time-activity were used to model tracer influx, egress, and AUC.

RESULTS

LOF and All-turbinate 60 min AUC were positively associated, thus suggesting a connection between the brain and the nose. Further, the Aβ+ subgroup demonstrated impaired tracer kinetics, marked by reduced tracer influx and slower egress.

CONCLUSION

The data show that tracer kinetics for brain and nasal turbinates are related to each other and both reflect the amyloid status of the brain. As such, these data add to evidence that the nasal pathway is a potential CSF drainage site in humans. These data warrant further investigation of brain and nasal contributions to protein clearance in neurodegenerative disease.

Free water in gray matter linked to gut microbiota changes with decreased butyrate producers in Alzheimer's disease and mild cognitive impairment

Neuroinflammation contributes to the pathology and progression of Alzheimer's disease (AD), and it can be observed even with mild cognitive impairment (MCI), a prodromal phase of AD. Free water (FW) imaging estimates the extracellular water content and has been used to study neuroinflammation across several neurological diseases including AD. Recently, the role of gut microbiota has been implicated in the pathogenesis of AD. The relationship between FW imaging and gut microbiota was examined in patients with AD and MCI. Fifty-six participants underwent neuropsychological assessments, FW imaging, and gut microbiota analysis targeting the bacterial 16S rRNA gene. They were categorized into the cognitively normal control (NC) (n = 19), MCI (n = 19), and AD (n = 18) groups according to the neuropsychological assessments. The correlations of FW values, neuropsychological assessment scores, and the relative abundance of gut microbiota were analyzed. FW was higher in several white matter tracts and in gray matter regions, predominantly the frontal, temporal, limbic and paralimbic regions in the AD/MCI group than in the NC group. In the AD/MCI group, higher FW values in the temporal (superior temporal and temporal pole), limbic and paralimbic (insula, hippocampus and amygdala) regions were the most associated with worse neuropsychological assessment scores. In the AD/MCI group, FW values in these regions were negatively correlated with the relative abundances of butyrate-producing genera Anaerostipes, Lachnospiraceae UCG-004, and [Ruminococcus] gnavus group, which showed a significant decreasing trend in the order of the NC, MCI, and AD groups. The present study showed that increased FW in the gray matter regions related to cognitive impairment was associated with low abundances of butyrate producers in the AD/MCI group. These findings suggest an association between neuroinflammation and decreased levels of the short-chain fatty acid butyrate that is one of the major gut microbial metabolites having a potentially beneficial role in brain homeostasis.

Graded and sharp transitions in semantic function in left temporal lobe

Recent work has focussed on how patterns of functional change within the temporal lobe relate to whole-brain dimensions of intrinsic connectivity variation (Margulies et al., 2016). We examined two such 'connectivity gradients' reflecting the separation of (i) unimodal versus heteromodal and (ii) visual versus auditory-motor cortex, examining visually presented verbal associative and feature judgments, plus picture-based context and emotion generation. Functional responses along the first dimension sometimes showed graded change between modality-tuned and heteromodal cortex (in the verbal matching task), and other times showed sharp functional transitions, with deactivation at the extremes and activation in the middle of this gradient (internal generation). The second gradient revealed more visual than auditory-motor activation, regardless of content (associative, feature, context, emotion) or task process (matching/generation). We also uncovered subtle differences across each gradient for content type, which predominantly manifested as differences in relative magnitude of activation or deactivation.

Romer-EPTI: rotating-view motion-robust super-resolution EPTI for SNR-efficient distortion-free in-vivo mesoscale dMRI and microstructure imaging

Purpose

To overcome the major challenges in dMRI acquisition, including low SNR, distortion/blurring, and motion vulnerability.

Methods

A novel Romer-EPTI technique is developed to provide distortion-free dMRI with significant SNR gain, high motion-robustness, sharp spatial resolution, and simultaneous multi-TE imaging. It introduces a ROtating-view Motion-robust supEr-Resolution technique (Romer) combined with a distortion/blurring-free EPTI encoding. Romer enhances SNR by a simultaneous multi-thick-slice acquisition with rotating-view encoding, while providing high motion-robustness through a motion-aware super-resolution reconstruction, which also incorporates slice-profile and real-value diffusion, to resolve high-isotropic-resolution volumes. The in-plane encoding is performed using distortion/blurring-free EPTI, which further improves effective spatial resolution and motion robustness by preventing not only T2/T2*-blurring but also additional blurring resulting from combining encoded volumes with inconsistent geometries caused by dynamic distortions. Self-navigation was incorporated to enable efficient phase correction. Additional developments include strategies to address slab-boundary artifacts, achieve minimal TE for SNR gain at 7T, and achieve high robustness to strong phase variations at high b-values.

Results

Using Romer-EPTI, we demonstrate distortion-free whole-brain mesoscale in-vivo dMRI at both 3T (500-μm-iso) and 7T (485-μm-iso) for the first time, with high SNR efficiency (e.g., 25 × ), and high image quality free from distortion and slab-boundary artifacts with minimal blurring. Motion experiments demonstrate Romer-EPTI's high motion-robustness and ability to recover sharp images in the presence of motion. Romer-EPTI also demonstrates significant SNR gain and robustness in high b-value (b=5000s/mm2) and time-dependent dMRI.

Conclusion

Romer-EPTI significantly improves SNR, motion-robustness, and image quality, providing a highly efficient acquisition for high-resolution dMRI and microstructure imaging.

Magnetoencephalography Atlas Viewer for Dipole Localization and Viewing

Magnetoencephalography (MEG) is a noninvasive neuroimaging technique widely recognized for epilepsy and tumor mapping. MEG clinical reporting requires a multidisciplinary team, including expert input regarding each dipole's anatomic localization. Here, we introduce a novel tool, the "Magnetoencephalography Atlas Viewer" (MAV), which streamlines this anatomical analysis. The MAV normalizes the patient's Magnetic Resonance Imaging (MRI) to the Montreal Neurological Institute (MNI) space, reverse-normalizes MNI atlases to the native MRI, identifies MEG dipole files, and matches dipoles' coordinates to their spatial location in atlas files. It offers a user-friendly and interactive graphical user interface (GUI) for displaying individual dipoles, groups, coordinates, anatomical labels, and a tri-planar MRI view of the patient with dipole overlays. It evaluated over 273 dipoles obtained in clinical epilepsy subjects. Consensus-based ground truth was established by three neuroradiologists, with a minimum agreement threshold of two. The concordance between the ground truth and MAV labeling ranged from 79% to 84%, depending on the normalization method. Higher concordance rates were observed in subjects with minimal or no structural abnormalities on the MRI, ranging from 80% to 90%. The MAV provides a straightforward MEG dipole anatomic localization method, allowing a nonspecialist to prepopulate a report, thereby facilitating and reducing the time of clinical reporting.

Dynamic parallel transmit diffusion MRI at 7T

Diffusion MRI (dMRI) is inherently limited by SNR. Scanning at 7 T increases intrinsic SNR but 7 T MRI scans suffer from regions of signal dropout, especially in the temporal lobes and cerebellum. We applied dynamic parallel transmit (pTx) to allow whole-brain 7 T dMRI and compared with circularly polarized (CP) pulses in 6 subjects. Subject-specific 2-spoke dynamic pTx pulses were designed offline for 8 slabs covering the brain. We used vendor-provided B0 and B1+ mapping. Spokes positions were set using the Fourier difference approach, and RF coefficients optimized with a Jacobi-matrix high-flip-angle optimizer. Diffusion data were analyzed with FSL. Comparing whole-brain averages for pTx against CP scans: mean flip angle error improved by 15% for excitation (2-spoke-VERSE 15.7° vs CP 18.4°, P = 0.012) and improved by 14% for refocusing (2-spoke-VERSE 39.7° vs CP 46.2°, P = 0.008). Computed spin-echo signal standard deviation improved by 14% (2-spoke-VERSE 0.185 vs 0.214 CP, P = 0.025). Temporal SNR increased by 5.4% (2-spoke-VERSE 8.47 vs CP 8.04, P = 0.004) especially in the inferior temporal lobes. Diffusion fitting uncertainty decreased by 6.2% for first fibers (2-spoke VERSE 0.0655 vs CP 0.0703, P < 0.001) and 1.3% for second fibers (2-spoke VERSE 0.139 vs CP 0.141, P = 0.01). In conclusion, dynamic parallel transmit improves the uniformity of 7 T diffusion-weighted imaging. In future, less restrictive SAR limits for parallel transmit scans are expected to allow further improvements.