Functional integration between brain regions at rest occurs in multiple-frequency bands

Atypical brain function hierarchy in autism spectrum disorder: insights from a novel analytical approach based on neuronal oscillation pattern

Hierarchy is the basic character of the human brain. Neuronal oscillation is one of the fundamental features of brain function, revealing abnormal hierarchical structures in psychiatric disorders from a system-level perspective. However, to date, no research has yet quantified the normal and abnormal brain functional hierarchy based on oscillation patterns. Therefore, this study aimed to quantify brain hierarchy based on neuronal oscillation patterns using the wide-scale information across multiple frequency bands of functional magnetic resonance imaging (fMRI) data and further investigate atypical oscillation patterns in autism spectrum disorder (ASD) at the system level. We analyzed resting-state fMRI data from the Autism Brain Imaging Data Exchange II, including 132 participants with ASD and 132 healthy controls. The energy distribution patterns (EDPs) across frequency bands were calculated for different brain networks using multivariate empirical mode decomposition and Hilbert Transform to represent oscillation patterns. The gradient analysis was applied to quantify the EDP segregation among networks, and the network median distance of gradients was compared between the two groups. The k-means clustering was applied to intuitively verify the atypical EDP in ASD. Across all participants, we observed that the EDPs of different brain regions were spatially coupled to the brain hierarchy. Compared to healthy controls, the ASD exhibited reduced segregation between unimodal and transmodal regions on both energy gradient and clustering analyses, correlating with social deficits. Our results quantitatively confirm that oscillation patterns can reflect the functional segregation among networks and provide novel evidence of the system-level imbalances in neuronal oscillations in ASD.

Mapping Dynamic Metabolic Energy Distribution in Brain Networks using fMRI: A Novel Dynamic Time Warping Framework

Understanding how metabolic energy is distributed across brain networks is essential for elucidating healthy brain function and neurological disorders. Research has established the link between blood flow changes and glucose metabolic processes that fuel neural activity. Here, we introduce a novel framework based on the normalized dynamic time warping algorithm robust to neural temporal variability, enabling reliable insights into metabolic energy demands using functional magnetic resonance imaging data. Our findings indicate that healthy brains maintain balanced energy distribution, whereas imbalances are more pronounced in schizophrenia with links to both positive and negative symptoms, particularly during rapid neural processes. Additionally, we identified a dynamic state that supports the brain criticality theory and is associated with higher-order cognitive abilities, demonstrating our framework's functional and clinical relevance. By linking metabolic energy distribution to neural dynamics, this framework provides a novel way to estimate and quantify the brain's maintenance of functional balance in a broadly applicable manner for studying brain health and disorders.

Dark brain energy: Toward an integrative model of spontaneous slow oscillations

Neural oscillations facilitate the functioning of the human brain in spatial and temporal dimensions at various frequencies. These oscillations feature a universal frequency architecture that is governed by brain anatomy, ensuring frequency specificity remains invariant across different measurement techniques. Initial magnetic resonance imaging (MRI) methodology constrained functional MRI (fMRI) investigations to a singular frequency range, thereby neglecting the frequency characteristics inherent in blood oxygen level-dependent oscillations. With advancements in MRI technology, it has become feasible to decode intricate brain activities via multi-band frequency analysis (MBFA). During the past decade, the utilization of MBFA in fMRI studies has surged, unveiling frequency-dependent characteristics of spontaneous slow oscillations (SSOs) believed to base dark energy in the brain. There remains a dearth of conclusive insights and hypotheses pertaining to the properties and functionalities of SSOs in distinct bands. We surveyed the SSO MBFA studies during the past 15 years to delineate the attributes of SSOs and enlighten their correlated functions. We further proposed a model to elucidate the hierarchical organization of multi-band SSOs by integrating their function, aimed at bridging theoretical gaps and guiding future MBFA research endeavors.

Comparing Intra- and Inter-individual Correlational Brain Connectivity from Functional and Structural Neuroimaging Data

Inferring brain connectivity from inter-individual correlations has been applied across various neuroimaging modalities, including positron emission tomography (PET) and MRI. The variability underlying these inter-individual correlations is generally attributed to factors such as genetics, life experiences, and long-term influences like aging. This study leveraged two unique longitudinal datasets to examine intra-individual correlations of structural and functional brain measures across an extended time span. By focusing on intra-individual correlations, we aimed to minimize individual differences and investigate how aging and state-like effects contribute to brain connectivity patterns. Additionally, we compared intra-individual correlations with inter-individual correlations to better understand their relationship. In the first dataset, which included repeated scans from a single individual over 15 years, we found that intra-individual correlations in both regional homogeneity (ReHo) during resting-state and gray matter volumes (GMV) from structural MRI closely resembled resting-state functional connectivity. However, ReHo correlations were primarily driven by state-like variability, whereas GMV correlations were mainly influenced by aging. The second dataset, comprising multiple participants with longitudinal Fludeoxyglucose (18F) FDG-PET and MRI scans, replicated these findings. Both intra- and inter-individual correlations were strongly associated with resting-state functional connectivity, with functional measures (i.e., ReHo and FDG-PET) exhibiting greater similarity to resting-state connectivity than structural measures. This study demonstrated that controlling for various factors can enhance the interpretability of brain correlation structures. While inter- and intra-individual correlation patterns showed similarities, accounting for additional variables may improve our understanding of inter-individual connectivity measures.

The multi-frequency decomposition entropy learning for nonlinear fMRI data analysis

Functional magnetic resonance imaging (fMRI) have been widely adopted to explore the underlying neural mechanisms between psychiatric disorders which share common neurobiology and clinical manifestations. However, the existing studies mainly focus on linear relationships and ignore nonlinear contributions. To address the above issues, we propose a new method named multi-frequency decomposition entropy (MDE) learning for inferring nonlinear functional connectivity between brain regions. Firstly, the variational mode decomposition was used to divide fMRI data into five groups of frequency. Next, the copula entropy was used to calculate the nonlinear relationship between brain regions in each frequency group, and then the best important nonlinear relationships were screen out by using statistical t-test. Lastly, a gyrus importance index was proposed to reflect the distribution trend of gyri in different frequency groups. The results of applying MDE for the fMRI data analysis of schizophrenia, bipolar disorder, and attention-deficit hyperactivity disorder showed that the difference between the three groups of patient and healthy control is large at the hub nodes, and the nonlinear relationship between the patient groups is weak when they are at the same hub node. In addition, each disease exhibits unique characteristics compared with other diseases and healthy control. In a word, the nonlinear functional connectivity of different frequency groups reflect the differences and commonalities between diseases and reveal possible discriminating biomarkers among mental diseases.

Processing, evaluating, and understanding FMRI data with afni_proc.py

FMRI data are noisy, complicated to acquire, and typically go through many steps of processing before they are used in a study or clinical practice. Being able to visualize and understand the data from the start through the completion of processing, while being confident that each intermediate step was successful, is challenging. AFNI's afni_proc.py is a tool to create and run a processing pipeline for FMRI data. With its flexible features, afni_proc.py allows users to both control and evaluate their processing at a detailed level. It has been designed to keep users informed about all processing steps: it does not just process the data, but also first outputs a fully commented processing script that the users can read, query, interpret, and refer back to. Having this full provenance is important for being able to understand each step of processing; it also promotes transparency and reproducibility by keeping the record of individual-level processing and modeling specifics in a single, shareable place. Additionally, afni_proc.py creates pipelines that contain several automatic self-checks for potential problems during runtime. The output directory contains a dictionary of relevant quantities that can be programmatically queried for potential issues and a systematic, interactive quality control (QC) HTML. All of these features help users evaluate and understand their data and processing in detail. We describe these and other aspects of afni_proc.py here using a set of task-based and resting-state FMRI example commands.

Task-based default mode network connectivity predicts cognitive impairment and negative symptoms in first-episode schizophrenia

Individuals diagnosed with schizophrenia (SZ) demonstrate difficulty distinguishing between internally and externally generated stimuli. These aberrations in "source monitoring" have been theorized as contributing to symptoms of the disorder, including hallucinations and delusions. Altered connectivity within the default mode network (DMN) of the brain has been proposed as a mechanism through which discrimination between self-generated and externally generated events is disrupted. Source monitoring abnormalities in SZ have additionally been linked to impairments in selective attention and inhibitory processing, which are reliably observed via the N100 component of the event-related brain potential elicited during an auditory paired-stimulus paradigm. Given overlapping constructs associated with DMN connectivity and N100 in SZ, the present investigation evaluated relationships between these measures of disorder-related dysfunction and sought to clarify the nature of task-based DMN function in SZ. DMN connectivity and N100 measures were assessed using EEG recorded from SZ during their first episode of illness (N = 52) and demographically matched healthy comparison participants (N = 25). SZ demonstrated less evoked theta-band connectivity within DMN following presentation of pairs of identical auditory stimuli than HC. Greater DMN connectivity among SZ was associated with better performance on measures of sustained attention (p = .03) and working memory (p = .09), as well as lower severity of negative symptoms, though it was not predictive of N100 measures. Together, present findings provide EEG evidence of lower task-based connectivity among first-episode SZ, reflecting disruptions of DMN functions that support cognitive processes. Attentional processes captured by N100 appear to be supported by different neural mechanisms.

From convergence insufficiency to functional reorganization: A longitudinal randomized controlled trial of treatment-induced connectivity plasticity

INTRODUCTION

Convergence Insufficiency (CI) is the most prevalent oculomotor dysfunction of binocular vision that negatively impacts quality of life when performing visual near tasks. Decreased resting-state functional connectivity (RSFC) is reported in the CI participants compared to binocularly normal control participants. Studies report that therapeutic interventions such as office-based vergence and accommodative therapy (OBVAT) can improve CI participants' clinical signs, visual symptoms, and task-related functional activity. However, longitudinal studies investigating the RSFC changes after such treatments in participants with CI have not been conducted. This study aimed to investigate the neural basis of OBVAT using RSFC in CI participants compared to the placebo treatment to understand how OBVAT improves visual function and symptoms.

METHODS

A total of 51 CI participants between 18 and 35 years of age were included in the study and randomly allocated to receive either 12 one-hour sessions of OBVAT or placebo treatment for 6 to 8 weeks (1 to 2 sessions per week). Resting-state functional magnetic resonance imaging and clinical assessments were evaluated at baseline and outcome for each treatment group. Region of interest (ROI) analysis was conducted in nine ROIs of the oculomotor vergence network, including the following: cerebellar vermis (CV), frontal eye fields (FEF), supplementary eye fields (SEF), parietal eye fields (PEF), and primary visual cortices (V1). Paired t-tests assessed RSFC changes in each group. A linear regression analysis was conducted for significant ROI pairs in the group-level analysis for correlations with clinical measures.

RESULTS

Paired t-test results showed increased RSFC in 10 ROI pairs after the OBVAT but not placebo treatment (p < 0.05, false discovery rate corrected). These ROI pairs included the following: Left (L)-SEF-Right (R)-V1, L-SEF-CV, R-SEF-R-PEF, R-SEF-L-V1, R-SEF-R-V1, R-SEF-CV, R-PEF-CV, L-V1-CV, R-V1-CV, and L-V1-R-V1. Significant correlations were observed between the RSFC strength of the R-SEF-R-PEF ROI pair and the following clinical visual function parameters: positive fusional vergence and near point of convergence (p < 0.05).

CONCLUSION

OBVAT, but not placebo treatment, increased the RSFC in the ROIs of the oculomotor vergence network, which was correlated with the improvements in the clinical measures of the CI participants.

Multiband Group Independent Component Analysis: Unveiling Frequency-Dependent Dynamics of Functional Connectivity in Group-Level fMRI Analyses

Group independent component analysis (ICA) offers a method for decomposing fMRI data from multiple subjects into spatially independent maps and associated time courses. Traditionally, group ICA is applied to full-band fMRI data, with typically sampling rates between 0.25-1 Hz. In this paper, we introduce a novel approach known as "multiband group ICA." This method involves the application of bandpass filters to segment the fMRI data into distinct subbands, followed by the application of blind group ICA to all subbands of multisubject fMRI data, followed by back-reconstruction of individual subband information.To assess the feasibility and efficacy of our method, we utilized a bandpass filter to divide the fMRI data into two specific subbands: low frequency and high frequency. Our results not only showcase a substantial distinction in spatial maps and time courses for task-related components but also provide noteworthy insights into network-specific functional network connectivity (FNC) patterns, particularly within the visual network. Notably, motor regions exhibit predominant power in the lower frequency range, while visual regions span both low and high frequencies. Additionally, certain regions, such as the bilateral insular regions, display exclusive engagement in the high-frequency domain. Furthermore, low-frequency task-related components exhibit anticorrelation, contrasting with the strong correlation observed in high-frequency task-related components. This approach enhances our understanding of the frequency-dependent dynamics of functional connectivity in group-level fMRI analyses, offering valuable insights for future studies.

Pathology of pain and its implications for therapeutic interventions

Pain is estimated to affect more than 20% of the global population, imposing incalculable health and economic burdens. Effective pain management is crucial for individuals suffering from pain. However, the current methods for pain assessment and treatment fall short of clinical needs. Benefiting from advances in neuroscience and biotechnology, the neuronal circuits and molecular mechanisms critically involved in pain modulation have been elucidated. These research achievements have incited progress in identifying new diagnostic and therapeutic targets. In this review, we first introduce fundamental knowledge about pain, setting the stage for the subsequent contents. The review next delves into the molecular mechanisms underlying pain disorders, including gene mutation, epigenetic modification, posttranslational modification, inflammasome, signaling pathways and microbiota. To better present a comprehensive view of pain research, two prominent issues, sexual dimorphism and pain comorbidities, are discussed in detail based on current findings. The status quo of pain evaluation and manipulation is summarized. A series of improved and innovative pain management strategies, such as gene therapy, monoclonal antibody, brain-computer interface and microbial intervention, are making strides towards clinical application. We highlight existing limitations and future directions for enhancing the quality of preclinical and clinical research. Efforts to decipher the complexities of pain pathology will be instrumental in translating scientific discoveries into clinical practice, thereby improving pain management from bench to bedside.

Causal Forest Machine Learning Analysis of Parkinson's Disease in Resting-State Functional Magnetic Resonance Imaging

In recent years, Artificial Intelligence has been used to assist healthcare professionals in detecting and diagnosing neurodegenerative diseases. In this study, we propose a methodology to analyze functional Magnetic Resonance Imaging signals and perform classification between Parkinson's disease patients and healthy participants using Machine Learning algorithms. In addition, the proposed approach provides insights into the brain regions affected by the disease. The functional Magnetic Resonance Imaging from the PPMI and 1000-FCP datasets were pre-processed to extract time series from 200 brain regions per participant, resulting in 11,600 features. Causal Forest and Wrapper Feature Subset Selection algorithms were used for dimensionality reduction, resulting in a subset of features based on their heterogeneity and association with the disease. We utilized Logistic Regression and XGBoost algorithms to perform PD detection, achieving 97.6% accuracy, 97.5% F1 score, 97.9% precision, and 97.7%recall by analyzing sets with fewer than 300 features in a population including men and women. Finally, Multiple Correspondence Analysis was employed to visualize the relationships between brain regions and each group (women with Parkinson, female controls, men with Parkinson, male controls). Associations between the Unified Parkinson's Disease Rating Scale questionnaire results and affected brain regions in different groups were also obtained to show another use case of the methodology. This work proposes a methodology to (1) classify patients and controls with Machine Learning and Causal Forest algorithm and (2) visualize associations between brain regions and groups, providing high-accuracy classification and enhanced interpretability of the correlation between specific brain regions and the disease across different groups.

Brain local stability and network flexibility of table tennis players: a 7T MRI study

Table tennis players have adaptive visual and sensorimotor networks, which are the key brain regions to acquire environmental information and generate motor output. This study examined 20 table tennis players and 21 control subjects through ultrahigh field 7 Tesla magnetic resonance imaging. First, we measured percentage amplitude of fluctuation across five different frequency bands and found that table tennis players had significantly lower percentage amplitude of fluctuation values than control subjects in 18 brain regions, suggesting enhanced stability of spontaneous brain fluctuation amplitudes in visual and sensorimotor networks. Functional connectional analyses revealed increased static functional connectivity between two sensorimotor nodes and other frontal-parietal regions among table tennis players. Additionally, these players displayed enhanced dynamic functional connectivity coupled with reduced static connectivity between five nodes processing visual and sensory information input, and other large-scale cross-regional areas. These findings highlight that table tennis players undergo neural adaptability through a dual mechanism, characterized by global stability in spontaneous brain fluctuation amplitudes and heightened flexibility in visual sensory networks. Our study offers novel insights into the mechanisms of neural adaptability in athletes, providing a foundation for future efforts to enhance cognitive functions in diverse populations, such as athletes, older adults, and individuals with cognitive impairments.

Neurovascular coupling in eye-open-eye-close task and resting state: Spectral correspondence between concurrent EEG and fMRI

Neurovascular coupling serves as an essential neurophysiological mechanism in functional neuroimaging, which is generally presumed to be robust and invariant across different physiological states, encompassing both task engagement and resting state. Nevertheless, emerging evidence suggests that neurovascular coupling may exhibit state dependency, even in normal human participants. To investigate this premise, we analyzed the cross-frequency spectral correspondence between concurrently recorded electroencephalography (EEG) and functional magnetic resonance imaging (fMRI) data, utilizing them as proxies for neurovascular coupling during the two conditions: an eye-open-eye-close (EOEC) task and a resting state. We hypothesized that given the state dependency of neurovascular coupling, EEG-fMRI spectral correspondences would change between the two conditions in the visual system. During the EOEC task, we observed a negative phase-amplitude-coupling (PAC) between EEG alpha-band and fMRI visual activity. Conversely, in the resting state, a pronounced amplitude-amplitude-coupling (AAC) emerged between EEG and fMRI signals, as evidenced by the spectral correspondence between the EEG gamma-band of the midline occipital channel (Oz) and the high-frequency fMRI signals (0.15-0.25 Hz) in the visual network. This study reveals distinct scenarios of EEG-fMRI spectral correspondence in healthy participants, corroborating the state-dependent nature of neurovascular coupling.

Spectral graph model for fMRI: a biophysical, connectivity-based generative model for the analysis of frequency-resolved resting state fMRI

Resting state functional MRI (rs-fMRI) is a popular and widely used technique to explore the brain's functional organization and to examine if it is altered in neurological or mental disorders. The most common approach for its analysis targets the measurement of the synchronized fluctuations between brain regions, characterized as functional connectivity (FC), typically relying on pairwise correlations in activity across different brain regions. While hugely successful in exploring state- and disease-dependent network alterations, these statistical graph theory tools suffer from two key limitations. First, they discard useful information about the rich frequency content of the fMRI signal. The rich spectral information now achievable from advances in fast multiband acquisitions is consequently being under-utilized. Second, the analyzed FCs are phenomenological without a direct neurobiological underpinning in the underlying structures and processes in the brain. There does not currently exist a complete generative model framework for whole brain resting fMRI that is informed by its underlying biological basis in the structural connectome. Here we propose that a different approach can solve both challenges at once: the use of an appropriately realistic yet parsimonious biophysical signal generation model followed by graph spectral (i.e. eigen) decomposition. We call this model a Spectral Graph Model (SGM) for fMRI, using which we can not only quantify the structure-function relationship in individual subjects, but also condense the variable and individual-specific repertoire of fMRI signal's spectral and spatial features into a small number of biophysically-interpretable parameters. We expect this model-based inference of rs-fMRI that seamlessly integrates with structure can be used to examine state and trait characteristics of structure-function relations in a variety of brain disorders.

Longitudinal Functional Connectome in Pediatric Concussion: An Advancing Concussion Assessment in Pediatrics Study

Advanced magnetic resonance imaging (MRI) techniques indicate that concussion (i.e., mild traumatic brain injury) disrupts brain structure and function in children. However, the functional connectivity of brain regions within global and local networks (i.e., functional connectome) is poorly understood in pediatric concussion. This prospective, longitudinal study addressed this gap using data from the largest neuroimaging study of pediatric concussion to date to study the functional connectome longitudinally after concussion as compared with mild orthopedic injury (OI). Children and adolescents (n = 967) 8-16.99 years with concussion or mild OI were recruited from pediatric emergency departments within 48 h post-injury. Pre-injury and 1-month post-injury symptom ratings were used to classify concussion with or without persistent symptoms based on reliable change. Subjects completed a post-acute (2-33 days) and chronic (3 or 6 months via random assignment) MRI scan. Graph theory metrics were derived from 918 resting-state functional MRI scans in 585 children (386 concussion/199 OI). Linear mixed-effects modeling was performed to assess group differences over time, correcting for multiple comparisons. Relative to OI, the global clustering coefficient was reduced at 3 months post-injury in older children with concussion and in females with concussion and persistent symptoms. Time post-injury and sex moderated group differences in local (regional) network metrics of several brain regions, including degree centrality, efficiency, and clustering coefficient of the angular gyrus, calcarine fissure, cuneus, and inferior occipital, lingual, middle occipital, post-central, and superior occipital gyrus. Relative to OI, degree centrality and nodal efficiency were reduced post-acutely, and nodal efficiency and clustering coefficient were reduced chronically after concussion (i.e., at 3 and 6 months post-injury in females; at 6 months post-injury in males). Functional network alterations were more robust and widespread chronically as opposed to post-acutely after concussion, and varied by sex, age, and symptom recovery at 1-month post-injury. Local network segregation reductions emerged globally (across the whole brain network) in older children and in females with poor recovery chronically after concussion. Reduced functioning between neighboring regions could negatively disrupt specialized information processing. Local network metric alterations were demonstrated in several posterior regions that are involved in vision and attention after concussion relative to OI. This indicates that functioning of superior parietal and occipital regions could be particularly susceptibile to the effects of concussion. Moreover, those regional alterations were especially apparent at later time periods post-injury, emerging after post-concussive symptoms resolved in most and persisted up to 6 months post-injury, and differed by biological sex. This indicates that neurobiological changes continue to occur up to 6 months after pediatric concussion, although changes emerge earlier in females than in males. Changes could reflect neural compensation mechanisms.

Distinct functional subnetworks of cognitive domains in older adults with minor cognitive deficits

Although past research has established a relationship between functional connectivity and cognitive function, less is known about which cognitive domains are associated with which specific functional networks. This study investigated associations between functional connectivity and global cognitive function and performance in the domains of memory, executive function and psychomotor speed in 166 older adults aged 75-91 years (mean = 80.3 ± 3.8) with minor cognitive deficits (Mini-Mental State Examination scores between 21 and 27). Functional connectivity was assessed within 10 standard large-scale resting-state networks and on a finer spatial resolution between 300 nodes in a functional connectivity matrix. No domain-specific associations with mean functional connectivity within large-scale resting-state networks were found. Node-level analysis revealed that associations between functional connectivity and cognitive performance differed across cognitive functions in strength, location and direction. Specific subnetworks of functional connections were found for each cognitive domain in which higher connectivity between some nodes but lower connectivity between other nodes were related to better cognitive performance. Our findings add to a growing body of literature showing differential sensitivity of functional connections to specific cognitive functions and may be a valuable resource for hypothesis generation of future studies aiming to investigate specific cognitive dysfunction with resting-state functional connectivity in people with beginning cognitive deficits.

Auditory inputs modulate intrinsic neuronal timescales during sleep

Functional magnetic resonance imaging (fMRI) studies have demonstrated that intrinsic neuronal timescales (INT) undergo modulation by external stimulation during consciousness. It remains unclear if INT keep the ability for significant stimulus-induced modulation during primary unconscious states, such as sleep. This fMRI analysis addresses this question via a dataset that comprises an awake resting-state plus rest and stimulus states during sleep. We analyzed INT measured via temporal autocorrelation supported by median frequency (MF) in the frequency-domain. Our results were replicated using a biophysical model. There were two main findings: (1) INT prolonged while MF decreased from the awake resting-state to the N2 resting-state, and (2) INT shortened while MF increased during the auditory stimulus in sleep. The biophysical model supported these results by demonstrating prolonged INT in slowed neuronal populations that simulate the sleep resting-state compared to an awake state. Conversely, under sine wave input simulating the stimulus state during sleep, the model's regions yielded shortened INT that returned to the awake resting-state level. Our results highlight that INT preserve reactivity to stimuli in states of unconsciousness like sleep, enhancing our understanding of unconscious brain dynamics and their reactivity to stimuli.

Altered Large-Scale Resting-State Functional Network Connectivity in Convergence Insufficiency Young Adults Compared With Binocularly Normal Controls

Purpose

To investigate the underlying resting-state functional connectivity (RSFC) of symptomatic convergence insufficiency (CI) compared with binocularly normal controls (BNC) using functional magnetic resonance imaging (fMRI) under The Convergence Insufficiency Neuro‑mechanism Adult Population Study (NCT03593031).

Methods

A total of 101 participants were eligible for this study. After removing datasets with motion artifacts, 49 CI and 47 BNC resting-state functional magnetic resonance imaging datasets were analyzed. CI was diagnosed with the following signs: (1) receded near point of convergence of 6 cm or greater, (2) decreased positive fusional vergence of less than 15∆ or failing Sheard's criteria of twice the near phoria, (3) near phoria of at least 4∆ more exophoric compared with the distance phoria, and (4) symptoms using the Convergence Insufficiency Symptom Survey (score of ≥21). RSFC was assessed using a group-level independent components analysis and dual regression. A behavioral correlation analysis using linear regression method was performed between clinical measures and RSFC using the significant difference between the CI and BNC.

Results

On average, a decreased RSFC was observed within the frontoparietal network, default mode network and visual network in patients with CI, compared with the participants with BNC (P < 0.05, corrected for multiple comparisons). The default mode network RSFC strength was significantly correlated with the PFV, near point of convergence, and difference between the horizontal phoria at near compared with far (P < 0.05).

Conclusions

Results support altered RSFC in patients with CI compared with participants with BNC and suggest that these differences in underlying neurophysiology may in part be in connection with the differences in optometric visual function used to diagnose CI.

Comparison of hemodynamic response functions obtained from resting-state functional MRI and invasive electrophysiological recordings in rats

Resting-state functional MRI (rs-fMRI) is a popular technology that has enriched our understanding of brain and spinal cord functioning, including how different regions communicate (connectivity). But fMRI is an indirect measure of neural activity capturing blood hemodynamics. The hemodynamic response function (HRF) interfaces between the unmeasured neural activity and measured fMRI time series. The HRF is variable across brain regions and individuals, and is modulated by non-neural factors. Ignoring this HRF variability causes errors in FC estimates. Hence, it is crucial to reliably estimate the HRF from rs-fMRI data. Robust techniques have emerged to estimate the HRF from fMRI time series. Although such techniques have been validated non-invasively using simulated and empirical fMRI data, thorough invasive validation using simultaneous electrophysiological recordings, the gold standard, has been elusive. This report addresses this gap in the literature by comparing HRFs derived from invasive intracranial electroencephalogram recordings with HRFs estimated from simultaneously acquired fMRI data in six epileptic rats. We found that the HRF shape parameters (HRF amplitude, latency and width) were not significantly different (p>0.05) between ground truth and estimated HRFs. In the single pathological region, the HRF width was marginally significantly different (p=0.03). Our study provides preliminary invasive validation for the efficacy of the HRF estimation technique in reliably estimating the HRF non-invasively from rs-fMRI data directly. This has a notable impact on rs-fMRI connectivity studies, and we recommend that HRF deconvolution be performed to minimize HRF variability and improve connectivity estimates.

Frequency Dependent Changes of Regional Homogeneity in Children with Growth Hormone Deficiency

Abnormal spontaneous neural activity in children with growth hormone deficiency (GHD) has been found in previous resting-state functional magnetic resonance imaging (rs-fMRI) studies. Nevertheless, the spontaneous neural activity of GHD in different frequency bands is still unclear. Here, we combined rs-fMRI and regional homogeneity (ReHo) methods to analyze the spontaneous neural activity of 26 GHD children and 15 healthy controls (HCs) with age- and sex-matching in four frequency bands: slow-5 (0.014-0.031 Hz), slow-4 (0.031-0.081 Hz), slow-3 (0.081-0.224 Hz), and slow-2 (0.224-0.25 Hz). In the slow-5 band, GHD children compared with HCs displayed higher ReHo in the left dorsolateral part of the superior frontal gyrus, triangular part of the inferior frontal gyrus, precentral gyrus and middle frontal gyrus, and right angular gyrus, while lower ReHo in the right precentral gyrus, and several medial orbitofrontal regions. In the slow-4 band, GHD children relative to HCs revealed increased ReHo in the right middle temporal gyrus, whereas reduced ReHo in the left superior parietal gyrus, right middle occipital gyrus, and bilateral medial parts of the superior frontal gyrus. In the slow-2 band, compared with HCs, GHD children showed increased ReHo in the right anterior cingulate gyrus, and several prefrontal regions, while decreased ReHo in the left middle occipital gyrus, and right fusiform gyrus and anterior cingulate gyrus. Our findings demonstrate that regional brain activity in GHD children exhibits extensive abnormalities, and these abnormalities are related to specific frequency bands, which may provide bases for understanding its pathophysiology significance.

Deficits of Visual Cortex Function in Acute Acquired Concomitant Esotropia Patients

Purpose

The purpose of this study was to explore the cortical deficits of patients with acquired concomitant esotropia (AACE) using the resting-state functional magnetic resonance imaging (rs-fMRI) technique.

Methods

Rs-fMRI signals from 25 patients with AACE and 25 matched controls were collected. The repeated-measures analysis of variance (RM-ANOVA) test and two-sample t-test were used to investigate statistical differences of the amplitudes of low-frequency fluctuation (ALFF) signals and correlation analysis was performed to validate the relationship of signal change and clinical features.

Results

The AACE group showed decreased ALFF in both hemispheres symmetrically (t = 0.38, P = 0.71), with peak t in both middle occipital gyrus. The ALFF signal from the upper left inferior frontal gyrus was negatively correlated with the age of onset (r = 0.62, P = 0.0008), and the ALFF signal from the right superior temporal gyrus was negatively correlated with the near work hours (r = 0.63, P = 0.0008). The ALFF signal in the left fusiform gyrus was positively correlated with both near (r = 0.48, P = 0.01) and far (r = 0.44, P = 0.03) deviation, whereas it was only positively correlated with far deviation (r = 0.44, P = 0.03) in the right. Besides, the age of onset and the near work hour were independent factors of signal changes.

Conclusions

Using the ALFF signal of rs-fMRI, we found functional deficits in the primary visual cortex and dorsal pathway in patients with AACE. There were functional changes in the fusiform gyrus, and the greater the deviation angle, the higher the changing level. These findings reveal the association of AACE and the visual center, giving us more clues about the treatment of AACE.

Task-Based and Resting-State Cortical Functional Differences After Spinal Cord Injury: A Pilot Functional Near-Infrared Spectroscopy Study

Brain reorganization following spinal cord injury (SCI) has been well-established using animal and human studies. Yet, much is unknown regarding functional recovery and adverse secondary outcomes after SCI. Functional near-infrared spectroscopy (fNIRS) is a neuroimaging technique that offers methodological flexibility in a real-world setting. We used fNIRS to examine the cortical functional differences between 12 males with thoracolumbar SCI (46.41 ± 11.09 years of age) and 12 healthy males (47.61 ± 11.94 years of age) during resting state and task conditions-bilateral finger tapping (FT), mental imagery of bilateral FT with action observation (FTI+AO), and bilateral ankle tapping (AT). We found an overall decrease in hemodynamic response of the SCI group during all three task conditions. Task modulated functional connectivity (FC) computed using beta series correlation technique was compared using independent sample t-tests at α = 0.05. Connectivity between the right mediolateral sensorimotor network (SMN) and the right medial SMN was reduced during the FT task in SCI. A mixed analysis of variance revealed that the FC within the right mediolateral SMN was reduced during FT but preserved during FTI+AO (i.e., comparable to controls) in the SCI group. Lower FC of these regions was associated with longer injury durations. Additionally, we found a general decrease in resting state FC of the SCI group, specifically in the Slow-3 frequency range (0.073 to 0.1 Hz). These results, though preliminary, are consistent with past studies and highlight the potential of fNIRS in SCI and rehabilitative research.

Connectivity gradients in spontaneous brain activity at multiple frequency bands

The intrinsic organizational structure of the brain is reflected in spontaneous brain oscillations. Its functional integration and segregation hierarchy have been discovered in space by leveraging gradient approaches to low-frequency functional connectivity. This hierarchy of brain oscillations has not yet been fully understood, since previous studies have mainly concentrated on the brain oscillations from a single limited frequency range (~ 0.01-0.1 Hz). In this work, we extended the frequency range and performed gradient analysis across multiple frequency bands of fast resting-state fMRI signals from the Human Connectome Project and condensed a frequency-rank cortical map of the highest gradient. We found that the coarse skeletons of the functional organization hierarchy are generalizable across the multiple frequency bands. Beyond that, the highest integration levels of connectivity vary in the frequency domain across different large-scale brain networks. These findings are replicated in another independent dataset and demonstrated that different brain networks can integrate information at varying rates, indicating the significance of examining the intrinsic architecture of spontaneous brain activity from the perspective of multiple frequency bands.

Power spectral analysis can determine language laterality from resting-state functional MRI data in healthy controls

BACKGROUND AND PURPOSE

Resting-state functional magnetic resonance imaging (rsfMRI) has been proposed as an alternative to task-based fMRI including clinical situations such as preoperative brain tumor planning, due to advantages including ease of performance and time savings. However, one of its drawbacks is the limited ability to accurately lateralize language function.

METHODS

Using the rsfMRI data of healthy controls, we carried out a power spectra analysis on three regions of interest (ROIs): Broca's area (BA) in the frontal cortex for language, hand motor (HM) area in the primary motor cortex, and the primary visual cortex (V1). Spike removal, motion correction, linear trend removal, and spatial smoothing were applied. Spontaneous low-frequency fluctuations (0.01-0.1 Hz) were filtered to enable functional integration.

RESULTS

BA showed greater power on the left hemisphere relative to the right (p = .0055), while HM (p = .1563) and V1 (p = .4681) were not statistically significant. A novel index, termed the power laterality index (PLI), computed to estimate the degree of power lateralization for each brain region, revealed a statistically significant difference between BA and V1 (p < .00001), where V1 was used as a control since the primary visual cortex does not lateralize. Validation studies used to compare PLI to a laterality index computed using phonemic fluency, a task-based, language fMRI paradigm, demonstrated good correlation.

CONCLUSIONS

The power spectra for BA revealed left language lateralization, which was not replicated in HM or V1. This work demonstrates the feasibility and validity of an ROI-based power spectra analysis on rsfMRI data for language lateralization.

Resting State Network Segregation Modulates Age-Related Differences in Language Production

Older adults typically exhibit decline in language production. However, how the brain supports or fails to support these processes is unclear. Moreover, there are competing hypotheses about the nature of age-related neural changes and whether age-related increases in neural activity reflect compensation or a decline in neural efficiency. In the current study, we investigated the neural bases of language production focusing on resting state functional connectivity. We hypothesized that language production performance, functional connectivity, and their relationship would differ as a function of age. Consistent with prior work, older age was associated with worse language production performance. Functional connectivity analyses showed that network segregation within the left hemisphere language network was maintained across adulthood. However, increased age was associated with lower whole brain network segregation. Moreover, network segregation was related to language production ability. In both network analyses, there were significant interactions with age-higher network segregation was associated with better language production abilities for younger and middle-aged adults, but not for older adults. Interestingly, there was a stronger relationship between language production and the whole brain network segregation than between production and the language network. These results highlight the utility of network segregation measures as an index of brain function, with higher network segregation associated with better language production ability. Moreover, these results are consistent with stability in the left hemisphere language network across adulthood and suggest that dedifferentiation among brain networks, outside of the language network, is a hallmark of aging and may contribute to age-related language production difficulties.

Functional ultrasound detects frequency-specific acute and delayed S-ketamine effects in the healthy mouse brain

Introduction

S-ketamine has received great interest due to both its antidepressant effects and its potential to induce psychosis when administered subchronically. However, no studies have investigated both its acute and delayed effects using in vivo small-animal imaging. Recently, functional ultrasound (fUS) has emerged as a powerful alternative to functional magnetic resonance imaging (fMRI), outperforming it in sensitivity and in spatiotemporal resolution. In this study, we employed fUS to thoroughly characterize acute and delayed S-ketamine effects on functional connectivity (FC) within the same cohort at slow frequency bands ranging from 0.01 to 1.25 Hz, previously reported to exhibit FC.

Methods

We acquired fUS in a total of 16 healthy C57/Bl6 mice split in two cohorts (n = 8 received saline, n = 8 S-ketamine). One day after the first scans, performed at rest, the mice received the first dose of S-ketamine during the second measurement, followed by four further doses administered every 2 days. First, we assessed FC reproducibility and reliability at baseline in six frequency bands. Then, we investigated the acute and delayed effects at day 1 after the first dose and at day 9, 1 day after the last dose, for all bands, resulting in a total of four fUS measurements for every mouse.

Results

We found reproducible (r > 0.9) and reliable (r > 0.9) group-average readouts in all frequency bands, only the 0.01-0.27 Hz band performing slightly worse. Acutely, S-ketamine induced strong FC increases in five of the six bands, peaking in the 0.073-0.2 Hz band. These increases comprised both cortical and subcortical brain areas, yet were of a transient nature, FC almost returning to baseline levels towards the end of the scan. Intriguingly, we observed robust corticostriatal FC decreases in the fastest band acquired (0.75 Hz-1.25  Hz). These changes persisted to a weaker extent after 1 day and at this timepoint they were accompanied by decreases in the other five bands as well. After 9 days, the decreases in the 0.75-1.25 Hz band were maintained, however no changes between cohorts could be detected in any other bands.

Discussion

In summary, the study reports that acute and delayed ketamine effects in mice are not only dissimilar but have different directionalities in most frequency bands. The complementary readouts of the employed frequency bands recommend the use of fUS for frequency-specific investigation of pharmacological effects on FC.

A schizophrenia study based on multi-frequency dynamic functional connectivity analysis of fMRI

At present, fMRI studies mainly focus on the entire low-frequency band (0. 01-0.08 Hz). However, the neuronal activity is dynamic, and different frequency bands may contain different information. Therefore, a novel multi-frequency-based dynamic functional connectivity (dFC) analysis method was proposed in this study, which was then applied to a schizophrenia study. First, three frequency bands (Conventional: 0.01-0.08 Hz, Slow-5: 0.0111-0.0302 Hz, and Slow-4: 0.0302-0.0820 Hz) were obtained using Fast Fourier Transform. Next, the fractional amplitude of low-frequency fluctuations was used to identify abnormal regions of interest (ROIs) of schizophrenia, and dFC among these abnormal ROIs was implemented by the sliding time window method at four window-widths. Finally, recursive feature elimination was employed to select features, and the support vector machine was applied for the classification of patients with schizophrenia and healthy controls. The experimental results showed that the proposed multi-frequency method (Combined: Slow-5 and Slow-4) had a better classification performance compared with the conventional method at shorter sliding window-widths. In conclusion, our results revealed that the dFCs among the abnormal ROIs varied at different frequency bands and the efficiency of combining multiple features from different frequency bands can improve classification performance. Therefore, it would be a promising approach for identifying brain alterations in schizophrenia.

PhiPipe: A multi-modal MRI data processing pipeline with test-retest reliability and predicative validity assessments

Magnetic resonance imaging (MRI) has been one of the primary instruments to measure the properties of the human brain non-invasively in vivo. MRI data generally needs to go through a series of processing steps (i.e., a pipeline) before statistical analysis. Currently, the processing pipelines for multi-modal MRI data are still rare, in contrast to single-modal pipelines. Furthermore, the reliability and validity of the output of the pipelines are critical for the MRI studies. However, the reliability and validity measures are not available or adequate for almost all pipelines. Here, we present PhiPipe, a multi-modal MRI processing pipeline. PhiPipe could process T1-weighted, resting-state BOLD, and diffusion-weighted MRI data and generate commonly used brain features in neuroimaging. We evaluated the test-retest reliability of PhiPipe's brain features by computing intra-class correlations (ICC) in four public datasets with repeated scans. We further evaluated the predictive validity by computing the correlation of brain features with chronological age in three public adult lifespan datasets. The multivariate reliability and predictive validity of the PhiPipe results were also evaluated. The results of PhiPipe were consistent with previous studies, showing comparable or better reliability and validity when compared with two popular single-modality pipelines, namely DPARSF and PANDA. The publicly available PhiPipe provides a simple-to-use solution to multi-modal MRI data processing. The accompanied reliability and validity assessments could help researchers make informed choices in experimental design and statistical analysis. Furthermore, this study provides a framework for evaluating the reliability and validity of image processing pipelines.

Resting-State Functional MRI/PET Profile as a Potential Alternative to Tri-Modality EEG-MR/PET Imaging: An Exploratory Study in Drug-Refractory Epilepsy

Objective The study explores whether the epileptic networks associate with predetermined seizure onset zone (SOZ) identified from other modalities such as electroencephalogram/video electroencephalogram/structural MRI (EEG/VEEG/sMRI) and with the degree of resting-state functional MRI/positron emission tomography (RS-fMRI/PET) coupling. Here, we have analyzed the subgroup of patients who reported having a seizure on the day of scan as postictal cases and compared the findings with interictal cases (seizure-free interval).

Methods We performed independent component analysis (ICA) on RS-fMRI and 20 ICA were hand-labeled as large scale, noise, downstream, and epilepsy networks (Epinets) based on their profile in spatial, time series, and power spectrum domains. We had a total of 43 cases, with 4 cases in the postictal group (100%). Of 39 cases, 14 cases did not yield any Epinet and 25 cases (61%) were analyzed for the final study. The analysis was done patient-wise and correlated with predetermined SOZ.

Results The yield of finding Epinets on RS-fMRI is more during the postictal period than in the interictal period, although PET and RS-fMRI spatial, time series, and power spectral patterns were similar in both these subgroups. Overlaps between large-scale and downstream networks were noted, indicating that epilepsy propagation can involve large-scale cognition networks. Lateralization to SOZ was noted as blood oxygen level–dependent activation and correlated with sMRI/PET findings. Postoperative surgical failure cases showed residual Epinet profile.

Conclusion RS-fMRI may be a viable option for trimodality imaging to obtain simultaneous physiological information at the functional network and metabolic level.

Frequency-specific brain network architecture in resting-state fMRI

The analysis of brain function in resting-state network (RSN) models, ascertained through the functional connectivity pattern of resting-state functional magnetic resonance imaging (rs-fMRI), is sufficiently powerful for studying large-scale functional integration of the brain. However, in RSN-based research, the network architecture has been regarded as the same through different frequency bands. Thus, here, we aimed to examined whether the network architecture changes with frequency. The blood oxygen level-dependent (BOLD) signal was decomposed into four frequency bands-ranging from 0.007 to 0.438 Hz-and the clustering algorithm was applied to each of them. The best clustering number was selected for each frequency band based on the overlap ratio with task activation maps. The results demonstrated that resting-state BOLD signals exhibited frequency-specific network architecture; that is, the networks finely subdivided in the lower frequency bands were integrated into fewer networks in higher frequency bands rather than reconfigured, and the default mode network and networks related to perception had sufficiently strong architecture to survive in an environment with a lower signal-to-noise ratio. These findings provide a novel framework to enable improved understanding of brain function through the multiband frequency analysis of ultra-slow rs-fMRI data.

Quality control practices in FMRI analysis: Philosophy, methods and examples using AFNI

Quality control (QC) is a necessary, but often an under-appreciated, part of FMRI processing. Here we describe procedures for performing QC on acquired or publicly available FMRI datasets using the widely used AFNI software package. This work is part of the Research Topic, "Demonstrating Quality Control (QC) Procedures in fMRI." We used a sequential, hierarchical approach that contained the following major stages: (1) GTKYD (getting to know your data, esp. its basic acquisition properties), (2) APQUANT (examining quantifiable measures, with thresholds), (3) APQUAL (viewing qualitative images, graphs, and other information in systematic HTML reports) and (4) GUI (checking features interactively with a graphical user interface); and for task data, and (5) STIM (checking stimulus event timing statistics). We describe how these are complementary and reinforce each other to help researchers stay close to their data. We processed and evaluated the provided, publicly available resting state data collections (7 groups, 139 total subjects) and task-based data collection (1 group, 30 subjects). As specified within the Topic guidelines, each subject's dataset was placed into one of three categories: Include, exclude or uncertain. The main focus of this paper, however, is the detailed description of QC procedures: How to understand the contents of an FMRI dataset, to check its contents for appropriateness, to verify processing steps, and to examine potential quality issues. Scripts for the processing and analysis are freely available.

Robust dynamic brain coactivation states estimated in individuals

A confluence of evidence indicates that brain functional connectivity is not static but rather dynamic. Capturing transient network interactions in the individual brain requires a technology that offers sufficient within-subject reliability. Here, we introduce an individualized network-based dynamic analysis technique and demonstrate that it is reliable in detecting subject-specific brain states during both resting state and a cognitively challenging language task. We evaluate the extent to which brain states show hemispheric asymmetries and how various phenotypic factors such as handedness and gender might influence network dynamics, discovering a right-lateralized brain state that occurred more frequently in men than in women and more frequently in right-handed versus left-handed individuals. Longitudinal brain state changes were also shown in 42 patients with subcortical stroke over 6 months. Our approach could quantify subject-specific dynamic brain states and has potential for use in both basic and clinical neuroscience research.

Alcohol use is associated with affective and interoceptive network alterations in bipolar disorder

INTRODUCTION

Alcohol use in bipolar disorder (BD) is associated with mood lability and negative illness trajectory, while also impacting functional networks related to emotion, cognition, and introspection. The adverse impact of alcohol use in BD may be explained by its additive effects on these networks, thereby contributing to a poorer clinical outcome.

METHODS

Forty BD-I (DSM-IV-TR) and 46 psychiatrically healthy controls underwent T1 and resting state functional MRI scanning and the Alcohol Use Disorders Identification Test-Consumption (AUDIT-C) to assess alcohol use. Functional images were decomposed using spatial independent component analysis into 14 resting state networks (RSN), which were examined for effect of alcohol use and diagnosis-by-alcohol use accounting for age, sex, and diagnosis.

RESULTS

Despite the groups consuming similar amounts of alcohol (BD: mean score ± SD 3.63 ± 3; HC 4.72 ± 3, U = 713, p = .07), for BD participants, greater alcohol use was associated with increased connectivity of the paracingulate gyrus within a default mode network (DMN) and reduced connectivity within an executive control network (ECN) relative to controls. Independently, greater alcohol use was associated with increased connectivity within an ECN and reduced connectivity within a DMN. A diagnosis of BD was associated with increased connectivity of a DMN and reduced connectivity of an ECN.

CONCLUSION

Affective symptomatology in BD is suggested to arise from the aberrant functionality of networks subserving emotive, cognitive, and introspective processes. Taken together, our results suggest that during euthymic periods, alcohol can contribute to the weakening of emotional regulation and response, potentially explaining the increased lability of mood and vulnerability to relapse within the disorder.

Alternations of neuromagnetic activity across neurocognitive core networks among benign childhood epilepsy with centrotemporal spikes: A multi-frequency MEG study

Objective

We aimed to investigate the alternations of neuromagnetic activity across neurocognitive core networks among early untreated children having benign childhood epilepsy with centrotemporal spikes (BECTS).

Methods

We recorded the Magnetoencephalography (MEG) resting-state data from 48 untreated children having BECTS and 24 healthy children. The fourth edition of the Wechsler Intelligence Scale for Children (WISC-IV) was utilized to divide the children with BECTS into two groups: the cognitive impairment (CI) group with a full-scale intelligence quotient (FSIQ) of < 90 and the cognitive non-impairment (CNI) group with an FSIQ of > 90. We selected 26 bilateral cognitive-related regions of interest based on the triple network model. The neurocognitive core network spectral power was estimated using a minimum norm estimate (MNE).

Results

In the CNI group, the spectral power inside the bilateral anterior cingulate cortex (ACC) and the bilateral caudal middle frontal cortex (CMF) enhanced within the delta band and reduced within the alpha band. Both the CI and the CNI group demonstrated enhanced spectral power inside the bilateral posterior cingulate cortex (PCC), bilateral precuneus (PCu) region, bilateral superior and middle temporal cortex, bilateral inferior parietal lobe (IPL), and bilateral supramarginal cortex (SM) region in the delta band. Moreover, there was decreased spectral power in the alpha band. In addition, there were consistent changes in the high-frequency spectrum (> 90 Hz). The spectral power density within the insula cortex (IC), superior temporal cortex (ST), middle temporal cortex (MT), and parahippocampal cortex (PaH) also decreased. Therefore, studying high-frequency activity could lead to a new understanding of the pathogenesis of BECTS.

Conclusion

The alternations of spectral power among neurocognitive core networks could account for CI among early untreated children having BECTS. The dynamic properties of spectral power in different frequency bands could behave as biomarkers for diagnosing new BECTS.

Characterization of resting functional MRI activity alterations across epileptic foci and networks

Brain network alterations have been studied extensively in patients with mesial temporal lobe epilepsy (mTLE) and other focal epilepsies using resting-state functional magnetic resonance imaging (fMRI). However, little has been done to characterize the basic fMRI signal alterations caused by focal epilepsy. Here, we characterize how mTLE affects the fMRI signal in epileptic foci and networks. Resting-state fMRI and diffusion MRI were collected from 47 unilateral mTLE patients and 96 healthy controls. FMRI activity, quantified by amplitude of low-frequency fluctuations, was increased in the epileptic focus and connected regions in mTLE. Evidence for spread of this epileptic fMRI activity was found through linear relationships of regional activity across subjects, the association of these relationships with functional connectivity, and increased activity along white matter tracts. These fMRI activity increases were found to be dependent on the epileptic focus, where the activity was related to disease severity, suggesting the focus to be the origin of these pathological alterations. Furthermore, we found fMRI activity decreases in the default mode network of right mTLE patients with different properties than the activity increases found in the epileptic focus. This work provides insights into basic fMRI signal alterations and their potential spread across networks in focal epilepsy.

Effect of Levodopa Medication on Human Brain Connectome in Parkinson's Disease-A Combined Graph Theory and EEG Study

Background. Levodopa-based drugs are widely used for mitigating the complications induced by Parkinson's disease (PD). Despite the positive effects, several issues regarding the way that levodopa changes brain activities have remained unclear. Methods. A combined strategy using EEG data and graph theory was used for investigating how levodopa changed connectome and processing hubs of the brain during resting-state. Obtained results were subjected to ANOVA test and multiple-comparison post-hoc correction procedure. Results. Outcomes showed that graph topology was not significantly different between PD and healthy groups during the eyes-closed condition, while in the eyes-open condition, statistically significant differences were found. The main effect of levodopa medication was observed for gamma-band activity in which levodopa changed the brain connectome toward a star-like topology. Considering the beta subband of EEG data, graph leaf number increased following levodopa medication in PD patients. Enhanced brain connectivity in the gamma band and reduced beta band connections in the basal ganglia were also observed after levodopa medication. Furthermore, source localization using dipole fitting showed that levodopa suppressed the activity of collateral trigone. Conclusion. Our combined EEG and graph analysis showed that levodopa medication changed the brain connectome, especially in the high-frequency range of brain electrical activities (beta and gamma).

Hemodynamic functional connectivity optimization of frequency EEG microstates enables attention LSTM framework to classify distinct temporal cortical communications of different cognitive tasks

Temporal analysis of global cortical communication of cognitive tasks in coarse EEG information is still challenging due to the underlying complex neural mechanisms. This study proposes an attention-based time-series deep learning framework that processes fMRI functional connectivity optimized quasi-stable frequency microstates for classifying distinct temporal cortical communications of the cognitive task. Seventy volunteers were subjected to visual target detection tasks, and their electroencephalogram (EEG) and functional MRI (fMRI) were acquired simultaneously. At first, the acquired EEG information was preprocessed and bandpass to delta, theta, alpha, beta, and gamma bands and then subjected to quasi-stable frequency-microstate estimation. Subsequently, time-series elicitation of each frequency microstates is optimized with graph theory measures of simultaneously eliciting fMRI functional connectivity between frontal, parietal, and temporal cortices. The distinct neural mechanisms associated with each optimized frequency-microstate were analyzed using microstate-informed fMRI. Finally, these optimized, quasi-stable frequency microstates were employed to train and validate the attention-based Long Short-Term Memory (LSTM) time-series architecture for classifying distinct temporal cortical communications of the target from other cognitive tasks. The temporal, sliding input sampling windows were chosen between 180 to 750 ms/segment based on the stability of transition probabilities of the optimized microstates. The results revealed 12 distinct frequency microstates capable of deciphering target detections' temporal cortical communications from other task engagements. Particularly, fMRI functional connectivity measures of target engagement were observed significantly correlated with the right-diagonal delta (r = 0.31), anterior-posterior theta (r = 0.35), left-right theta (r = - 0.32), alpha (r = - 0.31) microstates. Further, neuro-vascular information of microstate-informed fMRI analysis revealed the association of delta/theta and alpha/beta microstates with cortical communications and local neural processing, respectively. The classification accuracies of the attention-based LSTM were higher than the traditional LSTM architectures, particularly the frameworks that sampled the EEG data with a temporal width of 300 ms/segment. In conclusion, the study demonstrates reliable temporal classifications of global cortical communication of distinct tasks using an attention-based LSTM utilizing fMRI functional connectivity optimized quasi-stable frequency microstates.

External drivers of BOLD signal's non-stationarity

A fundamental challenge in neuroscience is to uncover the principles governing how the brain interacts with the external environment. However, assumptions about external stimuli fundamentally constrain current computational models. We show in silico that unknown external stimulation can produce error in the estimated linear time-invariant dynamical system. To address these limitations, we propose an approach to retrieve the external (unknown) input parameters and demonstrate that the estimated system parameters during external input quiescence uncover spatiotemporal profiles of external inputs over external stimulation periods more accurately. Finally, we unveil the expected (and unexpected) sensory and task-related extra-cortical input profiles using functional magnetic resonance imaging data acquired from 96 subjects (Human Connectome Project) during the resting-state and task scans. This dynamical systems model of the brain offers information on the structure and dimensionality of the BOLD signal's external drivers and shines a light on the likely external sources contributing to the BOLD signal's non-stationarity. Our findings show the role of exogenous inputs in the BOLD dynamics and highlight the importance of accounting for external inputs to unravel the brain's time-varying functional dynamics.

Dynamic spectral signatures of mirror movements in the sensorimotor functional connectivity network of patients with Kallmann syndrome

In Kallmann syndrome (KS), the peculiar phenomenon of bimanual synkinesis or mirror movement (MM) has been associated with a spectral shift, from lower to higher frequencies, of the resting-state fMRI signal of the large-scale sensorimotor brain network (SMN). To possibly determine whether a similar frequency specificity exists across different functional connectivity SMN states, and to capture spontaneous transitions between them, we investigated the dynamic spectral changes of the SMN functional connectivity in KS patients with and without MM symptom. Brain MRI data were acquired at 3 Tesla in 39 KS patients (32 without MM, KSMM-, seven with MM, KSMM+) and 26 age- and sex-matched healthy control (HC) individuals. The imaging protocol included 20-min rs-fMRI scans enabling detailed spectro-temporal analyses of large-scale functional connectivity brain networks. Group independent component analysis was used to extract the SMN. A sliding window approach was used to extract the dynamic spectral power of the SMN functional connectivity within the canonical physiological frequency range of slow rs-fMRI signal fluctuations (0.01–0.25 Hz). K-means clustering was used to determine (and count) the most recurrent dynamic states of the SMN and detect the number of transitions between them. Two most recurrent states were identified, for which the spectral power peaked at a relatively lower (state 1) and higher (state 2) frequency. Compared to KS patients without MM and HC subjects, the SMN of KS patients with MM displayed significantly larger spectral power changes in the slow 3 canonical sub-band (0.073–0.198 Hz) and significantly fewer transitions between state 1 (less recurrent) and state 2 (more recurrent). These findings demonstrate that the presence of MM in KS patients is associated with reduced spontaneous transitions of the SMN between dynamic functional connectivity states and a higher recurrence and an increased spectral power change of the high-frequency state. These results provide novel information about the large-scale brain functional dynamics that could help to understand the pathologic mechanisms of bimanual synkinesis in KS syndrome and, potentially, other neurological disorders where MM may also occur.

Cognitive decline is associated with frequency-specific resting state functional changes in normal aging

Resting state low-frequency brain activity may aid in our understanding of the mechanisms of aging-related cognitive decline. Our purpose was to explore the characteristics of the amplitude of low-frequency fluctuations (ALFF) in different frequency bands of fMRI to better understand cognitive aging. Thirty-seven cognitively normal older individuals underwent a battery of neuropsychological tests and MRI scans at baseline and four years later. ALFF from five different frequency bands (typical band, slow-5, slow-4, slow-3, and slow-2) were calculated and analyzed. A two-way ANOVA was used to explore the interaction effects in voxel-wise whole brain ALFF of the time and frequency bands. Paired-sample t-test was used to explore within-group changes over four years. Partial correlation analysis was performed to assess associations between the altered ALFF and cognitive function. Significant interaction effects of time × frequency were distributed over inferior frontal gyrus, superior frontal gyrus, right rolandic operculum, left thalamus, and right putamen. Significant ALFF reductions in all five frequency bands were mainly found in the right hemisphere and the posterior cerebellum; whereas localization of the significantly increased ALFF were mainly found in the cerebellum at typical band, slow-5 and slow-4 bands, and left hemisphere and the cerebellum at slow-3, slow-2 bands. In addition, ALFF changes showed frequency-specific correlations with changes in cognition. These results suggest that changes of local brain activity in cognitively normal aging should be investigated in multiple frequency bands. The association between ALFF changes and cognitive function can potentially aid better understanding of the mechanisms underlying normal cognitive aging.

Frequency-Specific Analysis of the Dynamic Reconfiguration of the Brain in Patients with Schizophrenia

The analysis of resting-state fMRI signals usually focuses on the low-frequency range/band (0.01−0.1 Hz), which does not cover all aspects of brain activity. Studies have shown that distinct frequency bands can capture unique fluctuations in brain activity, with high-frequency signals (>0.1 Hz) providing valuable information for the diagnosis of schizophrenia. We hypothesized that it is meaningful to study the dynamic reconfiguration of schizophrenia through different frequencies. Therefore, this study used resting-state functional magnetic resonance (RS-fMRI) data from 42 schizophrenia and 40 normal controls to investigate dynamic network reconfiguration in multiple frequency bands (0.01−0.25 Hz, 0.01−0.027 Hz, 0.027−0.073 Hz, 0.073−0.198 Hz, 0.198−0.25 Hz). Based on the time-varying dynamic network constructed for each frequency band, we compared the dynamic reconfiguration of schizophrenia and normal controls by calculating the recruitment and integration. The experimental results showed that the differences between schizophrenia and normal controls are observed in the full frequency, which is more significant in slow3. In addition, as visual network, attention network, and default mode network differ a lot from each other, they can show a high degree of connectivity, which indicates that the functional network of schizophrenia is affected by the abnormal brain state in these areas. These shreds of evidence provide a new perspective and promote the current understanding of the characteristics of dynamic brain networks in schizophrenia.

Frequency-specific coactivation patterns in resting-state and their alterations in schizophrenia: An fMRI study

The resting‐state human brain is a dynamic system that shows frequency‐dependent characteristics. Recent studies demonstrate that coactivation pattern (CAP) analysis can identify recurring brain states with similar coactivation configurations. However, it is unclear whether and how CAPs depend on the frequency bands. The current study investigated the spatial and temporal characteristics of CAPs in the four frequency sub‐bands from slow‐5 (0.01–0.027 Hz), slow‐4 (0.027–0.073 Hz), slow‐3 (0.073–0.198 Hz), to slow‐2 (0.198–0.25 Hz), in addition to the typical low‐frequency range (0.01–0.08 Hz). In the healthy subjects, six CAP states were obtained at each frequency band in line with our prior study. Similar spatial patterns with the typical range were observed in slow‐5, 4, and 3, but not in slow‐2. While the frequency increased, all CAP states displayed shorter persistence, which caused more between‐state transitions. Specifically, from slow‐5 to slow‐4, the coactivation not only changed significantly in distributed cortical networks, but also increased in the basal ganglia as well as the amygdala. Schizophrenia patients showed significant alteration in the persistence of CAPs of slow‐5. Using leave‐one‐pair‐out, hold‐out and resampling validations, the highest classification accuracy (84%) was achieved by slow‐4 among different frequency bands. In conclusion, our findings provide novel information about spatial and temporal characteristics of CAP states at different frequency bands, which contributes to a better understanding of the frequency aspect of biomarkers for schizophrenia and other disorders.

Alteration in Functional Magnetic Resonance Imaging Signal Complexity Across Multiple Time Scales in Patients With Migraine Without Aura

Background

Migraine is a primary neurological disorder associated with complex brain activity. Recently, mounting evidence has suggested that migraine is underpinned by aberrant dynamic brain activity characterized by linear and non-linear changes across a variety of time scales. However, the abnormal dynamic brain activity at different time scales is still unknown in patients with migraine without aura (MWoA). This study aimed to assess the altered patterns of brain activity dynamics over different time scales and the potential pathophysiological mechanisms of alterations in patients with MWoA.

Methods

Multiscale entropy in 50 patients and 20 healthy controls (HCs) was calculated to investigate the patterns and altered brain complexity (BC) across five different time scales. Spearman rank correlation analysis between BC in regions showing significant intergroup differences and clinical scores (i.e., frequency of migraine attacks, duration, headache impact test) was conducted in patients with MWoA.

Results

The spatial distribution of BC varied across different time scales. At time scale1, BC was higher in the posterior default mode network (DMN) across participants. Compared with HCs, patients with MWoA had higher BC in the DMN and sensorimotor network. At time scale2, BC was mainly higher in the anterior DMN across participants. Patients with MWoA had higher BC in the sensorimotor network. At time scale3, BC was mainly higher in the frontoparietal network across participants. Patients with MWoA had increased BC in the parietal gyrus. At time scale4, BC is mainly higher in the sensorimotor network. Patients with MWoA had higher BC in the postcentral gyrus. At time scale5, BC was mainly higher in the DMN. Patients with MWoA had lower BC in the posterior DMN. In particular, BC values in the precuneus and paracentral lobule significantly correlated with clinical symptoms.

Conclusion

Migraine is associated with alterations in dynamic brain activity in the sensorimotor network and DMN over multiple time scales. Time-varying BC within these regions could be linked to instability in pain transmission and modulation. Our findings provide new evidence for the hypothesis of abnormal dynamic brain activity in migraine.

Brain Imaging Biomarkers for Chronic Pain

The prevalence of chronic pain has reached epidemic levels. In addition to personal suffering chronic pain is associated with psychiatric and medical co-morbidities, notably substance misuse, and a huge a societal cost amounting to hundreds of billions of dollars annually in medical cost, lost wages, and productivity. Chronic pain does not have a cure or quantitative diagnostic or prognostic tools. In this manuscript we provide evidence that this situation is about to change. We first start by summarizing our current understanding of the role of the brain in the pathogenesis of chronic pain. We particularly focus on the concept of learning in the emergence of chronic pain, and the implication of the limbic brain circuitry and dopaminergic signaling, which underly emotional learning and decision making, in this process. Next, we summarize data from our labs and from other groups on the latest brain imaging findings in different chronic pain conditions focusing on results with significant potential for translation into clinical applications. The gaps in the study of chronic pain and brain imaging are highlighted in throughout the overview. Finally, we conclude by discussing the costs and benefits of using brain biomarkers of chronic pain and compare to other potential markers.

Predicting behavior through dynamic modes in resting-state fMRI data

Dynamic properties of resting-state functional connectivity (FC) provide rich information on brain-behavior relationships. Dynamic mode decomposition (DMD) has been used as a method to characterize FC dynamics. However, it remains unclear whether dynamic modes (DMs), spatial-temporal coherent patterns computed by DMD, provide information about individual behavioral differences. This study established a methodological approach to predict individual differences in behavior using DMs. Furthermore, we investigated the contribution of DMs within each of seven specific frequency bands (0-0.1,...,0.6-0.7 Hz) for prediction. To validate our approach, we tested whether each of 59 behavioral measures could be predicted by performing multivariate pattern analysis on a Gram matrix, which was created using subject-specific DMs computed from resting-state functional magnetic resonance imaging (rs-fMRI) data of individuals. DMD successfully predicted behavior and outperformed temporal and spatial independent component analysis, which is the conventional data decomposition method for extracting spatial activity patterns. Most of the behavioral measures that were predicted with significant accuracy in a permutation test were related to cognition. We found that DMs within frequency bands <0.2 Hz primarily contributed to prediction and had spatial structures similar to several common resting-state networks. Our results indicate that DMD is efficient in extracting spatiotemporal features from rs-fMRI data.

Kuramoto Model-Based Analysis Reveals Oxytocin Effects on Brain Network Dynamics

The oxytocin effects on large-scale brain networks such as Default Mode Network (DMN) and Frontoparietal Network (FPN) have been largely studied using fMRI data. However, these studies are mainly based on the statistical correlation or Bayesian causality inference, lacking interpretability at the physical and neuroscience level. Here, we propose a physics-based framework of the Kuramoto model to investigate oxytocin effects on the phase dynamic neural coupling in DMN and FPN. Testing on fMRI data of 59 participants administrated with either oxytocin or placebo, we demonstrate that oxytocin changes the topology of brain communities in DMN and FPN, leading to higher synchronization in the FPN and lower synchronization in the DMN, as well as a higher variance of the coupling strength within the DMN and more flexible coupling patterns at group level. These results together indicate that oxytocin may increase the ability to overcome the corresponding internal oscillation dispersion and support the flexibility in neural synchrony in various social contexts, providing new evidence for explaining the oxytocin modulated social behaviors. Our proposed Kuramoto model-based framework can be a potential tool in network neuroscience and offers physical and neural insights into phase dynamics of the brain.

Multi-Band Brain Network Analysis for Functional Neuroimaging Biomarker Identification

The functional connectomic profile is one of the non-invasive imaging biomarkers in the computer-assisted diagnostic system for many neuro-diseases. However, the diagnostic power of functional connectivity is challenged by mixed frequency-specific neuronal oscillations in the brain, which makes the single Functional Connectivity Network (FCN) often underpowered to capture the disease-related functional patterns. To address this challenge, we propose a novel functional connectivity analysis framework to conduct joint feature learning and personalized disease diagnosis, in a semi-supervised manner, aiming at focusing on putative multi-band functional connectivity biomarkers from functional neuroimaging data. Specifically, we first decompose the Blood Oxygenation Level Dependent (BOLD) signals into multiple frequency bands by the discrete wavelet transform, and then cast the alignment of all fully-connected FCNs derived from multiple frequency bands into a parameter-free multi-band fusion model. The proposed fusion model fuses all fully-connected FCNs to obtain a sparsely-connected FCN (sparse FCN for short) for each individual subject, as well as lets each sparse FCN be close to its neighbored sparse FCNs and be far away from its furthest sparse FCNs. Furthermore, we employ the l1 -SVM to conduct joint brain region selection and disease diagnosis. Finally, we evaluate the effectiveness of our proposed framework on various neuro-diseases, i.e., Fronto-Temporal Dementia (FTD), Obsessive-Compulsive Disorder (OCD), and Alzheimer's Disease (AD), and the experimental results demonstrate that our framework shows more reasonable results, compared to state-of-the-art methods, in terms of classification performance and the selected brain regions. The source code can be visited by the url https://github.com/reynard-hu/mbbna.

Discrete spectral eigenmode-resonance network of brain dynamics and connectivity

The problem of finding a compact natural representation of brain dynamics and connectivity is addressed using an expansion in terms of physical spatial eigenmodes and their frequency resonances. It is demonstrated that this discrete expansion via the system transfer function enables linear and nonlinear dynamics to be analyzed in compact form in terms of natural dynamic "atoms," each of which is a frequency resonance of an eigenmode. Because these modal resonances are determined by the system dynamics, not the investigator, they are privileged over widely used phenomenological patterns, and obviate the need for artificial discretizations and thresholding in coordinate space. It is shown that modal resonances participate as nodes of a discrete spectral network, are noninteracting in the linear regime, but are linked nonlinearly by wave-wave coalescence and decay processes. The modal resonance formulation is shown to be capable of speeding numerical calculations of strongly nonlinear interactions. Recent work in brain dynamics, especially based on neural field theory (NFT) approaches, allows eigenmodes and their resonances to be estimated from data without assuming a specific brain model. This means that dynamic equations can be inferred using system identification methods from control theory, rather than being assumed, and resonances can be interpreted as control-systems data filters. The results link brain activity and connectivity with control-systems functions such as prediction and attention via gain control and can also be linked to specific NFT predictions if desired, thereby providing a convenient bridge between physiologically based theories and experiment. Amplitudes of modes and resonances can also be tracked to provide a more direct and temporally localized representation of the dynamics than correlations and covariances, which are widely used in the field. By synthesizing many different lines of research, this work provides a way to link quantitative electrophysiological and imaging measurements, connectivity, brain dynamics, and function. This underlines the need to move between coordinate and spectral representations as required. Moreover, standard theoretical-physics approaches and mathematical methods can be used in place of ad hoc statistical measures such as those based on graph theory of artificially discretized and decimated networks, which are highly prone to selection effects and artifacts.