Brain charts for the human lifespan.

Over the past few decades, neuroimaging has become a ubiquitous tool in basic research and clinical studies of the human brain. However, no reference standards currently exist to quantify individual differences in neuroimaging metrics over time, in contrast to growth charts for anthropometric traits such as height and weight1. Here we assemble an interactive open resource to benchmark brain morphology derived from any current or future sample of MRI data ( http://www.brainchart.io/ ). With the goal of basing these reference charts on the largest and most inclusive dataset available, acknowledging limitations due to known biases of MRI studies relative to the diversity of the global population, we aggregated 123,984 MRI scans, across more than 100 primary studies, from 101,457 human participants between 115 days post-conception to 100 years of age. MRI metrics were quantified by centile scores, relative to non-linear trajectories2 of brain structural changes, and rates of change, over the lifespan. Brain charts identified previously unreported neurodevelopmental milestones3, showed high stability of individuals across longitudinal assessments, and demonstrated robustness to technical and methodological differences between primary studies. Centile scores showed increased heritability compared with non-centiled MRI phenotypes, and provided a standardized measure of atypical brain structure that revealed patterns of neuroanatomical variation across neurological and psychiatric disorders. In summary, brain charts are an essential step towards robust quantification of individual variation benchmarked to normative trajectories in multiple, commonly used neuroimaging phenotypes.

Identifying canonical and replicable multi-scale intrinsic connectivity networks in 100k+ resting-state fMRI datasets.

Despite the known benefits of data-driven approaches, the lack of approaches for identifying functional neuroimaging patterns that capture both individual variations and inter-subject correspondence limits the clinical utility of rsfMRI and its application to single-subject analyses. Here, using rsfMRI data from over 100k individuals across private and public datasets, we identify replicable multi-spatial-scale canonical intrinsic connectivity network (ICN) templates via the use of multi-model-order independent component analysis (ICA). We also study the feasibility of estimating subject-specific ICNs via spatially constrained ICA. The results show that the subject-level ICN estimations vary as a function of the ICN itself, the data length, and the spatial resolution. In general, large-scale ICNs require less data to achieve specific levels of (within- and between-subject) spatial similarity with their templates. Importantly, increasing data length can reduce an ICN's subject-level specificity, suggesting longer scans may not always be desirable. We also find a positive linear relationship between data length and spatial smoothness (possibly due to averaging over intrinsic dynamics), suggesting studies examining optimized data length should consider spatial smoothness. Finally, consistency in spatial similarity between ICNs estimated using the full data and subsets across different data lengths suggests lower within-subject spatial similarity in shorter data is not wholly defined by lower reliability in ICN estimates, but may be an indication of meaningful brain dynamics which average out as data length increases.

Moving beyond the 'CAP' of the Iceberg: Intrinsic connectivity networks in fMRI are continuously engaging and overlapping.

Resting-state functional magnetic resonance imaging is currently the mainstay of functional neuroimaging and has allowed researchers to identify intrinsic connectivity networks (aka functional networks) at different spatial scales. However, little is known about the temporal profiles of these networks and whether it is best to model them as continuous phenomena in both space and time or, rather, as a set of temporally discrete events. Both categories have been supported by series of studies with promising findings. However, a critical question is whether focusing only on time points presumed to contain isolated neural events and disregarding the rest of the data is missing important information, potentially leading to misleading conclusions. In this work, we argue that brain networks identified within the spontaneous blood oxygenation level-dependent (BOLD) signal are not limited to temporally sparse burst moments and that these event present time points (EPTs) contain valuable but incomplete information about the underlying functional patterns. We focus on the default mode and show evidence that is consistent with its continuous presence in the BOLD signal, including during the event absent time points (EATs), i.e., time points that exhibit minimum activity and are the least likely to contain an event. Moreover, our findings suggest that EPTs may not contain all the available information about their corresponding networks. We observe distinct default mode connectivity patterns obtained from all time points (AllTPs), EPTs, and EATs. We show evidence of robust relationships with schizophrenia symptoms that are both common and unique to each of the sets of time points (AllTPs, EPTs, EATs), likely related to transient patterns of connectivity. Together, these findings indicate the importance of leveraging the full temporal data in functional studies, including those using event-detection approaches.

The Developmental Chronnecto-Genomics (Dev-CoG) study: A multimodal study on the developing brain.

Brain development has largely been studied through unimodal analysis of neuroimaging data, providing independent results for structural and functional data. However, structure clearly impacts function and vice versa, pointing to the need for performing multimodal data collection and analysis to improve our understanding of brain development, and to further inform models of typical and atypical brain development across the lifespan. Ultimately, such models should also incorporate genetic and epigenetic mechanisms underlying brain structure and function, although currently this area is poorly specified. To this end, we are reporting here a multi-site, multi-modal dataset that captures cognitive function, brain structure and function, and genetic and epigenetic measures to better quantify the factors that influence brain development in children originally aged 9-14 years. Data collection for the Developmental Chronnecto-Genomics (Dev-CoG) study (http://devcog.mrn.org/) includes cognitive, emotional, and social performance scales, structural and functional MRI, diffusion MRI, magnetoencephalography (MEG), and saliva collection for DNA analysis of single nucleotide polymorphisms (SNPs) and DNA methylation patterns. Across two sites (The Mind Research Network and the University of Nebraska Medical Center), data from over 200 participants were collected and these children were re-tested annually for at least 3 years. The data collection protocol, sample demographics, and data quality measures for the dataset are presented here. The sample will be made freely available through the collaborative informatics and neuroimaging suite (COINS) database at the conclusion of the study.

Transient increased thalamic-sensory connectivity and decreased whole-brain dynamism in autism.

Autism spectrum disorder (ASD) is a neurodevelopmental disorder associated with social communication deficits and restricted/repetitive behaviors and is characterized by large-scale atypical subcortical-cortical connectivity, including impaired resting-state functional connectivity between thalamic and sensory regions. Previous studies have typically focused on the abnormal static connectivity in ASD and overlooked potential valuable dynamic patterns in brain connectivity. However, resting-state brain connectivity is indeed highly dynamic, and abnormalities in dynamic brain connectivity have been widely identified in psychiatric disorders. In this study, we investigated the dynamic functional network connectivity (dFNC) between 51 intrinsic connectivity networks in 170 individuals with ASD and 195 age-matched typically developing (TD) controls using independent component analysis and a sliding window approach. A hard clustering state analysis and a fuzzy meta-state analysis were conducted respectively, for the exploration of local and global aberrant dynamic connectivity patterns in ASD. We examined the group difference in dFNC between thalamic and sensory networks in each functional state and group differences in four high-dimensional dynamic measures. The results showed that compared with TD controls, individuals with ASD show an increase in transient connectivity between hypothalamus/subthalamus and some sensory networks (right postcentral gyrus, bi paracentral lobule, and lingual gyrus) in certain functional states, and diminished global meta-state dynamics of the whole-brain functional network. In addition, these atypical dynamic patterns are significantly associated with autistic symptoms indexed by the Autism Diagnostic Observation Schedule. These converging results support and extend previous observations regarding hyperconnectivity between thalamic and sensory regions and stable whole-brain functional configuration in ASD. Dynamic brain connectivity may serve as a potential biomarker of ASD and further investigation of these dynamic patterns might help to advance our understanding of behavioral differences in this complex neurodevelopmental disorder.

Characterizing dynamic amplitude of low-frequency fluctuation and its relationship with dynamic functional connectivity: An application to schizophrenia.

The human brain is a highly dynamic system with non-stationary neural activity and rapidly-changing neural interaction. Resting-state dynamic functional connectivity (dFC) has been widely studied during recent years, and the emerging aberrant dFC patterns have been identified as important features of many mental disorders such as schizophrenia (SZ). However, only focusing on the time-varying patterns in FC is not enough, since the local neural activity itself (in contrast to the inter-connectivity) is also found to be highly fluctuating from research using high-temporal-resolution imaging techniques. Exploring the time-varying patterns in brain activity and their relationships with time-varying brain connectivity is important for advancing our understanding of the co-evolutionary property of brain network and the underlying mechanism of brain dynamics. In this study, we introduced a framework for characterizing time-varying brain activity and exploring its associations with time-varying brain connectivity, and applied this framework to a resting-state fMRI dataset including 151 SZ patients and 163 age- and gender matched healthy controls (HCs). In this framework, 48 brain regions were first identified as intrinsic connectivity networks (ICNs) using group independent component analysis (GICA). A sliding window approach was then adopted for the estimation of dynamic amplitude of low-frequency fluctuation (dALFF) and dFC, which were used to measure time-varying brain activity and time-varying brain connectivity respectively. The dALFF was further clustered into six reoccurring states by the k-means clustering method and the group difference in occurrences of dALFF states was explored. Lastly, correlation coefficients between dALFF and dFC were calculated and the group difference in these dALFF-dFC correlations was explored. Our results suggested that 1) ALFF of brain regions was highly fluctuating during the resting-state and such dynamic patterns are altered in SZ, 2) dALFF and dFC were correlated in time and their correlations are altered in SZ. The overall results support and expand prior work on abnormalities of brain activity, static FC (sFC) and dFC in SZ, and provide new evidence on aberrant time-varying brain activity and its associations with brain connectivity in SZ, which might underscore the disrupted brain cognitive functions in this mental disorder.

Resting-state thalamic dysconnectivity in schizophrenia and relationships with symptoms.

BACKGROUND: Schizophrenia (SZ) is a severe neuropsychiatric disorder associated with disrupted connectivity within the thalamic-cortico-cerebellar network. Resting-state functional connectivity studies have reported thalamic hypoconnectivity with the cerebellum and prefrontal cortex as well as thalamic hyperconnectivity with sensory cortical regions in SZ patients compared with healthy comparison participants (HCs). However, fundamental questions remain regarding the clinical significance of these connectivity abnormalities. METHOD: Resting state seed-based functional connectivity was used to investigate thalamus to whole brain connectivity using multi-site data including 183 SZ patients and 178 matched HCs. Statistical significance was based on a voxel-level FWE-corrected height threshold of p < 0.001. The relationships between positive and negative symptoms of SZ and regions of the brain demonstrating group differences in thalamic connectivity were examined. RESULTS: HC and SZ participants both demonstrated widespread positive connectivity between the thalamus and cortical regions. Compared with HCs, SZ patients had reduced thalamic connectivity with bilateral cerebellum and anterior cingulate cortex. In contrast, SZ patients had greater thalamic connectivity with multiple sensory-motor regions, including bilateral pre- and post-central gyrus, middle/inferior occipital gyrus, and middle/superior temporal gyrus. Thalamus to middle temporal gyrus connectivity was positively correlated with hallucinations and delusions, while thalamus to cerebellar connectivity was negatively correlated with delusions and bizarre behavior. CONCLUSIONS: Thalamic hyperconnectivity with sensory regions and hypoconnectivity with cerebellar regions in combination with their relationship to clinical features of SZ suggest that thalamic dysconnectivity may be a core neurobiological feature of SZ that underpins positive symptoms.

EEG Signatures of Dynamic Functional Network Connectivity States.

The human brain operates by dynamically modulating different neural populations to enable goal directed behavior. The synchrony or lack thereof between different brain regions is thought to correspond to observed functional connectivity dynamics in resting state brain imaging data. In a large sample of healthy human adult subjects and utilizing a sliding windowed correlation method on functional imaging data, earlier we demonstrated the presence of seven distinct functional connectivity states/patterns between different brain networks that reliably occur across time and subjects. Whether these connectivity states correspond to meaningful electrophysiological signatures was not clear. In this study, using a dataset with concurrent EEG and resting state functional imaging data acquired during eyes open and eyes closed states, we demonstrate the replicability of previous findings in an independent sample, and identify EEG spectral signatures associated with these functional network connectivity changes. Eyes open and eyes closed conditions show common and different connectivity patterns that are associated with distinct EEG spectral signatures. Certain connectivity states are more prevalent in the eyes open case and some occur only in eyes closed state. Both conditions exhibit a state of increased thalamocortical anticorrelation associated with reduced EEG spectral alpha power and increased delta and theta power possibly reflecting drowsiness. This state occurs more frequently in the eyes closed state. In summary, we find a link between dynamic connectivity in fMRI data and concurrently collected EEG data, including a large effect of vigilance on functional connectivity. As demonstrated with EEG and fMRI, the stationarity of connectivity cannot be assumed, even for relatively short periods.

Subcortical brain volume abnormalities in 2028 individuals with schizophrenia and 2540 healthy controls via the ENIGMA consortium.

The profile of brain structural abnormalities in schizophrenia is still not fully understood, despite decades of research using brain scans. To validate a prospective meta-analysis approach to analyzing multicenter neuroimaging data, we analyzed brain MRI scans from 2028 schizophrenia patients and 2540 healthy controls, assessed with standardized methods at 15 centers worldwide. We identified subcortical brain volumes that differentiated patients from controls, and ranked them according to their effect sizes. Compared with healthy controls, patients with schizophrenia had smaller hippocampus (Cohen's d=-0.46), amygdala (d=-0.31), thalamus (d=-0.31), accumbens (d=-0.25) and intracranial volumes (d=-0.12), as well as larger pallidum (d=0.21) and lateral ventricle volumes (d=0.37). Putamen and pallidum volume augmentations were positively associated with duration of illness and hippocampal deficits scaled with the proportion of unmedicated patients. Worldwide cooperative analyses of brain imaging data support a profile of subcortical abnormalities in schizophrenia, which is consistent with that based on traditional meta-analytic approaches. This first ENIGMA Schizophrenia Working Group study validates that collaborative data analyses can readily be used across brain phenotypes and disorders and encourages analysis and data sharing efforts to further our understanding of severe mental illness.

Lateralization of resting state networks and relationship to age and gender.

Brain lateralization is a widely studied topic, however there has been little work focused on lateralization of intrinsic networks (regions showing similar patterns of covariation among voxels) in the resting brain. In this study, we evaluate resting state network lateralization in an age and gender-balanced functional magnetic resonance imaging (fMRI) dataset comprising over 600 healthy subjects ranging in age from 12 to 71. After establishing sample-wide network lateralization properties, we continue with an investigation of age and gender effects on network lateralization. All data was gathered on the same scanner and preprocessed using an automated pipeline (Scott et al., 2011). Networks were extracted via group independent component analysis (gICA) (Calhoun et al., 2001). Twenty-eight resting state networks discussed in previous (Allen et al., 2011) work were re-analyzed with a focus on lateralization. We calculated homotopic voxelwise measures of laterality in addition to a global lateralization measure, called the laterality cofactor, for each network. As expected, many of the intrinsic brain networks were lateralized. For example, the visual network was strongly right lateralized, auditory network and default mode networks were mostly left lateralized. Attentional and frontal networks included nodes that were left lateralized and other nodes that were right lateralized. Age was strongly related to lateralization in multiple regions including sensorimotor network regions precentral gyrus, postcentral gyrus and supramarginal gyrus; and visual network regions lingual gyrus; attentional network regions inferior parietal lobule, superior parietal lobule and middle temporal gyrus; and frontal network regions including the inferior frontal gyrus. Gender showed significant effects mainly in two regions, including visual and frontal networks. For example, the inferior frontal gyrus was more right lateralized in males. Significant effects of age were found in sensorimotor and visual networks on the global measure. In summary, we report a large-sample of lateralization study that finds intrinsic functional brain networks to be highly lateralized, with regions that are strongly related to gender and age locally, and with age a strong factor in lateralization, and gender exhibiting a trend-level effect on global measures of laterality.

Distinct intrinsic network connectivity patterns of post-traumatic stress disorder symptom clusters.

OBJECTIVE: Post-traumatic stress disorder (PTSD) is considered a multidimensional disorder, with distinct symptom clusters including re-experiencing, avoidance/numbing, hyperarousal, and most recently depersonalization/derealization. However, the extent of differing intrinsic network connectivity underlying these symptoms has not been fully investigated. We therefore investigated the degree of association between resting connectivity of the salience (SN), default mode (DMN), and central executive (CEN) networks and PTSD symptom severity. METHOD: Using resting-state functional MRI data from PTSD participants (n = 21), we conducted multivariate analyses to test whether connectivity of extracted independent components varied as a function of re-experiencing, avoidance/numbing, hyperarousal, and depersonalization/derealization. RESULTS: Hyperarousal symptoms were associated with reduced connectivity of posterior insula/superior temporal gyrus within SN [peak Montréal Neurological Institute (MNI): -44, -8, 0, t = -4.2512, k = 40]. Depersonalization/derealization severity was associated with decreased connectivity of perigenual anterior cingulate/ventromedial prefrontal cortex within ventral anterior DMN (peak MNI: 8, 40, -4; t = -3.8501; k = 15) and altered synchrony between two DMN components and between DMN and CEN. CONCLUSION: Our results are consistent with prior research showing intrinsic network disruptions in PTSD and imply heterogeneous connectivity patterns underlying PTSD symptom dimensions. These findings suggest possible biomarkers for PTSD and its dissociative subtype.

Functional connectivity in the developing brain: a longitudinal study from 4 to 9months of age.

We characterize the development of intrinsic connectivity networks (ICNs) from 4 to 9months of age with resting state magnetic resonance imaging performed on sleeping infants without sedative medication. Data is analyzed with independent component analysis (ICA). Using both low (30 components) and high (100 components) ICA model order decompositions, we find that the functional network connectivity (FNC) map is largely similar at both 4 and 9months. However at 9months the connectivity strength decreases within local networks and increases between more distant networks. The connectivity within the default-mode network, which contains both local and more distant nodes, also increases in strength with age. The low frequency power spectrum increases with age only in the posterior cingulate cortex and posterior default mode network. These findings are consistent with a general developmental pattern of increasing longer distance functional connectivity over the first year of life and raise questions regarding the developmental importance of the posterior cingulate at this age.

The impact of parent socio-economic status on executive functioning and cortical morphology in individuals with schizophrenia and healthy controls.

BACKGROUND: Relatively lower executive functioning is characteristic of individuals with schizophrenia. As low socio-economic status (SES) early in life (i.e. parent SES) has been linked with lower executive skills in healthy children, we hypothesized that parental SES (pSES) would be more strongly related to executive functioning in individuals with schizophrenia than in controls and have a greater impact on prefrontal cortical morphology. METHOD: Healthy controls (n = 125) and individuals with schizophrenia (n = 102) completed tests assessing executive functioning and intelligence. The groups were matched on pSES, which was evaluated with the Hollingshead-Redlich scale. A principal components analysis (PCA) was conducted on 10 variables from six executive tests, yielding three specific components (fluency, planning and response inhibition). Voxel-based morphometry (VBM) was used to evaluate effects of pSES on gray matter (GM) concentration. RESULTS: Lower pSES was associated with lower scores across the three executive functioning components, and a significant group by pSES interaction was observed such that low pSES, in particular, affected individuals with schizophrenia. These effects remained significant when intellectual ability, education and self-SES (sSES) were added as covariates. VBM revealed that lower pSES was associated with reduced GM volume in several anterior brain regions, especially the superior frontal gyrus, in patients but not in controls. CONCLUSIONS: These findings suggest that individuals with schizophrenia may be particularly vulnerable to the adverse impact of low pSES, in terms of both lower executive skills and reduced anterior GM volumes.

Plastic modulation of PTSD resting-state networks and subjective wellbeing by EEG neurofeedback.

OBJECTIVE: Electroencephalographic (EEG) neurofeedback training has been shown to produce plastic modulations in salience network and default mode network functional connectivity in healthy individuals. In this study, we investigated whether a single session of neurofeedback training aimed at the voluntary reduction of alpha rhythm (8-12 Hz) amplitude would be related to differences in EEG network oscillations, functional MRI (fMRI) connectivity, and subjective measures of state anxiety and arousal in a group of individuals with post-traumatic stress disorder (PTSD). METHOD: Twenty-one individuals with PTSD related to childhood abuse underwent 30 min of EEG neurofeedback training preceded and followed by a resting-state fMRI scan. RESULTS: Alpha desynchronizing neurofeedback was associated with decreased alpha amplitude during training, followed by a significant increase ('rebound') in resting-state alpha synchronization. This rebound was linked to increased calmness, greater salience network connectivity with the right insula, and enhanced default mode network connectivity with bilateral posterior cingulate, right middle frontal gyrus, and left medial prefrontal cortex. CONCLUSION: Our study represents a first step in elucidating the potential neurobehavioural mechanisms mediating the effects of neurofeedback treatment on regulatory systems in PTSD. Moreover, it documents for the first time a spontaneous EEG 'rebound' after neurofeedback, pointing to homeostatic/compensatory mechanisms operating in the brain.

MIR137 polygenic risk for schizophrenia and ephrin-regulated pathway: Role in lateral ventricles and corpus callosum volume.

Background/Objective. Enlarged lateral ventricle (LV) volume and decreased volume in the corpus callosum (CC) are hallmarks of schizophrenia (SZ). We previously showed an inverse correlation between LV and CC volumes in SZ, with global functioning decreasing with increased LV volume. This study investigates the relationship between LV volume, CC abnormalities, and the microRNA MIR137 and its regulated genes in SZ, because of MIR137's essential role in neurodevelopment. Methods. Participants were 1224 SZ probands and 1466 unaffected controls from the GENUS Consortium. Brain MRI scans, genotype, and clinical data were harmonized across cohorts and employed in the analyses. Results. Increased LV volumes and decreased CC central, mid-anterior, and mid-posterior volumes were observed in SZ probands. The MIR137-regulated ephrin pathway was significantly associated with CC:LV ratio, explaining a significant proportion (3.42 %) of CC:LV variance, and more than for LV and CC separately. Other pathways explained variance in either CC or LV, but not both. CC:LV ratio was also positively correlated with Global Assessment of Functioning, supporting previous subsample findings. SNP-based heritability estimates were higher for CC central:LV ratio (0.79) compared to CC or LV separately. Discussion. Our results indicate that the CC:LV ratio is highly heritable, influenced in part by variation in the MIR137-regulated ephrin pathway. Findings suggest that the CC:LV ratio may be a risk indicator in SZ that correlates with global functioning.

Using joint ICA to link function and structure using MEG and DTI in schizophrenia.

In this study we employed joint independent component analysis (jICA) to perform a novel multivariate integration of magnetoencephalography (MEG) and diffusion tensor imaging (DTI) data to investigate the link between function and structure. This model-free approach allows one to identify covariation across modalities with different temporal and spatial scales [temporal variation in MEG and spatial variation in fractional anisotropy (FA) maps]. Healthy controls (HC) and patients with schizophrenia (SP) participated in an auditory/visual multisensory integration paradigm to probe cortical connectivity in schizophrenia. To allow direct comparisons across participants and groups, the MEG data were registered to an average head position and regional waveforms were obtained by calculating the local field power of the planar gradiometers. Diffusion tensor images obtained in the same individuals were preprocessed to provide FA maps for each participant. The MEG/FA data were then integrated using the jICA software (http://mialab.mrn.org/software/fit). We identified MEG/FA components that demonstrated significantly different (p<0.05) covariation in MEG/FA data between diagnostic groups (SP vs. HC) and three components that captured the predominant sensory responses in the MEG data. Lower FA values in bilateral posterior parietal regions, which include anterior/posterior association tracts, were associated with reduced MEG amplitude (120-170 ms) of the visual response in occipital sensors in SP relative to HC. Additionally, increased FA in a right medial frontal region was linked with larger amplitude late MEG activity (300-400 ms) in bilateral central channels for SP relative to HC. Step-wise linear regression provided evidence that right temporal, occipital and late central components were significant predictors of reaction time and cognitive performance based on the Measurement and Treatment Research to Improve Cognition in Schizophrenia (MATRICS) cognitive assessment battery. These results point to dysfunction in a posterior visual processing network in schizophrenia, with reduced MEG amplitude, reduced FA and poorer overall performance on the MATRICS. Interestingly, the spatial location of the MEG activity and the associated FA regions are spatially consistent with white matter regions that subserve these brain areas. This novel approach provides evidence for significant pairing between function (neurophysiology) and structure (white matter integrity) and demonstrates that this multivariate, multimodal integration technique is sensitive to group differences in function and structure.

A method for estimating and characterizing explicitly nonlinear dynamic functional network connectivity in resting-state fMRI data.

The past 10 years have seen an explosion of approaches that focus on the study of time-resolved change in functional connectivity (FC). FC characterization among networks at a whole-brain level is frequently termed functional network connectivity (FNC). Time-resolved or dynamic functional network connectivity (dFNC) focuses on the estimation of transient, recurring, whole-brain patterns of FNC. While most approaches in this area have attempted to capture dynamic linear correlation, we are particularly interested in whether explicitly nonlinear relationships, above and beyond linear, are present and contain unique information. This study thus proposes an approach to assess explicitly nonlinear dynamic functional network connectivity (EN dFNC) derived from the relationship among independent component analysis time courses. Linear relationships were removed at each time point to evaluate, typically ignored, explicitly nonlinear dFNC using normalized mutual information (NMI). Simulations showed the proposed method estimated explicitly nonlinearity over time, even within relatively short windows of data. We then, applied our approach on 151 schizophrenia patients, and 163 healthy controls fMRI data and found three unique, highly structured, mostly long-range, functional states that also showed significant group differences. In particular, explicitly nonlinear relationships tend to be more widespread than linear ones. Results also highlighted a state with long range connections to the visual domain, which were significantly reduced in schizophrenia. Overall, this work suggests that quantifying EN dFNC may provide a complementary and potentially valuable tool for studying brain function by exposing relevant variation that is typically ignored.

Spatial Dynamic Subspaces Encode Sex-Specific Schizophrenia Disruptions in Transient Network Overlap and its Links to Genetic Risk.

Background: Recent advances in resting-state fMRI allow us to study spatial dynamics, the phenomenon of brain networks spatially evolving over time. However, most dynamic studies still use subject-specific, spatially-static nodes. As recent studies have demonstrated, incorporating time-resolved spatial properties is crucial for precise functional connectivity estimation and gaining unique insights into brain function. Nevertheless, estimating time-resolved networks poses challenges due to the low signal-to-noise ratio, limited information in short time segments, and uncertain identification of corresponding networks within and between subjects. Methods: We adapt a reference-informed network estimation technique to capture time-resolved spatial networks and their dynamic spatial integration and segregation. We focus on time-resolved spatial functional network connectivity (spFNC), an estimate of network spatial coupling, to study sex-specific alterations in schizophrenia and their links to multi-factorial genomic data. Results: Our findings are consistent with the dysconnectivity and neurodevelopment hypotheses and align with the cerebello-thalamo-cortical, triple-network, and frontoparietal dysconnectivity models, helping to unify them. The potential unification offers a new understanding of the underlying mechanisms. Notably, the posterior default mode/salience spFNC exhibits sex-specific schizophrenia alteration during the state with the highest global network integration and correlates with genetic risk for schizophrenia. This dysfunction is also reflected in high-dimensional (voxel-level) space in regions with weak functional connectivity to corresponding networks. Conclusions: Our method can effectively capture spatially dynamic networks, detect nuanced SZ effects, and reveal the intricate relationship of dynamic information to genomic data. The results also underscore the potential of dynamic spatial dependence and weak connectivity in the clinical landscape.