Classification of schizophrenia and bipolar patients using static and dynamic resting-state fMRI brain connectivity

Machine learning with multiple modalities of brain magnetic resonance imaging data to identify the presence of bipolar disorder

BACKGROUND

Bipolar disorder (BD) is a chronic psychiatric mood disorder that is solely diagnosed based on clinical symptoms. These symptoms often overlap with other psychiatric disorders. Efforts to use machine learning (ML) to create predictive models for BD based on data from brain imaging are expanding but have often been limited using only a single modality and the exclusion of the cerebellum, which may be relevant in BD.

METHODS

In this study, we sought to improve ML classification of BD by combining information from structural, functional, and diffusion-weighted imaging. Participants (108 BD I, 78 control) with BD type I and matched controls were recruited into an imaging study. This dataset was randomly divided into training and testing sets. For each of the three modalities, a separate ML model was selected, trained, and then used to generate a prediction of the class of each test subject. Majority voting was used to combine results from the three models to make a final prediction of whether a subject had BD. An independent replication sample was used to evaluate the ability of the ML classification to generalize to data collected at other sites.

RESULTS

Combining the three machine learning models through majority voting resulted in an accuracy of 89.5 % for classification of the test subjects as being in the BD or control group. Bootstrapping resulted in a 95 % confidence interval of 78.9 %-97.4 % for test accuracy. Performance was reduced when only using 2 of the 3 modalities. Analysis of feature importance revealed that the cerebellum and nodes of the emotional control network were among the most important regions for classification. The machine learning model performed at chance on the independent replication sample.

CONCLUSION

BD I could be identified with high accuracy in our relatively small sample by combining structural, functional, and diffusion-weighted imaging data within a single site but not generalize well to an independent replication sample. Future studies using harmonized imaging protocols may facilitate generalization of ML models.

Enhancing Prediction of Human Traits and Behaviors through Ensemble Learning of Traditional and Novel Resting-State fMRI Connectivity Analyses

Recent advances in cognitive neuroscience have focused on using resting-state functional connectivity (RSFC) data from fMRI scans to more accurately predict human traits and behaviors. Traditional approaches generally analyze RSFC by correlating averaged time-series data across regions of interest (ROIs) or networks, which may overlook important spatial signal patterns. To address this limitation, we introduced a novel linear regression technique that estimates RSFC by predicting spatial brain activity patterns in a target ROI from those in a seed ROI. We applied both traditional and our novel RSFC estimation methods to a large-scale dataset from the Human Connectome Project and the Brain Genomics Superstruct Project, analyzing resting-state fMRI data to predict sex, age, personality traits, and psychological task performance. To enhance prediction accuracy, we developed an ensemble learner that combines these qualitatively different methods using a weighted average approach. Our findings revealed that hierarchical clustering of RSFC patterns using our novel method displays distinct whole-brain grouping patterns compared to the traditional approach. Importantly, the ensemble model, integrating these diverse weak learners, outperformed the traditional RSFC method in predicting human traits and behaviors. Notably, the predictions from the traditional and novel methods showed relatively low similarity, indicating that our novel approach captures unique and previously undetected information about human traits and behaviors through fine-grained local spatial patterns of neural activation. These results highlight the potential of combining traditional and innovative RSFC analysis techniques to enrich our understanding of the neural basis of human traits and behaviors.

Functional dysconnectivity of visual and somatomotor networks yields a simple and robust biomarker for psychosis

People with psychosis exhibit thalamo-cortical hyperconnectivity and cortico-cortical hypoconnectivity with sensory networks, however, it remains unclear if this applies to all sensory networks, whether it arises from other illness factors, or whether such differences could form the basis of a viable biomarker. To address the foregoing, we harnessed data from the Human Connectome Early Psychosis Project and computed resting-state functional connectivity (RSFC) matrices for 54 healthy controls and 105 psychosis patients. Primary visual, secondary visual ("visual2"), auditory, and somatomotor networks were defined via a recent brain network partition. RSFC was determined for 718 regions via regularized partial correlation. Psychosis patients-both affective and non-affective-exhibited cortico-cortical hypoconnectivity and thalamo-cortical hyperconnectivity in somatomotor and visual2 networks but not in auditory or primary visual networks. When we averaged and normalized the visual2 and somatomotor network connections, and subtracted the thalamo-cortical and cortico-cortical connectivity values, a robust psychosis biomarker emerged (p = 2e-10, Hedges' g = 1.05). This "somato-visual" biomarker was present in antipsychotic-naive patients and did not depend on confounds such as psychiatric comorbidities, substance/nicotine use, stress, anxiety, or demographics. It had moderate test-retest reliability (ICC = 0.62) and could be recovered in five-minute scans. The marker could discriminate groups in leave-one-site-out cross-validation (AUC = 0.79) and improve group classification upon being added to a well-known neurocognition task. Finally, it could differentiate later-stage psychosis patients from healthy or ADHD controls in two independent data sets. These results introduce a simple and robust RSFC biomarker that can distinguish psychosis patients from controls by the early illness stages.

Dynamic functional connectivity in anorexia nervosa: alterations in states of low connectivity and state transitions

BACKGROUND

The onset of anorexia nervosa (AN) frequently occurs during adolescence and is associated with preoccupation with body weight and shape and extreme underweight. Altered resting state functional connectivity in the brain has been described in individuals with AN, but only from a static perspective. The current study investigated the temporal dynamics of functional connectivity in adolescents with AN and how it relates to clinical features.

METHOD

99 female patients acutely ill with AN and 99 pairwise age-matched female healthy control (HC) participants were included in the study. Using resting-state functional MRI data and an established sliding-window analytic approach, we identified dynamic resting-state functional connectivity states and extracted dynamic indices such as dwell time (the duration spent in a state), fraction time (the proportion of the total time occupied by a state), and number of transitions (number of switches) from one state to another, to test for group differences.

RESULTS

Individuals with AN had relatively reduced fraction time in a mildly connected state with pronounced connectivity within the default mode network (DMN) and an overall reduced number of transitions between states.

CONCLUSIONS

These findings revealed by a dynamic, but not static analytic approach might hint towards a more "rigid" connectivity, a phenomenon commonly observed in internalizing mental disorders, and in AN possibly related to a reduction in energetic costs as a result of nutritional deprivation.

Abnormal temporal variability of thalamo-cortical circuit in patients with moderate-to-severe obstructive sleep apnea

This study investigated the abnormal dynamic functional connectivity (dFC) variability of the thalamo-cortical circuit in patients with obstructive sleep apnea (OSA) and explored the relationship between these changes and the clinical characteristics of patients with OSA. A total of 91 newly diagnosed patients with moderate-to-severe OSA and 84 education-matched healthy controls (HCs) were included. All participants underwent neuropsychological testing and a functional magnetic resonance imaging scan. We explored the thalamo-cortical dFC changes by dividing the thalamus into 16 subregions and combining them using a sliding-window approach. Correlation analysis assessed the relationship between dFC variability and clinical features, and the support vector machine method was used for classification. The OSA group exhibited increased dFC variability between the thalamic subregions and extensive cortical areas, compared with the HCs group. Decreased dFC variability was observed in some frontal-occipital-temporal cortical regions. These dFC changes positively correlated with daytime sleepiness, disease severity, and cognitive scores. Altered dFC variability contributed to the discrimination between patients with OSA and HCs, with a classification accuracy of 77.8%. Our findings show thalamo-cortical overactivation and disconnection in patients with OSA, disrupting information flow within the brain networks. These results enhance understanding of the temporal variability of thalamo-cortical circuits in patients with OSA.

Functional dysconnectivity of visual and somatomotor networks yields a simple and robust biomarker for psychosis

People with psychosis exhibit thalamo-cortical hyperconnectivity and cortico-cortical hypoconnectivity with sensory networks, however, it remains unclear if this applies to all sensory networks, whether it arises from other illness factors, or whether such differences could form the basis of a viable biomarker. To address the foregoing, we harnessed data from the Human Connectome Early Psychosis Project and computed resting-state functional connectivity (RSFC) matrices for 54 healthy controls and 105 psychosis patients. Primary visual, secondary visual ("visual2"), auditory, and somatomotor networks were defined via a recent brain network partition. RSFC was determined for 718 regions via regularized partial correlation. Psychosis patients- both affective and non-affective-exhibited cortico-cortical hypoconnectivity and thalamo-cortical hyperconnectivity in somatomotor and visual2 networks but not in auditory or primary visual networks. When we averaged and normalized the visual2 and somatomotor network connections, and subtracted the thalamo-cortical and cortico-cortical connectivity values, a robust psychosis biomarker emerged (p=2e-10, Hedges' g=1.05). This "somato-visual" biomarker was present in antipsychotic-naive patients and did not depend on confounds such as psychiatric comorbidities, substance/nicotine use, stress, anxiety, or demographics. It had moderate test-retest reliability (ICC=.61) and could be recovered in five-minute scans. The marker could discriminate groups in leave-one-site-out cross-validation (AUC=.79) and improve group classification upon being added to a well-known neurocognition task. Finally, it could differentiate later-stage psychosis patients from healthy or ADHD controls in two independent data sets. These results introduce a simple and robust RSFC biomarker that can distinguish psychosis patients from controls by the early illness stages.

Classifying Unconscious, Psychedelic, and Neuropsychiatric Brain States with Functional Connectivity, Graph Theory, and Cortical Gradient Analysis

Accurate and generalizable classification of brain states is essential for understanding their neural underpinnings and improving clinical diagnostics. Traditionally, functional connectivity patterns and graph-theoretic metrics have been utilized. However, cortical gradient features, which reflect global brain organization, offer a complementary approach. We hypothesized that a machine learning model integrating these three feature sets would effectively discriminate between baseline and atypical brain states across a wide spectrum of conditions, even though the underlying neural mechanisms vary. To test this, we extracted features from brain states associated with three meta-conditions including unconsciousness (NREM2 sleep, propofol deep sedation, and propofol general anesthesia), psychedelic states induced by hallucinogens (subanesthetic ketamine, lysergic acid diethylamide, and nitrous oxide), and neuropsychiatric disorders (attention-deficit hyperactivity disorder, bipolar disorder, and schizophrenia). We used support vector machine with nested cross-validation to construct our models. The soft voting ensemble model marked the average balanced accuracy (average of specificity and sensitivity) of 79% (62-98% across all conditions), outperforming individual base models (70-76%). Notably, our models exhibited varying degrees of transferability across different datasets, with performance being dependent on the specific brain states and feature sets used. Feature importance analysis across meta-conditions suggests that the underlying neural mechanisms vary significantly, necessitating tailored approaches for accurate classification of specific brain states. This finding underscores the value of our feature-integrated ensemble models, which leverage the strengths of multiple feature types to achieve robust performance across a broader range of brain states. While our approach offers valuable insights into the neural signatures of different brain states, future work is needed to develop and validate even more generalizable models that can accurately classify brain states across a wider array of conditions.

Joint multi-site domain adaptation and multi-modality feature selection for the diagnosis of psychiatric disorders

Identifying biomarkers for computer-aided diagnosis (CAD) is crucial for early intervention of psychiatric disorders. Multi-site data have been utilized to increase the sample size and improve statistical power, while multi-modality classification offers significant advantages over traditional single-modality based approaches for diagnosing psychiatric disorders. However, inter-site heterogeneity and intra-modality heterogeneity present challenges to multi-site and multi-modality based classification. In this paper, brain functional and structural networks (BFNs/BSNs) from multiple sites were constructed to establish a joint multi-site multi-modality framework for psychiatric diagnosis. To do this we developed a hypergraph based multi-source domain adaptation (HMSDA) which allowed us to transform source domain subjects into a target domain. A local ordinal structure based multi-task feature selection (LOSMFS) approach was developed by integrating the transformed functional and structural connections (FCs/SCs). The effectiveness of our method was validated by evaluating diagnosis of both schizophrenia (SZ) and autism spectrum disorder (ASD). The proposed method obtained accuracies of 92.2 %±2.22 % and 84.8 %±2.68 % for the diagnosis of SZ and ASD, respectively. We also compared with 6 DA, 10 multi-modality feature selection, and 8 multi-site and multi-modality methods. Results showed the proposed HMSDA+LOSMFS effectively integrated multi-site and multi-modality data to enhance psychiatric diagnosis and identify disorder-specific diagnostic brain connections.

Diverging functional connectivity timescales: Capturing distinct aspects of cognitive performance in early psychosis

BACKGROUND

Psychosis spectrum disorders (PSDs) are marked by cognitive impairments, the neurobiological correlates of which remain poorly understood. Here, we investigate the entropy of time-varying functional connectivity (TVFC) patterns from resting-state functional magnetic resonance imaging (rs-fMRI) as potential biomarker for cognitive performance in PSDs. By combining our results with multimodal reference data, we hope to generate new insights into the mechanisms underlying cognitive dysfunction in PSDs. We hypothesized that low-entropy TVFC patterns (LEN) would be more behaviorally informative than high-entropy TVFC patterns (HEN), especially for tasks that require extensive integration across diverse cognitive subdomains.

METHODS

rs-fMRI and behavioral data from 97 patients in the early phases of psychosis and 53 controls were analyzed. Positron emission tomography (PET) and magnetoencephalography (MEG) data were taken from a public repository (Hansen et al., 2022). Multivariate analyses were conducted to examine relationships between TVFC patterns at multiple spatial scales and cognitive performance in patients.

RESULTS

Compared to HEN, LEN explained significantly more cognitive variance on average in PSD patients, driven by superior encoding of information on psychometrically more integrated tasks. HEN better captured information in specific subdomains of executive functioning. Nodal HEN-LEN transitions were spatially aligned with neurobiological gradients reflecting monoaminergic transporter densities and MEG beta-power. Exploratory analyses revealed a close statistical relationship between LEN and positive symptom severity in patients.

CONCLUSION

Our entropy-based analysis of TVFC patterns dissociates distinct aspects of cognition in PSDs. By linking topographies of neurotransmission and oscillatory dynamics with cognitive performance, it enhances our understanding of the mechanisms underlying cognitive deficits in PSDs.

Continuous Dictionary of Nodes Model and Bilinear-Diffusion Representation Learning for Brain Disease Analysis

Brain networks based on functional magnetic resonance imaging (fMRI) provide a crucial perspective for diagnosing brain diseases. Representation learning has recently attracted tremendous attention due to its strong representation capability, which can be naturally applied to brain disease analysis. However, traditional representation learning only considers direct and local node interactions in original brain networks, posing challenges in constructing higher-order brain networks to represent indirect and extensive node interactions. To address this problem, we propose the Continuous Dictionary of Nodes model and Bilinear-Diffusion (CDON-BD) network for brain disease analysis. The CDON model is innovatively used to learn the original brain network, with its encoder weights directly regarded as latent features. To fully integrate latent features, we further utilize Bilinear Pooling to construct higher-order brain networks. The Diffusion Module is designed to capture extensive node interactions in higher-order brain networks. Compared to state-of-the-art methods, CDON-BD demonstrates competitive classification performance on two real datasets. Moreover, the higher-order representations learned by our method reveal brain regions relevant to the diseases, contributing to a better understanding of the pathology of brain diseases.

Neural signatures of default mode network subsystems in first-episode, drug-naive patients with major depressive disorder after 6-week thought induction psychotherapy treatment

Evidence indicates that the default mode network (DMN) plays a crucial role in the neuropathology of major depressive disorder (MDD). However, the neural signatures of DMN subsystems in MDD after low resistance Thought Induction Psychotherapy (TIP) remain incompletely understood. We collected functional magnetic resonance imaging data from 20 first-episode, drug-naive MDD and 20 healthy controls (HCs). The DMN was segmented into three subsystems and seed-based functional connectivity (FC) was computed. After 6-week treatment, the significantly reduced FCs with the medial temporal lobe memory subsystem in MDD at baseline were enhanced and were comparable to that in HCs. Changed Hamilton Depression Rating Scale scores were significantly related with changed FC between the posterior cingulate cortex (PCC) and the right precuneus (PCUN). Further, changed serotonin 5-hydroxytryptamine levels were significantly correlated with changed FCs between the PCC and the left PCUN, between the posterior inferior parietal lobule and the left inferior temporal gyrus, and between the retrosplenial cortex and the right inferior frontal gyrus, opercular part. Finally, the support vector machine obtained an accuracy of 67.5% to distinguish between MDD at baseline and HCs. These findings may deepen our understanding of the neural basis of the effects of TIP on DMN subsystems in MDD.

Abnormal large-scale brain functional network dynamics in social anxiety disorder

AIMS

Although static abnormalities of functional brain networks have been observed in patients with social anxiety disorder (SAD), the brain connectome dynamics at the macroscale network level remain obscure. We therefore used a multivariate data-driven method to search for dynamic functional network connectivity (dFNC) alterations in SAD.

METHODS

We conducted spatial independent component analysis, and used a sliding-window approach with a k-means clustering algorithm, to characterize the recurring states of brain resting-state networks; then state transition metrics and FNC strength in the different states were compared between SAD patients and healthy controls (HC), and the relationship to SAD clinical characteristics was explored.

RESULTS

Four distinct recurring states were identified. Compared with HC, SAD patients demonstrated higher fractional windows and mean dwelling time in the highest-frequency State 3, representing "widely weaker" FNC, but lower in States 2 and 4, representing "locally stronger" and "widely stronger" FNC, respectively. In State 1, representing "widely moderate" FNC, SAD patients showed decreased FNC mainly between the default mode network and the attention and perceptual networks. Some aberrant dFNC signatures correlated with illness duration.

CONCLUSION

These aberrant patterns of brain functional synchronization dynamics among large-scale resting-state networks may provide new insights into the neuro-functional underpinnings of SAD.

Dynamic properties in functional connectivity changes and striatal dopamine deficiency in Parkinson's disease

Recent studies in Parkinson's disease (PD) patients reported disruptions in dynamic functional connectivity (dFC, i.e., a characterization of spontaneous fluctuations in functional connectivity over time). Here, we assessed whether the integrity of striatal dopamine terminals directly modulates dFC metrics in two separate PD cohorts, indexing dopamine-related changes in large-scale brain network dynamics and its implications in clinical features. We pooled data from two disease-control cohorts reflecting early PD. From the Parkinson's Progression Marker Initiative (PPMI) cohort, resting-state functional magnetic resonance imaging (rsfMRI) and dopamine transporter (DaT) single-photon emission computed tomography (SPECT) were available for 63 PD patients and 16 age- and sex-matched healthy controls. From the clinical research group 219 (KFO) cohort, rsfMRI imaging was available for 52 PD patients and 17 age- and sex-matched healthy controls. A subset of 41 PD patients and 13 healthy control subjects additionally underwent 18F-DOPA-positron emission tomography (PET) imaging. The striatal synthesis capacity of 18F-DOPA PET and dopamine terminal quantity of DaT SPECT images were extracted for the putamen and the caudate. After rsfMRI pre-processing, an independent component analysis was performed on both cohorts simultaneously. Based on the derived components, an individual sliding window approach (44 s window) and a subsequent k-means clustering were conducted separately for each cohort to derive dFC states (reemerging intra- and interindividual connectivity patterns). From these states, we derived temporal metrics, such as average dwell time per state, state attendance, and number of transitions and compared them between groups and cohorts. Further, we correlated these with the respective measures for local dopaminergic impairment and clinical severity. The cohorts did not differ regarding age and sex. Between cohorts, PD groups differed regarding disease duration, education, cognitive scores and L-dopa equivalent daily dose. In both cohorts, the dFC analysis resulted in three distinct states, varying in connectivity patterns and strength. In the PPMI cohort, PD patients showed a lower state attendance for the globally integrated (GI) state and a lower number of transitions than controls. Significantly, worse motor scores (Unified Parkinson's Disease Rating Scale Part III) and dopaminergic impairment in the putamen and the caudate were associated with low average dwell time in the GI state and a low total number of transitions. These results were not observed in the KFO cohort: No group differences in dFC measures or associations between dFC variables and dopamine synthesis capacity were observed. Notably, worse motor performance was associated with a low number of bidirectional transitions between the GI and the lesser connected (LC) state across the PD groups of both cohorts. Hence, in early PD, relative preservation of motor performance may be linked to a more dynamic engagement of an interconnected brain state. Specifically, those large-scale network dynamics seem to relate to striatal dopamine availability. Notably, most of these results were obtained only for one cohort, suggesting that dFC is impacted by certain cohort features like educational level, or disease severity. As we could not pinpoint these features with the data at hand, we suspect that other, in our case untracked, demographical features drive connectivity dynamics in PD. PRACTITIONER POINTS: Exploring dopamine's role in brain network dynamics in two Parkinson's disease (PD) cohorts, we unraveled PD-specific changes in dynamic functional connectivity. Results in the Parkinson's Progression Marker Initiative (PPMI) and the KFO cohort suggest motor performance may be linked to a more dynamic engagement and disengagement of an interconnected brain state. Results only in the PPMI cohort suggest striatal dopamine availability influences large-scale network dynamics that are relevant in motor control.

Applications of MRI in Schizophrenia: Current Progress in Establishing Clinical Utility

Schizophrenia is a severe mental illness that significantly impacts the lives of affected individuals and with increasing mortality rates. Early detection and intervention are crucial for improving outcomes but the lack of validated biomarkers poses great challenges in such efforts. The use of magnetic resonance imaging (MRI) in schizophrenia enables the investigation of the disorder's etiological and neuropathological substrates in vivo. After decades of research, promising findings of MRI have been shown to aid in screening high-risk individuals and predicting illness onset, and predicting symptoms and treatment outcomes of schizophrenia. The integration of machine learning and deep learning techniques makes it possible to develop intelligent diagnostic and prognostic tools with extracted or selected imaging features. In this review, we aimed to provide an overview of current progress and prospects in establishing clinical utility of MRI in schizophrenia. We first provided an overview of MRI findings of brain abnormalities that might underpin the symptoms or treatment response process in schizophrenia patients. Then, we summarized the ongoing efforts in the computer-aided utility of MRI in schizophrenia and discussed the gap between MRI research findings and real-world applications. Finally, promising pathways to promote clinical translation were provided. EVIDENCE LEVEL: 2 TECHNICAL EFFICACY: Stage 3.

The overlap across psychotic disorders: A functional network connectivity analysis

Functional network connectivity (FNC) has previously been shown to distinguish patient groups from healthy controls (HC). However, the overlap across psychiatric disorders such as schizophrenia (SZ), bipolar (BP), and schizoaffective disorder (SAD) is not evident yet. This study focuses on studying the overlap across these three psychotic disorders in both dynamic and static FNC (dFNC/sFNC). We used resting-state fMRI, demographics, and clinical information from the Bipolar-Schizophrenia Network on Intermediate Phenotypes cohort (BSNIP). The data includes three groups of patients with schizophrenia (SZ, N = 181), bipolar (BP, N = 163), and schizoaffective (SAD, N = 130) and HC (N = 238) groups. After estimating each individual's dFNC, we group them into three distinct states. We evaluated two dFNC features, including occupancy rate (OCR) and distance travelled over time. Finally, the extracted features, including both sFNC and dFNC, are tested statistically across patients and HC groups. In addition, we explored the link between the clinical scores and the extracted features. We evaluated the connectivity patterns and their overlap among SZ, BP, and SAD disorders (false discovery rate or FDR corrected p < 0.05). Results showed dFNC captured unique information about overlap across disorders where all disorder groups showed similar pattern of activity in state 2. Moreover, the results showed similar patterns between SZ and SAD in state 1 which was different than BP. Finally, the distance travelled feature of SZ (average R = 0.245, p < 0.01) and combined distance travelled from all disorders was predictive of the PANSS symptoms scores (average R = 0.147, p < 0.01).

Bright light therapy-induced improvements of mood, cognitive functions and cerebellar functional connectivity in subthreshold depression: A randomized controlled trial

Background

The efficacy of bright light therapy (BLT) in ameliorating depression has been validated. The present study is to investigate the changes of depressive symptoms, cognitive function and cerebellar functional connectivity (FC) following BLT in individuals with subthreshold depression (StD).

Method

Participants were randomly assigned to BLT group (N = 47) or placebo (N = 41) in this randomized controlled trial between March 2020 and June 2022. Depression severity and cognitive function were assessed, as well as resting-state functional MRI scan was conducted before and after 8-weeks treatment. Seed-based whole-brain static FC (sFC) and dynamic FC (dFC) analyses of the bilateral cerebellar subfields were conducted. Besides, a multivariate regression model examined whether baseline brain FC was associated with changes of depression severity and cognitive function during BLT treatment.

Results

After 8-week BLT treatment, individuals with StD showed improved depressive symptoms and attention/vigilance cognitive function. BLT also increased sFC between the right cerebellar lobule IX and left temporal pole, and decreased sFC within the cerebellum, and dFC between the right cerebellar lobule IX and left medial prefrontal cortex. Moreover, the fusion of sFC and dFC at baseline could predict the improvement of attention/vigilance in response to BLT.

Conclusions

The current study identified that BLT improved depressive symptoms and attention/vigilance, as well as changed cerebellum-DMN connectivity, especially in the cerebellar-frontotemporal and cerebellar internal FC. In addition, the fusion features of sFC and dFC at pre-treatment could serve as an imaging biomarker for the improvement of attention/vigilance cognitive function after BLT in StD.

Estimation of static and dynamic functional connectivity in resting-state fMRI using zero-frequency resonator

Resting-state functional magnetic resonance imaging (rs-fMRI) is increasingly being used to infer the functional organization of the brain. Blood oxygen level-dependent (BOLD) features related to spontaneous neuronal activity, are yet to be clearly understood. Prior studies have hypothesized that rs-fMRI is spontaneous event-related and these events convey crucial information about the neuronal activity in estimating resting state functional connectivity (FC). Attempts have been made to extract these temporal events using a predetermined threshold. However, the thresholding methods in addition to being very sensitive to noise, may consider redundant events or exclude the low-valued inflection points. Here, we extract the event-related temporal onsets from the rs-fMRI time courses using a zero-frequency resonator (ZFR). The ZFR reflects the transient behavior of the BOLD events at its output. The conditional rate (CR) of the BOLD events occurring in a time course with respect to a seed time course is used to derive static FC. The temporal activity around the estimated events called high signal-to-noise ratio (SNR) segments are also obtained in the rs-fMRI time course and are then used to compute static and dynamic FCs during rest. Coactivation pattern (CAP) is the dynamic FC obtained using the high SNR segments driven by the ZFR. The static FC demonstrates that the ZFR-based CR distinguishes the coactivation and non-coactivation scores well in the distribution. CAP analysis demonstrated the stable and longer dwell time dominant resting state functional networks with high SNR segments driven by the ZFR. Static and dynamic FC analysis underpins that the ZFR-driven temporal onsets of BOLD events derive reliable and consistent FCs in the resting brain using a subset of the time points.

A Demographic-Conditioned Variational Autoencoder for fMRI Distribution Sampling and Removal of Confounds

Objective

fMRI and derived measures such as functional connectivity (FC) have been used to predict brain age, general fluid intelligence, psychiatric disease status, and preclinical neurodegenerative disease. However, it is not always clear that all demographic confounds, such as age, sex, and race, have been removed from fMRI data. Additionally, many fMRI datasets are restricted to authorized researchers, making dissemination of these valuable data sources challenging.

Methods

We create a variational autoencoder (VAE)-based model, DemoVAE, to decorrelate fMRI features from demographics and generate high-quality synthetic fMRI data based on user-supplied demographics. We train and validate our model using two large, widely used datasets, the Philadelphia Neurodevel-opmental Cohort (PNC) and Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP).

Results

We find that DemoVAE recapitulates group differences in fMRI data while capturing the full breadth of individual variations. Significantly, we also find that most clinical and computerized battery fields that are correlated with fMRI data are not correlated with DemoVAE latents. An exception are several fields related to schizophrenia medication and symptom severity.

Conclusion

Our model generates fMRI data that captures the full distribution of FC better than traditional VAE or GAN models. We also find that most prediction using fMRI data is dependent on correlation with, and prediction of, demographics.

Significance

Our DemoVAE model allows for generation of high quality synthetic data conditioned on subject demographics as well as the removal of the confounding effects of demographics. We identify that FC-based prediction tasks are highly influenced by demographic confounds.

A Demographic-Conditioned Variational Autoencoder for fMRI Distribution Sampling and Removal of Confounds

Objective

fMRI and derived measures such as functional connectivity (FC) have been used to predict brain age, general fluid intelligence, psychiatric disease status, and preclinical neurodegenerative disease. However, it is not always clear that all demographic confounds, such as age, sex, and race, have been removed from fMRI data. Additionally, many fMRI datasets are restricted to authorized researchers, making dissemination of these valuable data sources challenging.

Methods

We create a variational autoencoder (VAE)-based model, DemoVAE, to decorrelate fMRI features from demographics and generate high-quality synthetic fMRI data based on user-supplied demographics. We train and validate our model using two large, widely used datasets, the Philadelphia Neurodevelopmental Cohort (PNC) and Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP).

Results

We find that DemoVAE recapitulates group differences in fMRI data while capturing the full breadth of individual variations. Significantly, we also find that most clinical and computerized battery fields that are correlated with fMRI data are not correlated with DemoVAE latents. An exception are several fields related to schizophrenia medication and symptom severity.

Conclusion

Our model generates fMRI data that captures the full distribution of FC better than traditional VAE or GAN models. We also find that most prediction using fMRI data is dependent on correlation with, and prediction of, demographics.

Significance

Our DemoVAE model allows for generation of high quality synthetic data conditioned on subject demographics as well as the removal of the confounding effects of demographics. We identify that FC-based prediction tasks are highly influenced by demographic confounds.

Diverging functional connectivity timescales: Capturing distinct aspects of cognitive performance in early psychosis

Background

Psychosis spectrum disorders (PSDs) are marked by cognitive impairments, the neurobiological correlates of which remain poorly understood. Here, we investigate the entropy of time-varying functional connectivity (TVFC) patterns from resting-state fMRI (rfMRI) as potential biomarker for cognitive performance in PSDs. By combining our results with multimodal reference data, we hope to generate new insights into the mechanisms underlying cognitive dysfunction in PSDs. We hypothesized that low-entropy TVFC patterns (LEN) would be more behaviorally informative than high-entropy TVFC patterns (HEN), especially for tasks that require extensive integration across diverse cognitive subdomains.

Methods

rfMRI and behavioral data from 97 patients in the early phases of psychosis and 53 controls were analyzed. Positron-Emission Tomography (PET) and magnetoencephalography (MEG) data were taken from a public repository (Hansen et al., 2022). Multivariate analyses were conducted to examine relationships between TVFC patterns at multiple spatial scales and cognitive performance in patients.

Results

Compared to HEN, LEN explained significantly more cognitive variance on average in PSD patients, driven by superior encoding of information on psychometrically more integrated tasks. HEN better captured information in specific subdomains of executive functioning. Nodal HEN-LEN transitions were spatially aligned with neurobiological gradients reflecting monoaminergic transporter densities and MEG beta power. Exploratory analyses revealed a close statistical relationship between LEN and positive PSD symptoms.

Conclusion

Our entropy-based analysis of TVFC patterns dissociates distinct aspects of cognition in PSDs. By linking topographies of neurotransmission and oscillatory dynamics with cognitive performance, it enhances our understanding of the mechanisms underlying cognitive deficits in PSDs.

CRediT Authorship Contribution Statement

Fabian Hirsch: Conceptualization, Methodology, Software, Formal analysis, Writing - Original Draft, Writing - Review & Editing, Visualization; Ângelo Bumanglag: Methodology, Software, Formal analysis, Writing - Review & Editing; Yifei Zhang: Methodology, Software, Formal analysis, Writing - Review & Editing; Afra Wohlschlaeger: Methodology, Writing - Review & Editing, Supervision, Project administration.

Resting-state EEG dynamic functional connectivity distinguishes non-psychotic major depression, psychotic major depression and schizophrenia

This study aims to identify dynamic patterns within the spatiotemporal feature space that are specific to nonpsychotic major depression (NPMD), psychotic major depression (PMD), and schizophrenia (SCZ). The study also evaluates the effectiveness of machine learning algorithms based on these network manifestations in differentiating individuals with NPMD, PMD, and SCZ. A total of 579 participants were recruited, including 152 patients with NPMD, 45 patients with PMD, 185 patients with SCZ, and 197 healthy controls (HCs). A dynamic functional connectivity (DFC) approach was employed to estimate the principal FC states within each diagnostic group. Incremental proportions of data (ranging from 10% to 100%) within each diagnostic group were used for variability testing. DFC metrics, such as proportion, mean duration, and transition number, were examined among the four diagnostic groups to identify disease-related neural activity patterns. These patterns were then used to train a two-layer classifier for the four groups (HC, NPMD, PMD, and SCZ). The four principal brain states (i.e., states 1,2,3, and 4) identified by the DFC approach were highly representative within and across diagnostic groups. Between-group comparisons revealed significant differences in network metrics of state 2 and state 3, within delta, theta, and gamma frequency bands, between healthy individuals and patients in each diagnostic group (p < 0.01, FDR corrected). Moreover, the identified key dynamic network metrics achieved an accuracy of 73.1 ± 2.8% in the four-way classification of HC, NPMD, PMD, and SCZ, outperforming the static functional connectivity (SFC) approach (p < 0.001). These findings suggest that the proposed DFC approach can identify dynamic network biomarkers at the single-subject level. These biomarkers have the potential to accurately differentiate individual subjects among various diagnostic groups of psychiatric disorders or healthy controls. This work may contribute to the development of a valuable EEG-based diagnostic tool with enhanced accuracy and assistive capabilities.

Mapping the landscape: a bibliometric analysis of resting-state fMRI research on schizophrenia over the past 25 years

Schizophrenia, a multifaceted mental disorder characterized by disturbances in thought, perception, and emotion, has been extensively investigated through resting-state fMRI, uncovering changes in spontaneous brain activity among those affected. However, a bibliometric examination regarding publication trends in resting-state fMRI studies related to schizophrenia is lacking. This study obtained relevant publications from the Web of Science Core Collection spanning the period from 1998 to 2022. Data extracted from these publications included information on countries/regions, institutions, authors, journals, and keywords. The collected data underwent analysis and visualization using VOSviewer software. The primary analyses included examination of international and institutional collaborations, authorship patterns, co-citation analyses of authors and journals, as well as exploration of keyword co-occurrence and temporal trend networks. A total of 859 publications were retrieved, indicating an overall growth trend from 1998 to 2022. China and the United States emerged as the leading contributors in both publication outputs and citations, with Central South University and the University of New Mexico being identified as the most productive institutions. Vince D. Calhoun had the highest number of publications and citation counts, while Karl J. Friston was recognized as the most influential author based on co-citations. Key journals such as Neuroimage, Schizophrenia Research, Schizophrenia Bulletin, and Biological Psychiatry played pivotal roles in advancing this field. Recent popular keywords included support vector machine, antipsychotic medication, transcranial magnetic stimulation, and related terms. This study systematically synthesizes the historical development, current status, and future trends in resting-state fMRI research in schizophrenia, offering valuable insights for future research directions.

Predicting conversational satisfaction of face-to-face conversation through interpersonal similarity in resting-state functional connectivity

When conversing with an unacquainted person, if it goes well, we can obtain much satisfaction (referred to as conversational satisfaction). Can we predict how satisfied dyads will be with face-to-face conversation? To this end, we employed interpersonal similarity in whole-brain resting-state functional connectivity (RSFC), measured using functional magnetic resonance imaging before dyadic conversation. We investigated whether conversational satisfaction could be predicted from interpersonal similarity in RSFC using multivariate pattern analysis. Consequently, prediction was successful, suggesting that interpersonal similarity in RSFC is an effective neural biomarker predicting how much face-to-face conversation goes well. Furthermore, regression coefficients from predictive models suggest that both interpersonal similarity and dissimilarity contribute to good interpersonal relationships in terms of brain activity. The present study provides the potential of an interpersonal similarity approach using RSFC for understanding the foundations of human relationships and new neuroscientific insight into whether success in human interactions is predetermined.

Somatomotor-visual resting state functional connectivity increases after 2 years in the UK Biobank longitudinal cohort

Purpose

Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, in which high connectivity among all brain regions changes to a more modular structure with maturation. We examine FC changes in older adults after 2 years of aging in the UK Biobank (UKB) longitudinal cohort.

Approach

We process fMRI connectivity data using the Power264 atlas and then test whether the average internetwork FC changes in the 2722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t -test. We also compare the ability of Power264 and UKB-provided, independent component analysis (ICA)-based FC to determine which of a longitudinal scan pair is older. Finally, we investigate cross-sectional FC changes as well as differences due to differing scanner tasks in the UKB, Philadelphia Neurodevelopmental Cohort, and Alzheimer's Disease Neuroimaging Initiative datasets.

Results

We find a 6.8% average increase in somatomotor network (SMT)-visual network (VIS) connectivity from younger to older scans (corrected p < 10 - 15 ) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all internetwork connections, the average SMT-VIS connectivity is the best predictor of relative scan age. Using the full FC and a training set of 2000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC.

Conclusions

We conclude that SMT-VIS connectivity increases with age in the UKB longitudinal cohort and that resting state FC increases with age in the UKB cross-sectional cohort.

Dynamic functional network connectivity based on spatial source phase maps of complex-valued fMRI data: Application to schizophrenia

BACKGROUND

Dynamic spatial functional network connectivity (dsFNC) has shown advantages in detecting functional alterations impacted by mental disorders using magnitude-only fMRI data. However, complete fMRI data are complex-valued with unique and useful phase information.

METHODS

We propose dsFNC of spatial source phase (SSP) maps, derived from complex-valued fMRI data (named SSP-dsFNC), to capture the dynamics elicited by the phase. We compute mutual information for connectivity quantification, employ statistical analysis and Markov chains to assess dynamics, ultimately classifying schizophrenia patients (SZs) and healthy controls (HCs) based on connectivity variance and Markov chain state transitions across windows.

RESULTS

SSP-dsFNC yielded greater dynamics and more significant HC-SZ differences, due to the use of complete brain information from complex-valued fMRI data.

COMPARISON WITH EXISTING METHODS

Compared with magnitude-dsFNC, SSP-dsFNC detected additional and meaningful connections across windows (e.g., for right frontal parietal) and achieved 14.6% higher accuracy for classifying HCs and SZs.

CONCLUSIONS

This work provides new evidence about how SSP-dsFNC could be impacted by schizophrenia, and this information could be used to identify potential imaging biomarkers for psychotic diagnosis.

Investigating brain dynamics and their association with cognitive control in opioid use disorder using naturalistic and drug cue paradigms

Objectives

Opioid use disorder (OUD) impacts millions of people worldwide. The prevalence and debilitating effects of OUD present a pressing need to understand its neural mechanisms to provide more targeted interventions. Prior studies have linked altered functioning in large-scale brain networks with clinical symptoms and outcomes in OUD. However, these investigations often do not consider how brain responses change over time. Time-varying brain network engagement can convey clinically relevant information not captured by static brain measures.

Methods

We investigated brain dynamic alterations in individuals with OUD by applying a new multivariate computational framework to movie-watching (i.e., naturalistic; N=76) and task-based (N=70) fMRI. We further probed the associations between cognitive control and brain dynamics during a separate drug cue paradigm in individuals with OUD.

Results

Compared to healthy controls (N=97), individuals with OUD showed decreased variability in the engagement of recurring brain states during movie-watching. We also found that worse cognitive control was linked to decreased variability during the rest period when no opioid-related stimuli were present.

Conclusions

These findings suggest that individuals with OUD may experience greater difficulty in effectively engaging brain networks in response to evolving internal or external demands. Such inflexibility may contribute to aberrant response inhibition and biased attention toward opioid-related stimuli, two hallmark characteristics of OUD. By incorporating temporal information, the current study introduces novel information about how brain dynamics are altered in individuals with OUD and their behavioral implications.

Altered dynamic functional and effective connectivity in drug-naive children with Tourette syndrome

Tourette syndrome (TS) is a developmental neuropsychiatric disorder characterized by repetitive, stereotyped, involuntary tics, the neurological basis of which remains unclear. Although traditional resting-state MRI (rfMRI) studies have identified abnormal static functional connectivity (FC) in patients with TS, dynamic FC (dFC) remains relatively unexplored. The rfMRI data of 54 children with TS and 46 typically developing children (TDC) were analyzed using group independent component analysis to obtain independent components (ICs), and a sliding-window approach to generate dFC matrices. All dFC matrices were clustered into two reoccurring states, the state transition metrics were obtained. We conducted Granger causality and nodal topological analyses to further investigate the brain regions that may play the most important roles in driving whole-brain switching between different states. We found that children with TS spent more time in state 2 (PFDR < 0.001), a state characterized by strong connectivity between ICs, and switched more quickly between states (PFDR = 0.025) than TDC. The default mode network (DMN) may play an important role in abnormal state transitions because the FC that changed the most between the two states was between the DMN and other networks. Additionally, the DMN had increased degree centrality, efficiency and altered causal influence on other networks. Certain alterations related to executive function (r = -0.309, P < 0.05) and tic symptom ratings (r = 0.282; 0.413, P < 0.05) may represent important aspects of the pathophysiology of TS. These findings facilitate our understanding of the neural basis for the clinical presentation of TS.

Dynamic fusion of genomics and functional network connectivity in UK biobank reveals static and time-varying SNP manifolds

Many psychiatric and neurological disorders show significant heritability, indicating strong genetic influence. In parallel, dynamic functional network connectivity (dFNC) measures functional temporal coupling between brain networks in a time-varying manner and has proven to identify disease-related changes in the brain. However, it remains largely unclear how genetic risk contributes to brain dysconnectivity that further manifests into clinical symptoms. The current work aimed to address this gap by proposing a novel joint ICA (jICA)-based "dynamic fusion" framework to identify dynamically tuned SNP manifolds by linking static SNPs to dynamic functional information of the brain. The sliding window approach was utilized to estimate four dFNC states and compute subject-level state-specific dFNC features. Each state of dFNC features were then combined with 12946 SZ risk SNPs for jICA decomposition, resulting in four parallel fusions in 32861 European ancestry individuals within the UK Biobank cohort. The identified joint SNP-dFNC components were further validated for SZ relevance in an aggregated SZ cohort, and compared for across-state similarity to indicate level of dynamism. The results supported that dynamic fusion yielded "static" and "dynamic" components (i.e., high and low across-state similarity, respectively) for SNP and dFNC modalities. As expected, the SNP components presented a mixture of static and dynamic manifolds, with the latter largely driven by fusion with dFNC. We also showed that some of the dynamic SNP manifolds uniquely elicited by fusion with state-specific dFNC features complemented each other in terms of biological interpretation. This dynamic fusion framework thus allows expanding the SNP modality to manifolds in the time dimension, which provides a unique lens to elicit unique SNP correlates of dFNC otherwise unseen, promising additional insights on how genetic risk links to disease-related dysconnectivity.

On the variability of dynamic functional connectivity assessment methods

BACKGROUND

Dynamic functional connectivity (dFC) has become an important measure for understanding brain function and as a potential biomarker. However, various methodologies have been developed for assessing dFC, and it is unclear how the choice of method affects the results. In this work, we aimed to study the results variability of commonly used dFC methods.

METHODS

We implemented 7 dFC assessment methods in Python and used them to analyze the functional magnetic resonance imaging data of 395 subjects from the Human Connectome Project. We measured the similarity of dFC results yielded by different methods using several metrics to quantify overall, temporal, spatial, and intersubject similarity.

RESULTS

Our results showed a range of weak to strong similarity between the results of different methods, indicating considerable overall variability. Somewhat surprisingly, the observed variability in dFC estimates was found to be comparable to the expected functional connectivity variation over time, emphasizing the impact of methodological choices on the final results. Our findings revealed 3 distinct groups of methods with significant intergroup variability, each exhibiting distinct assumptions and advantages.

CONCLUSIONS

Overall, our findings shed light on the impact of dFC assessment analytical flexibility and highlight the need for multianalysis approaches and careful method selection to capture the full range of dFC variation. They also emphasize the importance of distinguishing neural-driven dFC variations from physiological confounds and developing validation frameworks under a known ground truth. To facilitate such investigations, we provide an open-source Python toolbox, PydFC, which facilitates multianalysis dFC assessment, with the goal of enhancing the reliability and interpretability of dFC studies.

Linking Polygenic Risk of Schizophrenia to Variation in Magnetic Resonance Imaging Brain Measures: A Comprehensive Systematic Review

BACKGROUND AND HYPOTHESIS

Schizophrenia is highly heritable, with a polygenic effect of many genes conferring risk. Evidence on whether cumulative risk also predicts alterations in brain morphology and function is inconsistent. This systematic review examined evidence for schizophrenia polygenic risk score (sczPRS) associations with commonly used magnetic resonance imaging (MRI) measures. We expected consistent evidence to emerge for significant sczPRS associations with variation in structure and function, specifically in frontal, temporal, and insula cortices that are commonly implicated in schizophrenia pathophysiology.

STUDY DESIGN

In accordance with Preferred Reporting Items for Systematic Reviews and Meta-Analyses (PRISMA) guidelines, we searched MEDLINE, Embase, and PsycINFO for peer-reviewed studies published between January 2013 and March 2022. Studies were screened against predetermined criteria and National Institutes of Health (NIH) quality assessment tools.

STUDY RESULTS

In total, 57 studies of T1-weighted structural, diffusion, and functional MRI were included (age range = 9-80 years, Nrange = 64-76 644). We observed moderate, albeit preliminary, evidence for higher sczPRS predicting global reductions in cortical thickness and widespread variation in functional connectivity, and to a lesser extent, region-specific reductions in frontal and temporal volume and thickness. Conversely, sczPRS does not predict whole-brain surface area or gray/white matter volume. Limited evidence emerged for sczPRS associations with diffusion tensor measures of white matter microstructure in a large community sample and smaller cohorts of children and young adults. These findings were broadly consistent across community and clinical populations.

CONCLUSIONS

Our review supports the hypothesis that schizophrenia is a disorder of disrupted within and between-region brain connectivity, and points to specific whole-brain and regional MRI metrics that may provide useful intermediate phenotypes.

Statistical Brain Network Analysis

The recent fusion of network science and neuroscience has catalyzed a paradigm shift in how we study the brain and led to the field of brain network analysis. Brain network analyses hold great potential in helping us understand normal and abnormal brain function by providing profound clinical insight into links between system-level properties and health and behavioral outcomes. Nonetheless, methods for statistically analyzing networks at the group and individual levels have lagged behind. We have attempted to address this need by developing three complementary statistical frameworks-a mixed modeling framework, a distance regression framework, and a hidden semi-Markov modeling framework. These tools serve as synergistic fusions of statistical approaches with network science methods, providing needed analytic foundations for whole-brain network data. Here we delineate these approaches, briefly survey related tools, and discuss potential future avenues of research. We hope this review catalyzes further statistical interest and methodological development in the field.

Dynamic functional connectivity changes associated with psychiatric traits and cognitive deficits in Cushing's disease

Cushing's disease is a rare neuroendocrine disorder with excessive endogenous cortisol, impaired cognition, and psychiatric symptoms. Evidence from resting-state fMRI revealed the abnormalities of static brain connectivity in patients with Cushing's disease (CD patients). However, it is unknown whether the CD patients' dynamic functional connectivity would be abnormal and whether the dynamic features are associated with deficits in cognition and psychopathological symptoms. Here, we evaluated 50 patients with Cushing's disease and 57 healthy participants by using resting-state fMRI and dynamic functional connectivity (dFNC) approach. We focused on the dynamic features of default mode network (DMN), salience network (SN), and central executive network (CEN) because these are binding sites for the cognitive-affective process, as well as vital in understanding the pathophysiology of psychiatric disorders. The dFNC was further clustered into four states by k-mean clustering. CD patients showed more dwell time in State 1 but less time in State 4. Intriguingly, group differences in dwell time in these two states can explain the cognitive deficits of CD patients. Moreover, the inter-network connections between DMN and SN and the engagement time in State 4 negatively correlated with anxiety and depression but positively correlated with cognitive performance. Finally, the classifier trained by the dynamic features of these networks successfully classified CD patients from healthy participants. Together, our study revealed the dynamic features of CD patients' brains and found their associations with impaired cognition and emotional symptoms, which may open new avenues for understanding the cognitive and affective deficits induced by Cushing's disease.

Individual differences in time-varying and stationary brain connectivity during movie watching from childhood to early adulthood: Age, sex, and behavioral associations

Spatially remote brain regions exhibit dynamic functional interactions across various task conditions. While time-varying functional connectivity during movie watching shows sensitivity to movie content, stationary functional connectivity remains relatively stable across videos. These findings suggest that dynamic and stationary functional interactions may represent different aspects of brain function. However, the relationship between individual differences in time-varying and stationary connectivity and behavioral phenotypes remains elusive. To address this gap, we analyzed an open-access functional MRI dataset comprising participants aged 5-22 years, who watched two cartoon movie clips. We calculated regional brain activity, time-varying connectivity, and stationary connectivity, examining associations with age, sex, and behavioral assessments. Model comparison revealed that time-varying connectivity was more sensitive to age and sex effects compared with stationary connectivity. The preferred age models exhibited quadratic log age or quadratic age effects, indicative of inverted-U shaped developmental patterns. In addition, females showed higher consistency in regional brain activity and time-varying connectivity than males. However, in terms of behavioral predictions, only stationary connectivity demonstrated the ability to predict full-scale intelligence quotient. These findings suggest that individual differences in time-varying and stationary connectivity may capture distinct aspects of behavioral phenotypes.

Altered Brain Dynamics Across Bipolar Disorder and Schizophrenia During Rest and Task Switching Revealed by Overlapping Brain States

BACKGROUND

Individuals with bipolar disorder (BD) and schizophrenia (SCZ) show aberrant brain dynamics (i.e., altered recruitment or traversal through different brain states over time). Existing investigations of brain dynamics typically assume that one dominant brain state characterizes each time point. However, as multiple brain states likely are engaged at any given moment, this approach can obscure alterations in less prominent but critical brain states. Here, we examined brain dynamics in BD and SCZ by implementing a novel framework that simultaneously assessed the engagement of multiple brain states.

METHODS

Four recurring brain states were identified by applying nonlinear manifold learning and k-means clustering to the Human Connectome Project task-based functional magnetic resonance imaging data. We then assessed moment-to-moment state engagement in 2 independent samples of healthy control participants and patients with BD or SCZ using resting-state (N = 336) or task-based (N = 217) functional magnetic resonance imaging data. Relative state engagement and state engagement variability were extracted and compared across groups using multivariate analysis of covariance, controlling for site, medication, age, and sex.

RESULTS

Our framework identified dynamic alterations in BD and SCZ, while a state discretization approach revealed no significant group differences. Participants with BD or SCZ showed reduced state engagement variability, but not relative state engagement, across multiple brain states during resting-state and task-based functional magnetic resonance imaging. We found decreased state engagement variability in older participants and preliminary evidence suggesting an association with avolition.

CONCLUSIONS

Assessing multiple brain states simultaneously can reflect the complexity of aberrant brain dynamics in BD and SCZ, providing a more comprehensive understanding of the neural mechanisms underpinning these conditions.

Causal inference in cognitive neuroscience

Causal inference is a key step in many research endeavors in cognitive science and neuroscience, and particularly cognitive neuroscience. Statistical knowledge is sufficient for prediction and diagnosis, but causal knowledge is required for action and intervention. Most statistics courses and textbooks emphasize the difficulty of causal inference, focusing on the maxim that "correlation does not mean causation": there can be multiple causal possibilities, often many of them, consistent with given observed statistics. This paper focuses instead on the conceptual issues and assumptions that confront causal and other kinds of inference, primarily focusing on cognitive neuroscience. We connect inference methods with goals and challenges, and provide concrete guidance about how to select appropriate tools for the scientific task. This article is categorized under: Psychology > Theory and Methods Philosophy > Foundations of Cognitive Science.

Somatomotor-Visual Resting State Functional Connectivity Increases After Two Years in the UK Biobank Longitudinal Cohort

Purpose

Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, where high connectivity among all brain regions changes to a more modular structure with maturation. In this work, we examine changes in FC in older adults after two years of aging in the UK Biobank longitudinal cohort.

Approach

We process data using the Power264 atlas, then test whether FC changes in the 2,722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t-test. We also compare the ability of Power264 and UKB-provided, ICA-based FC to determine which of a longitudinal scan pair is older.

Results

We find a 6.8% average increase in SMT-VIS connectivity from younger to older scan (from ρ = 0.39 to ρ = 0.42) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all inter-network connections, this average SMT-VIS connectivity is the best predictor of relative scan age, accurately predicting which scan is older 57% of the time. Using the full FC and a training set of 2,000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC.

Conclusions

We conclude that SMT-VIS connectivity increases in the longitudinal cohort, while resting state FC increases generally with age in the cross-sectional cohort. However, we consider the possibility of a change in resting state scanner task between UKB longitudinal data acquisitions.

Somatomotor-Visual Resting State Functional Connectivity Increases After Two Years in the UK Biobank Longitudinal Cohort

Purpose

Functional magnetic resonance imaging (fMRI) and functional connectivity (FC) have been used to follow aging in both children and older adults. Robust changes have been observed in children, where high connectivity among all brain regions changes to a more modular structure with maturation. In this work, we examine changes in FC in older adults after two years of aging in the UK Biobank longitudinal cohort.

Approach

We process data using the Power264 atlas, then test whether FC changes in the 2,722-subject longitudinal cohort are statistically significant using a Bonferroni-corrected t-test. We also compare the ability of Power264 and UKB-provided, ICA-based FC to determine which of a longitudinal scan pair is older.

Results

We find a 6.8% average increase in SMT-VIS connectivity from younger to older scan (from ρ = 0.39 to ρ = 0.42 ) that occurs in male, female, older subject (> 65 years old), and younger subject (< 55 years old) groups. Among all inter-network connections, this average SMT-VIS connectivity is the best predictor of relative scan age, accurately predicting which scan is older 57% of the time. Using the full FC and a training set of 2,000 subjects, one is able to predict which scan is older 82.5% of the time using either the full Power264 FC or the UKB-provided ICA-based FC.

Conclusions

We conclude that SMT-VIS connectivity increases in the longitudinal cohort, while resting state FC increases generally with age in the cross-sectional cohort. However, we consider the possibility of a change in resting state scanner task between UKB longitudinal data acquisitions.

Identifiability in Functional Connectivity May Unintentionally Inflate Prediction Results

Functional magnetic resonance (fMRI) is an invaluable tool in studying cognitive processes in vivo. Many recent studies use functional connectivity (FC), partial correlation connectivity (PC), or fMRI-derived brain networks to predict phenotypes with results that sometimes cannot be replicated. At the same time, FC can be used to identify the same subject from different scans with great accuracy. In this paper, we show a method by which one can unknowingly inflate classification results from 61% accuracy to 86% accuracy by treating longitudinal or contemporaneous scans of the same subject as independent data points. Using the UK Biobank dataset, we find one can achieve the same level of variance explained with 50 training subjects by exploiting identifiability as with 10,000 training subjects without double-dipping. We replicate this effect in four different datasets: the UK Biobank (UKB), the Philadelphia Neurodevelopmental Cohort (PNC), the Bipolar and Schizophrenia Network for Intermediate Phenotypes (BSNIP), and an OpenNeuro Fibromyalgia dataset (Fibro). The unintentional improvement ranges between 7% and 25% in the four datasets. Additionally, we find that by using dynamic functional connectivity (dFC), one can apply this method even when one is limited to a single scan per subject. One major problem is that features such as ROIs or connectivities that are reported alongside inflated results may confuse future work. This article hopes to shed light on how even minor pipeline anomalies may lead to unexpectedly superb results.

A meta-analysis and systematic review of single vs. multimodal neuroimaging techniques in the classification of psychosis

BACKGROUND

Psychotic disorders are characterized by structural and functional abnormalities in brain networks. Neuroimaging techniques map and characterize such abnormalities using unique features (e.g., structural integrity, coactivation). However, it is unclear if a specific method, or a combination of modalities, is particularly effective in identifying differences in brain networks of someone with a psychotic disorder.

METHODS

A systematic meta-analysis evaluated machine learning classification of schizophrenia spectrum disorders in comparison to healthy control participants using various neuroimaging modalities (i.e., T1-weighted imaging (T1), diffusion tensor imaging (DTI), resting state functional connectivity (rs-FC), or some combination (multimodal)). Criteria for manuscript inclusion included whole-brain analyses and cross-validation to provide a complete picture regarding the predictive ability of large-scale brain systems in psychosis. For this meta-analysis, we searched Ovid MEDLINE, PubMed, PsychInfo, Google Scholar, and Web of Science published between inception and March 13th 2023. Prediction results were averaged for studies using the same dataset, but parallel analyses were run that included studies with pooled sample across many datasets. We assessed bias through funnel plot asymmetry. A bivariate regression model determined whether differences in imaging modality, demographics, and preprocessing methods moderated classification. Separate models were run for studies with internal prediction (via cross-validation) and external prediction.

RESULTS

93 studies were identified for quantitative review (30 T1, 9 DTI, 40 rs-FC, and 14 multimodal). As a whole, all modalities reliably differentiated those with schizophrenia spectrum disorders from controls (OR = 2.64 (95%CI = 2.33 to 2.95)). However, classification was relatively similar across modalities: no differences were seen across modalities in the classification of independent internal data, and a small advantage was seen for rs-FC studies relative to T1 studies in classification in external datasets. We found large amounts of heterogeneity across results resulting in significant signs of bias in funnel plots and Egger's tests. Results remained similar, however, when studies were restricted to those with less heterogeneity, with continued small advantages for rs-FC relative to structural measures. Notably, in all cases, no significant differences were seen between multimodal and unimodal approaches, with rs-FC and unimodal studies reporting largely overlapping classification performance. Differences in demographics and analysis or denoising were not associated with changes in classification scores.

CONCLUSIONS

The results of this study suggest that neuroimaging approaches have promise in the classification of psychosis. Interestingly, at present most modalities perform similarly in the classification of psychosis, with slight advantages for rs-FC relative to structural modalities in some specific cases. Notably, results differed substantially across studies, with suggestions of biased effect sizes, particularly highlighting the need for more studies using external prediction and large sample sizes. Adopting more rigorous and systematized standards will add significant value toward understanding and treating this critical population.

Disrupted Dynamic Functional Network Connectivity Among Cognitive Control Networks in the Progression of Alzheimer's Disease

Background: Alzheimer's disease (AD) is the most common age-related dementia that promotes a decline in memory, thinking, and social skills. The initial stages of dementia can be associated with mild symptoms, and symptom progression to a more severe state is heterogeneous across patients. Recent work has demonstrated the potential for functional network mapping to assist in the prediction of symptomatic progression. However, this work has primarily used static functional connectivity (sFC) from resting-state functional magnetic resonance imaging. Recently, dynamic functional connectivity (dFC) has been recognized as a powerful advance in functional connectivity methodology to differentiate brain network dynamics between healthy and diseased populations. Methods: Group independent component analysis was applied to extract 17 components within the cognitive control network (CCN) from 1385 individuals across varying stages of AD symptomology. We estimated dFC among 17 components within the CCN, followed by clustering the dFCs into 3 recurring brain states, and then estimated a hidden Markov model and the occupancy rate for each subject. Then, we investigated the link between CCN dFC features and AD progression. Also, we investigated the link between sFC and AD progression and compared its results with dFC results. Results: Progression of AD symptoms was associated with increases in connectivity within the middle frontal gyrus. Also, the very mild AD (vmAD) showed less connectivity within the inferior parietal lobule (in both sFC and dFC) and between this region and the rest of CCN (in dFC analysis). Also, we found that within-middle frontal gyrus connectivity increases with AD progression in both sFC and dFC results. Finally, comparing with vmAD, we found that the normal brain spends significantly more time in a state with lower within-middle frontal gyrus connectivity and higher connectivity between the hippocampus and the rest of CCN, highlighting the importance of assessing the dynamics of brain connectivity in this disease. Conclusion: Our results suggest that AD progress not only alters the CCN connectivity strength but also changes the temporal properties in this brain network. This suggests the temporal and spatial pattern of CCN as a biomarker that differentiates different stages of AD.

Manifold learning for fMRI time-varying functional connectivity

Whole-brain functional connectivity (FC) measured with functional MRI (fMRI) evolves over time in meaningful ways at temporal scales going from years (e.g., development) to seconds [e.g., within-scan time-varying FC (tvFC)]. Yet, our ability to explore tvFC is severely constrained by its large dimensionality (several thousands). To overcome this difficulty, researchers often seek to generate low dimensional representations (e.g., 2D and 3D scatter plots) hoping those will retain important aspects of the data (e.g., relationships to behavior and disease progression). Limited prior empirical work suggests that manifold learning techniques (MLTs)-namely those seeking to infer a low dimensional non-linear surface (i.e., the manifold) where most of the data lies-are good candidates for accomplishing this task. Here we explore this possibility in detail. First, we discuss why one should expect tvFC data to lie on a low dimensional manifold. Second, we estimate what is the intrinsic dimension (ID; i.e., minimum number of latent dimensions) of tvFC data manifolds. Third, we describe the inner workings of three state-of-the-art MLTs: Laplacian Eigenmaps (LEs), T-distributed Stochastic Neighbor Embedding (T-SNE), and Uniform Manifold Approximation and Projection (UMAP). For each method, we empirically evaluate its ability to generate neuro-biologically meaningful representations of tvFC data, as well as their robustness against hyper-parameter selection. Our results show that tvFC data has an ID that ranges between 4 and 26, and that ID varies significantly between rest and task states. We also show how all three methods can effectively capture subject identity and task being performed: UMAP and T-SNE can capture these two levels of detail concurrently, but LE could only capture one at a time. We observed substantial variability in embedding quality across MLTs, and within-MLT as a function of hyper-parameter selection. To help alleviate this issue, we provide heuristics that can inform future studies. Finally, we also demonstrate the importance of feature normalization when combining data across subjects and the role that temporal autocorrelation plays in the application of MLTs to tvFC data. Overall, we conclude that while MLTs can be useful to generate summary views of labeled tvFC data, their application to unlabeled data such as resting-state remains challenging.

Sampling inequalities affect generalization of neuroimaging-based diagnostic classifiers in psychiatry

BACKGROUND

The development of machine learning models for aiding in the diagnosis of mental disorder is recognized as a significant breakthrough in the field of psychiatry. However, clinical practice of such models remains a challenge, with poor generalizability being a major limitation.

METHODS

Here, we conducted a pre-registered meta-research assessment on neuroimaging-based models in the psychiatric literature, quantitatively examining global and regional sampling issues over recent decades, from a view that has been relatively underexplored. A total of 476 studies (n = 118,137) were included in the current assessment. Based on these findings, we built a comprehensive 5-star rating system to quantitatively evaluate the quality of existing machine learning models for psychiatric diagnoses.

RESULTS

A global sampling inequality in these models was revealed quantitatively (sampling Gini coefficient (G) = 0.81, p < .01), varying across different countries (regions) (e.g., China, G = 0.47; the USA, G = 0.58; Germany, G = 0.78; the UK, G = 0.87). Furthermore, the severity of this sampling inequality was significantly predicted by national economic levels (β =  - 2.75, p < .001, R²_adj = 0.40; r =  - .84, 95% CI: - .41 to - .97), and was plausibly predictable for model performance, with higher sampling inequality for reporting higher classification accuracy. Further analyses showed that lack of independent testing (84.24% of models, 95% CI: 81.0-87.5%), improper cross-validation (51.68% of models, 95% CI: 47.2-56.2%), and poor technical transparency (87.8% of models, 95% CI: 84.9-90.8%)/availability (80.88% of models, 95% CI: 77.3-84.4%) are prevailing in current diagnostic classifiers despite improvements over time. Relating to these observations, model performances were found decreased in studies with independent cross-country sampling validations (all p < .001, BF₁₀ > 15). In light of this, we proposed a purpose-built quantitative assessment checklist, which demonstrated that the overall ratings of these models increased by publication year but were negatively associated with model performance.

CONCLUSIONS

Together, improving sampling economic equality and hence the quality of machine learning models may be a crucial facet to plausibly translating neuroimaging-based diagnostic classifiers into clinical practice.

Interpretable LSTM model reveals transiently-realized patterns of dynamic brain connectivity that predict patient deterioration or recovery from very mild cognitive impairment

Alzheimer's Disease (AZD) is a neurodegenerative disease for which there is now no known effective treatment. Mild cognitive impairment (MCI) is considered a precursor to AZD and affects cognitive abilities. Patients with MCI have the potential to recover cognitive health, can remain mildly cognitively impaired indefinitely or eventually progress to AZD. Identifying imaging-based predictive biomarkers for disease progression in patients presenting with evidence of very mild/questionable MCI (qMCI) can play an important role in triggering early dementia intervention. Dynamic functional network connectivity (dFNC) estimated from resting-state functional magnetic resonance imaging (rs-fMRI) has been increasingly studied in brain disorder diseases. In this work, employing a recent developed a time-attention long short-term memory (TA-LSTM) network to classify multivariate time series data. A gradient-based interpretation framework, transiently-realized event classifier activation map (TEAM) is introduced to localize the group-defining "activated" time intervals over the full time series and generate the class difference map. To test the trustworthiness of TEAM, we did a simulation study to validate the model interpretative power of TEAM. We then applied this simulation-validated framework to a well-trained TA-LSTM model which predicts the progression or recovery from questionable/mild cognitive impairment (qMCI) subjects after three years from windowless wavelet-based dFNC (WWdFNC). The FNC class difference map points to potentially important predictive dynamic biomarkers. Moreover, the more highly time-solved dFNC (WWdFNC) achieves better performance in both TA-LSTM and a multivariate CNN model than dFNC based on windowed correlations between timeseries, suggesting that better temporally resolved measures can enhance the model's performance.

Larger lesion volume in people with multiple sclerosis is associated with increased transition energies between brain states and decreased entropy of brain activity

Quantifying the relationship between the brain's functional activity patterns and its structural backbone is crucial when relating the severity of brain pathology to disability in multiple sclerosis (MS). Network control theory (NCT) characterizes the brain's energetic landscape using the structural connectome and patterns of brain activity over time. We applied NCT to investigate brain-state dynamics and energy landscapes in controls and people with MS (pwMS). We also computed entropy of brain activity and investigated its association with the dynamic landscape's transition energy and lesion volume. Brain states were identified by clustering regional brain activity vectors, and NCT was applied to compute the energy required to transition between these brain states. We found that entropy was negatively correlated with lesion volume and transition energy, and that larger transition energies were associated with pwMS with disability. This work supports the notion that shifts in the pattern of brain activity in pwMS without disability results in decreased transition energies compared to controls, but, as this shift evolves over the disease, transition energies increase beyond controls and disability occurs. Our results provide the first evidence in pwMS that larger lesion volumes result in greater transition energy between brain states and decreased entropy of brain activity.

Individual differences in time-varying and stationary brain connectivity during movie watching from childhood to early adulthood: age, sex, and behavioral associations

Spatially remote brain regions exhibit dynamic functional interactions across various task conditions. While time-varying functional connectivity during movie watching shows sensitivity to movie content, stationary functional connectivity remains relatively stable across videos. These findings suggest that dynamic and stationary functional interactions may represent different aspects of brain function. However, the relationship between individual differences in time-varying and stationary connectivity and behavioral phenotypes remains elusive. To address this gap, we analyzed an open-access functional MRI dataset comprising participants aged 5 to 22 years, who watched two cartoon movie clips. We calculated regional brain activity, time-varying connectivity, and stationary connectivity, examining associations with age, sex, and behavioral assessments. Model comparison revealed that time-varying connectivity was more sensitive to age and sex effects compared with stationary connectivity. The preferred age models exhibited quadratic log age or quadratic age effects, indicative of inverted-U shaped developmental patterns. In addition, females showed higher consistency in regional brain activity and time-varying connectivity than males. However, in terms of behavioral predictions, only stationary connectivity demonstrated the ability to predict full-scale intelligence quotient. These findings suggest that individual differences in time-varying and stationary connectivity may capture distinct aspects of behavioral phenotypes.

Evaluation of boundaries between mood and psychosis disorder using dynamic functional network connectivity (dFNC) via deep learning classification

The validity and reliability of diagnoses in psychiatry is a challenging topic in mental health. The current mental health categorization is based primarily on symptoms and clinical course and is not biologically validated. Among multiple ongoing efforts, neurological observations alongside clinical evaluations are considered to be potential solutions to address diagnostic problems. The Bipolar-Schizophrenia Network on Intermediate Phenotypes (B-SNIP) has published multiple papers attempting to reclassify psychotic illnesses based on biological rather than symptomatic measures. However, the effort to investigate the relationship between this new categorization approach and other neuroimaging techniques, including resting-state fMRI data, is still limited. This study focused on investigating the relationship between different psychotic disorders categorization methods and resting-state fMRI-based measures called dynamic functional network connectivity (dFNC) using state-of-the-art artificial intelligence (AI) approaches. We applied our method to 613 subjects, including individuals with psychosis and healthy controls, which were classified using both the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) and the B-SNIP biomarker-based (Biotype) approach. Statistical group differences and cross-validated classifiers were performed within each framework to assess how different categories. Results highlight interesting differences in occupancy in both DSM-IV and Biotype categorizations compared to healthy individuals, which are distributed across specific transient connectivity states. Biotypes tended to show less distinctiveness in occupancy level and included fewer cellwise differences. Classification accuracy obtained by DSM-IV and Biotype categories were both well above chance. Results provided new insights and highlighted the benefits of both DSM-IV and biology-based categories while also emphasizing the importance of future work in this direction, including employing further data types.

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.

Dysconnectivity in Schizophrenia Revisited: Abnormal Temporal Organization of Dynamic Functional Connectivity in Patients With a First Episode of Psychosis

BACKGROUND AND HYPOTHESIS

Abnormal functional connectivity between brain regions is a consistent finding in schizophrenia, including functional magnetic resonance imaging (fMRI) studies. Recent studies have highlighted that connectivity changes in time in healthy subjects. We here examined the temporal changes in functional connectivity in patients with a first episode of psychosis (FEP). Specifically, we analyzed the temporal order in which whole-brain organization states were visited.

STUDY DESIGN

Two case-control studies, including in each sample a subgroup scanned a second time after treatment. Chilean sample included 79 patients with a FEP and 83 healthy controls. Mexican sample included 21 antipsychotic-naïve FEP patients and 15 healthy controls. Characteristics of the temporal trajectories between whole-brain functional connectivity meta-states were examined via resting-state functional MRI using elements of network science. We compared the cohorts of cases and controls and explored their differences as well as potential associations with symptoms, cognition, and antipsychotic medication doses.

STUDY RESULTS

We found that the temporal sequence in which patients' brain dynamics visited the different states was more redundant and segregated. Patients were less flexible than controls in changing their network in time from different configurations, and explored the whole landscape of possible states in a less efficient way. These changes were related to the dose of antipsychotics the patients were receiving. We replicated the relationship with antipsychotic medication in the antipsychotic-naïve FEP sample scanned before and after treatment.

CONCLUSIONS

We conclude that psychosis is related to a temporal disorganization of the brain's dynamic functional connectivity, and this is associated with antipsychotic medication use.

Detecting individuals with severe mental illness using artificial intelligence applied to magnetic resonance imaging

BACKGROUND

Identifying individuals at risk for severe mental illness (SMI) is crucial for prevention and early intervention strategies. While MRI shows potential for case identification even before illness onset, no practical model for mental health risk monitoring has been developed. This study aims to develop an initial version of an efficient and practical model for mental health screening among at-risk populations.

METHODS

A deep learning model known as Multiple Instance Learning (MIL) was adopted to train and test a SMI detection model with clinical MRI scans of 14,915 patients with SMI (age 32.98 ± 12.01 years, 9102 women) and 4538 healthy controls (age 40.60 ± 10.95 years, 2424 women) in the primary dataset. Validation analysis was conducted in an independent dataset with 290 patients (age 28.08 ± 10.95 years, 169 women) and 310 healthy participants (age 33.55 ± 11.09 years, 165 women). Another three machine learning models of ResNet, DenseNet and EfficientNet were used for comparison. We also recruited 148 individuals receiving high-stress medical school education to characterize the potential real-world utility of the MIL model in detecting risk of mental illness.

FINDINGS

Similar performance of successful differentiation of individuals with SMI and healthy controls was observed for the MIL model (AUC: 0.82) and other models (ResNet, DenseNet, EfficientNet, 0.83, 0.81, and 0.80 respectively). MIL had better generalization in the validation test than other models (AUC: 0.82 vs 0.59, 0.66 and 0.59), and less drop-off in performance from 3.0T to 1.5T scanners. The MIL model did better in predicting clinician ratings of distress than self-ratings with questionnaires (84% vs 22%) in the medical student sample. Brain regions that contributed to SMI identification were mainly neocortical, including right precuneus, bilateral temporal regions, left precentral/postcentral gyrus, bilateral medial prefrontal cortex and right cerebellum.

INTERPRETATION

Our digital model based on brief clinical MRI protocols identified individual SMI patients with good accuracy and high sensitivity, suggesting that with incremental improvements the approach may offer potentially useful aid for early identification and intervention to prevent illness onset in vulnerable at-risk populations.

FUNDING

This study was supported by the National Natural Science Foundation of China, National Key Technologies R&D Program of China, and Sichuan Science and Technology Program.