SAN: Mitigating spatial covariance heterogeneity in cortical thickness data collected from multiple scanners or sites.

In neuroimaging studies, combining data collected from multiple study sites or scanners is becoming common to increase the reproducibility of scientific discoveries. At the same time, unwanted variations arise by using different scanners (inter-scanner biases), which need to be corrected before downstream analyses to facilitate replicable research and prevent spurious findings. While statistical harmonization methods such as ComBat have become popular in mitigating inter-scanner biases in neuroimaging, recent methodological advances have shown that harmonizing heterogeneous covariances results in higher data quality. In vertex-level cortical thickness data, heterogeneity in spatial autocorrelation is a critical factor that affects covariance heterogeneity. Our work proposes a new statistical harmonization method called spatial autocorrelation normalization (SAN) that preserves homogeneous covariance vertex-level cortical thickness data across different scanners. We use an explicit Gaussian process to characterize scanner-invariant and scanner-specific variations to reconstruct spatially homogeneous data across scanners. SAN is computationally feasible, and it easily allows the integration of existing harmonization methods. We demonstrate the utility of the proposed method using cortical thickness data from the Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) study. SAN is publicly available as an R package.

Trajectories of remitted psychotic depression: identification of predictors of worsening by machine learning.

BACKGROUND: Remitted psychotic depression (MDDPsy) has heterogeneity of outcome. The study's aims were to identify subgroups of persons with remitted MDDPsy with distinct trajectories of depression severity during continuation treatment and to detect predictors of membership to the worsening trajectory. METHOD: One hundred and twenty-six persons aged 18-85 years participated in a 36-week randomized placebo-controlled trial (RCT) that examined the clinical effects of continuing olanzapine once an episode of MDDPsy had remitted with sertraline plus olanzapine. Latent class mixed modeling was used to identify subgroups of participants with distinct trajectories of depression severity during the RCT. Machine learning was used to predict membership to the trajectories based on participant pre-trajectory characteristics. RESULTS: Seventy-one (56.3%) participants belonged to a subgroup with a stable trajectory of depression scores and 55 (43.7%) belonged to a subgroup with a worsening trajectory. A random forest model with high prediction accuracy (AUC of 0.812) found that the strongest predictors of membership to the worsening subgroup were residual depression symptoms at onset of remission, followed by anxiety score at RCT baseline and age of onset of the first lifetime depressive episode. In a logistic regression model that examined depression score at onset of remission as the only predictor variable, the AUC (0.778) was close to that of the machine learning model. CONCLUSIONS: Residual depression at onset of remission has high accuracy in predicting membership to worsening outcome of remitted MDDPsy. Research is needed to determine how best to optimize the outcome of psychotic MDDPsy with residual symptoms.

Effects of antipsychotic medication on functional connectivity in major depressive disorder with psychotic features.

The effect of antipsychotic medication on resting state functional connectivity in major depressive disorder (MDD) is currently unknown. To address this gap, we examined patients with MDD with psychotic features (MDDPsy) participating in the Study of the Pharmacotherapy of Psychotic Depression II. All participants were treated with sertraline plus olanzapine and were subsequently randomized to continue sertraline plus olanzapine or be switched to sertraline plus placebo. Participants completed an MRI at randomization and at study endpoint (study completion at Week 36, relapse, or early termination). The primary outcome was change in functional connectivity measured within and between specified networks and the rest of the brain. The secondary outcome was change in network topology measured by graph metrics. Eighty-eight participants completed a baseline scan; 73 completed a follow-up scan, of which 58 were usable for analyses. There was a significant treatment X time interaction for functional connectivity between the secondary visual network and rest of the brain (t = -3.684; p = 0.0004; pFDR = 0.0111). There was no significant treatment X time interaction for graph metrics. Overall, functional connectivity between the secondary visual network and the rest of the brain did not change in participants who stayed on olanzapine but decreased in those switched to placebo. There were no differences in changes in network topology measures when patients stayed on olanzapine or switched to placebo. This suggests that olanzapine may stabilize functional connectivity, particularly between the secondary visual network and the rest of the brain.

Enhanced conscious processing and blindsight-like detection of fear-conditioned stimuli under continuous flash suppression.

Emotional information, and specifically fear-related stimuli, have been shown to be preferentially processed at a nonconscious level and gain privileged access to awareness. However, recent evidence has emerged suggesting these findings are explained by low-level visual features rather than emotional salience. In this study, we tested the hypothesis that emotional salience increases both blindsight (i.e., detection with reduced awareness) and awareness of visually suppressed stimuli. We used fear conditioning to manipulate the emotional significance of neutral expressions presented under Continuous Flash Suppression. Fifty-two healthy participants were presented with perithreshold neutral faces, previously paired with an electric shock (CS+) or not (CS-), and asked to localise the quadrant wherein faces were presented and rate their level of confidence in the response. Results showed fear conditioning strength (indexed by skin conductance response to CS+ versus CS-) was positively associated with both increased "blindsight" and awareness of conditioned stimuli. These findings suggest emotional significance alone, and not merely low-level visual differences, can enhance pre-conscious and conscious processing of visual stimuli.

"Blindsight" and subjective awareness of fearful faces: Inversion reverses the deficits in fear perception associated with core psychopathic traits.

Though emotional faces preferentially reach awareness, the present study utilised both objective and subjective indices of awareness to determine whether they enhance subjective awareness and "blindsight". Under continuous flash suppression, participants localised a disgusted, fearful or neutral face (objective index), and rated their confidence (subjective index). Psychopathic traits were also measured to investigate their influence on emotion perception. As predicted, fear increased localisation accuracy, subjective awareness and "blindsight" of upright faces. Coldhearted traits were inversely related to subjective awareness, but not "blindsight", of upright fearful faces. In a follow-up experiment using inverted faces, increased localisation accuracy and awareness, but not "blindsight", were observed for fear. Surprisingly, awareness of inverted fearful faces was positively correlated with coldheartedness. These results suggest that emotion enhances both pre-conscious processing and the qualitative experience of awareness, but that pre-conscious and conscious processing of emotional faces rely on at least partially dissociable cognitive mechanisms.

Toxicology-based cancer causation analysis of CoCr-containing hip implants: a quantitative assessment of genotoxicity and tumorigenicity studies.

In this paper, quantitative methods were used to evaluate the weight of evidence regarding a causative relationship between cobalt-chromium (CoCr)-containing hip implants and increased cancer risk. We reviewed approximately 80 published papers and identified no-observed-adverse-effect level (NOAEL) and/or lowest-observed-adverse-effect level (LOAEL) values for specific endpoints of interest: genotoxic effects from in vitro studies with human cell lines as well as genotoxicity and tumor formation in animal bioassays. Test articles included Co particles and ions, Cr particles and ions, and CoCr alloy particles as well as CoCr alloy implants. The NOAEL/LOAEL values were compared with body burdens of Co/Cr particles and ions we calculated to exist in systemic tissues of hip implant patients under normal and excessive wear conditions. We found that approximately 40 tumor bioassays have been conducted with CoCr alloy implants or Co/Cr particles and ions at levels hundreds to thousands of times higher than those present in hip implant patients, and none reported a statistically significant increased incidence of systemic tumors. Results from in vitro and in vivo genotoxicity assays, which are relatively less informative owing to false positives and other factors, also indicated that DNA effects would be highly unlikely to occur as a result of wear debris from a CoCr implant. Hence, the toxicological weight of evidence suggests that CoCr-containing hip implants are unlikely to be associated with an increased risk of systemic cancers, which is consistent with published and ongoing cancer epidemiology studies involving patients with CoCr hip implants.

SAN: mitigating spatial covariance heterogeneity in cortical thickness data collected from multiple scanners or sites.

In neuroimaging studies, combining data collected from multiple study sites or scanners is becoming common to increase the reproducibility of scientific discoveries. At the same time, unwanted variations arise by using different scanners (inter-scanner biases), which need to be corrected before downstream analyses to facilitate replicable research and prevent spurious findings. While statistical harmonization methods such as ComBat have become popular in mitigating inter-scanner biases in neuroimaging, recent methodological advances have shown that harmonizing heterogeneous covariances results in higher data quality. In vertex-level cortical thickness data, heterogeneity in spatial autocorrelation is a critical factor that affects covariance heterogeneity. Our work proposes a new statistical harmonization method called SAN (Spatial Autocorrelation Normalization) that preserves homogeneous covariance vertex-level cortical thickness data across different scanners. We use an explicit Gaussian process to characterize scanner-invariant and scanner-specific variations to reconstruct spatially homogeneous data across scanners. SAN is computationally feasible, and it easily allows the integration of existing harmonization methods. We demonstrate the utility of the proposed method using cortical thickness data from the Social Processes Initiative in the Neurobiology of the Schizophrenia(s) (SPINS) study. SAN is publicly available as an R package.

Prefrontal cortex structural and developmental associations with callous-unemotional traits and aggression.

Youths with high levels of callous-unemotional (CU) traits and aggression are at an increased risk for developing antisocial behaviours into adulthood. In this population, neurostructural grey matter abnormalities have been observed in the prefrontal cortex. However, the directionality of these associations is inconsistent, prompting some to suggest they may vary across development. Although similar neurodevelopmental patterns have been observed for other disorders featuring emotional and behavioural dysregulation, few studies have tested this hypothesis for CU traits, and particularly not for aggression subtypes. The current study sought to examine grey matter correlates of CU traits and aggression (including its subtypes), and then determine whether these associations varied by age. Fifty-four youths (10-19 years old) who were characterized for CU traits and aggression underwent MRI. Grey matter volume and surface area within the anterior cingulate cortex was positively associated with CU traits. The correlation between CU traits and medial orbitofrontal cortex (mOFC) volume varied significantly as a function of age, as did the correlation between reactive aggression and mOFC surface area. These associations became more positive with age. There were no significant findings for proactive/total aggression. Results are interpreted considering the potential for delayed cortical maturation in youths with high CU traits/aggression.

The impact of acute violent videogame exposure on neurocognitive markers of empathic concern.

Research examining the purported association between violent gaming and aggression remains controversial due to concerns related to methodology, unclear neurocognitive mechanisms, and the failure to adequately consider the role of individual differences in susceptibility. To help address these concerns, we used fMRI and an emotional empathy task to examine whether acute and cumulative violent gaming exposure were associated with abnormalities in emotional empathy as a function of trait-empathy. Emotional empathy was targeted given its involvement in regulating not only aggression, but also other important social functions such as compassion and prosocial behaviour. We hypothesized that violent gaming exposure increases the risk of aberrant social behaviour by altering the aversive value of distress cues. Contrary to expectations, neither behavioural ratings nor empathy-related brain activity varied as a function of violent gaming exposure. Notably, however, activation patterns in somatosensory and motor cortices reflected an interaction between violent gaming exposure and trait empathy. Thus, our results are inconsistent with a straightforward relationship between violent gaming exposure and reduced empathy. Furthermore, they highlight the importance of considering both individual differences in susceptibility and other aspects of cognition related to social functioning to best inform public concern regarding safe gaming practices.

A systematic review of structural and functional magnetic resonance imaging studies on the neurobiology of depressive symptoms in schizophrenia spectrum disorders.

Depressive symptoms in Schizophrenia Spectrum Disorders (SSDs) negatively impact suicidality, prognosis, and quality of life. Despite this, efficacious treatments are limited, largely because the neural mechanisms underlying depressive symptoms in SSDs remain poorly understood. We conducted a systematic review to provide an overview of studies that investigated the neural correlates of depressive symptoms in SSDs using neuroimaging techniques. We searched MEDLINE, PsycINFO, EMBASE, Web of Science, and Cochrane Library databases from inception through June 19, 2023. Specifically, we focused on structural and functional magnetic resonance imaging (MRI), encompassing: (1) T1-weighted imaging measuring brain morphology; (2) diffusion-weighted imaging assessing white matter integrity; or (3) T2*-weighted imaging measures of brain function. Our search yielded 33 articles; 14 structural MRI studies, 18 functional (f)MRI studies, and 1 multimodal fMRI/MRI study. Reviewed studies indicate potential commonalities in the neurobiology of depressive symptoms between SSDs and major depressive disorders, particularly in subcortical and frontal brain regions, though confidence in this interpretation is limited. The review underscores a notable knowledge gap in our understanding of the neurobiology of depression in SSDs, marked by inconsistent approaches and few studies examining imaging metrics of depressive symptoms. Inconsistencies across studies' findings emphasize the necessity for more direct and comprehensive research focusing on the neurobiology of depression in SSDs. Future studies should go beyond "total score" depression metrics and adopt more nuanced assessment approaches considering distinct subdomains. This could reveal unique neurobiological profiles and inform investigations of targeted treatments for depression in SSDs.

Bayesian Optimization of Neurostimulation (BOONStim).

BACKGROUND: Transcranial magnetic stimulation (TMS) treatment response is influenced by individual variability in brain structure and function. Sophisticated, user-friendly approaches, incorporating both established functional magnetic resonance imaging (fMRI) and TMS simulation tools, to identify TMS targets are needed. OBJECTIVE: The current study presents the development and validation of the Bayesian Optimization of Neuro-Stimulation (BOONStim) pipeline. METHODS: BOONStim uses Bayesian optimization for individualized TMS targeting, automating interoperability between surface-based fMRI analytic tools and TMS electric field modeling. Bayesian optimization performance was evaluated in a sample dataset (N=10) using standard circular and functional connectivity-defined targets, and compared to grid optimization. RESULTS: Bayesian optimization converged to similar levels of total electric field stimulation across targets in under 30 iterations, converging within a 5% error of the maxima detected by grid optimization, and requiring less time. CONCLUSIONS: BOONStim is a scalable and configurable user-friendly pipeline for individualized TMS targeting with quick turnaround.

Neural Circuitry and Therapeutic Targeting of Depressive Symptoms in Schizophrenia Spectrum Disorders.

OBJECTIVE: Conceptual similarities between depressive and negative symptoms complicate biomarker and intervention development. This study employed a data-driven approach to delineate the neural circuitry underlying depressive and negative symptoms in schizophrenia spectrum disorders (SSDs). METHODS: Data from three studies were analyzed (157 participants with SSDs) to assess brain-behavior relationships: two neuroimaging studies and a randomized trial of repetitive transcranial magnetic stimulation (rTMS). Partial least squares correlation (PLSC) was used to investigate associations between resting-state functional connectivity and depressive and negative symptoms. Secondary analyses of rTMS trial data (active, N=37; sham, N=33) were used to assess relationships between PLSC-derived symptom profiles and treatment outcomes. RESULTS: PLSC identified three latent variables (LVs) relating functional brain circuitry with symptom profiles. LV1 related a general depressive symptom factor with positive associations between and within the default mode network (DMN), the frontoparietal network (FPN), and the cingulo-opercular network (CON). LV2 related negative symptoms (no depressive symptoms) via negative associations, especially between the FPN and the CON, but also between the DMN and the FPN and the CON. LV3 related a guilt and early wakening depression factor via negative rather than positive associations with the DMN, FPN, and CON. The secondary visual network had a positive association with general depressive symptoms and negative associations with guilt and negative symptoms. Active (but not sham) rTMS applied bilaterally to the dorsolateral prefrontal cortex (DLPFC) reduced general depressive but not guilt-related or negative symptoms. CONCLUSIONS: The results clearly differentiate the neural circuitry underlying depressive and negative symptoms, and segregated across the two-factor structure of depression in SSDs. These findings support divergent neurobiological pathways of depressive symptoms and negative symptoms in people with SSDs. As treatment options are currently limited, bilateral rTMS to the DLPFC is worth exploring further for general depressive symptoms in people with SSDs.

Transdiagnostic Neurobiology of Social Cognition and Individual Variability as Measured by Fractional Amplitude of Low-Frequency Fluctuation in Schizophrenia and Autism Spectrum Disorders.

Fractional amplitude of low-frequency fluctuation (fALFF) is a validated measure of resting-state spontaneous brain activity. Previous fALFF findings in autism and schizophrenia spectrum disorders (ASDs and SSDs) have been highly heterogeneous. We aimed to use fALFF in a large sample of typically developing control (TDC), ASD and SSD participants to explore group differences and relationships with inter-individual variability of fALFF maps and social cognition. fALFF from 495 participants (185 TDC, 68 ASD, and 242 SSD) was computed using functional magnetic resonance imaging as signal power within two frequency bands (i.e., slow-4 and slow-5), normalized by the power in the remaining frequency spectrum. Permutation analysis of linear models was employed to investigate the relationship of fALFF with diagnostic groups, higher-level social cognition, and lower-level social cognition. Each participant's average distance of fALFF map to all others was defined as a variability score, with higher scores indicating less typical maps. Lower fALFF in the visual and higher fALFF in the frontal regions were found in both SSD and ASD participants compared with TDCs. Limited differences were observed between ASD and SSD participants in the cuneus regions only. Associations between slow-4 fALFF and higher-level social cognitive scores across the whole sample were observed in the lateral occipitotemporal and temporoparietal junction. Individual variability within the ASD and SSD groups was also significantly higher compared with TDC. Similar patterns of fALFF and individual variability in ASD and SSD suggest some common neurobiological deficits across these related heterogeneous conditions.

Functional phenotypes in schizophrenia spectrum disorders: defining the constructs and identifying biopsychosocial correlates using data-driven methods.

Functional impairments contribute to poor quality of life in schizophrenia spectrum disorders (SSD). We sought to (Objective I) define the main functional phenotypes in SSD, then (Objective II) identify key biopsychosocial correlates, emphasizing interpretable data-driven methods. Objective I was tested on independent samples: Dataset I (N = 282) and Dataset II (N = 317), with SSD participants who underwent assessment of multiple functioning areas. Participants were clustered based on functioning. Objective II was evaluated in Dataset I by identifying key features for classifying functional phenotype clusters from among 65 sociodemographic, psychological, clinical, cognitive, and brain volume measures. Findings were replicated across latent discriminant analyses (LDA) and one-vs.-rest binomial regularized regressions to identify key predictors. We identified three clusters of participants in each dataset, demonstrating replicable functional phenotypes: Cluster 1-poor functioning across domains; Cluster 2-impaired Role Functioning, but partially preserved Independent and Social Functioning; Cluster 3-good functioning across domains. Key correlates were Avolition, anhedonia, left hippocampal volume, and measures of emotional intelligence and subjective social experience. Avolition appeared more closely tied to role functioning, and anhedonia to independent and social functioning. Thus, we found three replicable functional phenotypes with evidence that recovery may not be uniform across domains. Avolition and anhedonia were both critical but played different roles for different functional domains. It may be important to identify critical functional areas for individual patients and target interventions accordingly.

Multivariate Association between Functional Connectivity Gradients and Cognition in Schizophrenia Spectrum Disorders.

BACKGROUND: Schizophrenia Spectrum Disorders (SSDs), which are characterized by social cognitive deficits, have been associated with dysconnectivity in "unimodal" (e.g., visual, auditory) and "multimodal" (e.g., default-mode and frontoparietal) cortical networks. However, little is known regarding how such dysconnectivity relates to social and non-social cognition, and how such brain-behavioral relationships associate with clinical outcomes of SSDs. METHODS: We analyzed cognitive (non-social and social) measures and resting-state functional magnetic resonance imaging data from the 'Social Processes Initiative in Neurobiology of the Schizophrenia(s) (SPINS)' study (247 stable participants with SSDs and 172 healthy controls, ages 18-55). We extracted gradients from parcellated connectomes and examined the association between the first 3 gradients and the cognitive measures using partial least squares correlation (PLSC). We then correlated the PLSC dimensions with functioning and symptoms in the SSDs group. RESULTS: The SSDs group showed significantly lower differentiation on all three gradients. The first PLSC dimension explained 68.53% (p<.001) of the covariance and showed a significant difference between SSDs and Controls (bootstrap p<.05). PLSC showed that all cognitive measures were associated with gradient scores of unimodal and multimodal networks (Gradient 1), auditory, sensorimotor, and visual networks (Gradient 2), and perceptual networks and striatum (Gradient 3), which were less differentiated in SSDs. Furthermore, the first dimension was positively correlated with negative symptoms and functioning in the SSDs group. CONCLUSIONS: These results suggest a potential role of lower differentiation of brain networks in cognitive and functional impairments in SSDs.

Heterogeneity in functional connectivity: Dimensional predictors of individual variability during rest and task fMRI in psychosis.

BACKGROUND: Individuals with schizophrenia spectrum disorders (SSD) often demonstrate cognitive impairments, associated with poor functional outcomes. While neurobiological heterogeneity has posed challenges when examining social cognition in SSD, it provides a unique opportunity to explore brain-behavior relationships. The aim of this study was to investigate the relationship between individual variability in functional connectivity during resting state and the performance of a social task and social and non-social cognition in a large sample of controls and individuals diagnosed with SSD. METHODS: Neuroimaging and behavioral data were analyzed for 193 individuals with SSD and 155 controls (total n = 348). Individual variability was quantified through mean correlational distance (MCD) of functional connectivity between participants; MCD was defined as a global 'variability score'. Pairwise correlational distance was calculated as 1 - the correlation coefficient between a given pair of participants, and averaging distance from one participant to all other participants provided the mean correlational distance metric. Hierarchical regressions were performed on variability scores derived from resting state and Empathic Accuracy (EA) task functional connectivity data to determine potential predictors (e.g., age, sex, neurocognitive and social cognitive scores) of individual variability. RESULTS: Group comparison between SSD and controls showed greater SSD MCD during rest (p = 0.00038), while no diagnostic differences were observed during task (p = 0.063). Hierarchical regression analyses demonstrated the persistence of a significant diagnostic effect during rest (p = 0.008), contrasting with its non-significance during the task (p = 0.50), after social cognition was added to the model. Notably, social cognition exhibited significance in both resting state and task conditions (both p = 0.01). CONCLUSIONS: Diagnostic differences were more prevalent during unconstrained resting scans, whereas the task pushed participants into a more common pattern which better emphasized transdiagnostic differences in cognitive abilities. Focusing on variability may provide new opportunities for interventions targeting specific cognitive impairments to improve functional outcomes.

Task-based functional neural correlates of social cognition across autism and schizophrenia spectrum disorders.

BACKGROUND: Autism and schizophrenia spectrum disorders (SSDs) both feature atypical social cognition. Despite evidence for comparable group-level performance in lower-level emotion processing and higher-level mentalizing, limited research has examined the neural basis of social cognition across these conditions. Our goal was to compare the neural correlates of social cognition in autism, SSDs, and typically developing controls (TDCs). METHODS: Data came from two harmonized studies in individuals diagnosed with autism or SSDs and TDCs (aged 16-35 years), including behavioral social cognitive metrics and two functional magnetic resonance imaging (fMRI) tasks: a social mirroring Imitate/Observe (ImObs) task and the Empathic Accuracy (EA) task. Group-level comparisons, and transdiagnostic analyses incorporating social cognitive performance, were run using FSL's PALM for each task, covarying for age and sex (1000 permutations, thresholded at p < 0.05 FWE-corrected). Exploratory region of interest (ROI)-based analyses were also conducted. RESULTS: ImObs and EA analyses included 164 and 174 participants, respectively (autism N = 56/59, SSD N = 50/56, TDC N = 58/59). EA and both lower- and higher-level social cognition scores differed across groups. While canonical social cognitive networks were activated, no significant whole-brain or ROI-based group-level differences in neural correlates for either task were detected. Transdiagnostically, neural activity during the EA task, but not the ImObs task, was associated with lower- and higher-level social cognitive performance. LIMITATIONS: Despite attempting to match our groups on age, sex, and race, significant group differences remained. Power to detect regional brain differences is also influenced by sample size and multiple comparisons in whole-brain analyses. Our findings may not generalize to autism and SSD individuals with co-occurring intellectual disabilities. CONCLUSIONS: The lack of whole-brain and ROI-based group-level differences identified and the dimensional EA brain-behavior relationship observed across our sample suggest that the EA task may be well-suited to target engagement in novel intervention testing. Our results also emphasize the potential utility of cross-condition approaches to better understand social cognition across autism and SSDs.

Functional magnetic resonance imaging in schizophrenia: current evidence, methodological advances, limitations and future directions.

Functional neuroimaging emerged with great promise and has provided fundamental insights into the neurobiology of schizophrenia. However, it has faced challenges and criticisms, most notably a lack of clinical translation. This paper provides a comprehensive review and critical summary of the literature on functional neuroimaging, in particular functional magnetic resonance imaging (fMRI), in schizophrenia. We begin by reviewing research on fMRI biomarkers in schizophrenia and the clinical high risk phase through a historical lens, moving from case-control regional brain activation to global connectivity and advanced analytical approaches, and more recent machine learning algorithms to identify predictive neuroimaging features. Findings from fMRI studies of negative symptoms as well as of neurocognitive and social cognitive deficits are then reviewed. Functional neural markers of these symptoms and deficits may represent promising treatment targets in schizophrenia. Next, we summarize fMRI research related to antipsychotic medication, psychotherapy and psychosocial interventions, and neurostimulation, including treatment response and resistance, therapeutic mechanisms, and treatment targeting. We also review the utility of fMRI and data-driven approaches to dissect the heterogeneity of schizophrenia, moving beyond case-control comparisons, as well as methodological considerations and advances, including consortia and precision fMRI. Lastly, limitations and future directions of research in the field are discussed. Our comprehensive review suggests that, in order for fMRI to be clinically useful in the care of patients with schizophrenia, research should address potentially actionable clinical decisions that are routine in schizophrenia treatment, such as which antipsychotic should be prescribed or whether a given patient is likely to have persistent functional impairment. The potential clinical utility of fMRI is influenced by and must be weighed against cost and accessibility factors. Future evaluations of the utility of fMRI in prognostic and treatment response studies may consider including a health economics analysis.

Neuroimaging and Biosample Collection in the Toronto Adolescent and Youth Cohort Study: Rationale, Methods, and Early Data.

BACKGROUND: The Toronto Adolescent and Youth (TAY) Cohort Study will characterize the neurobiological trajectories of psychosis spectrum symptoms, functioning, and suicidality (i.e., suicidal thoughts and behaviors) in youth seeking mental health care. Here, we present the neuroimaging and biosample component of the protocol. We also present feasibility and quality control metrics for the baseline sample collected thus far. METHODS: The current study includes youths (ages 11-24 years) who were referred to child and youth mental health services within a large tertiary care center in Toronto, Ontario, Canada, with target recruitment of 1500 participants. Participants were offered the opportunity to provide any or all of the following: 1) 1-hour magnetic resonance imaging (MRI) scan (electroencephalography if ineligible for or declined MRI), 2) blood sample for genomic and proteomic data (or saliva if blood collection was declined or not feasible) and urine sample, and 3) heart rate recording to assess respiratory sinus arrhythmia. RESULTS: Of the first 417 participants who consented to participate between May 4, 2021, and February 2, 2023, 412 agreed to participate in the imaging and biosample protocol. Of these, 334 completed imaging, 341 provided a biosample, 338 completed respiratory sinus arrhythmia, and 316 completed all 3. Following quality control, data usability was high (MRI: T1-weighted 99%, diffusion-weighted imaging 99%, arterial spin labeling 90%, resting-state functional MRI 95%, task functional MRI 90%; electroencephalography: 83%; respiratory sinus arrhythmia: 99%). CONCLUSIONS: The high consent rates, good completion rates, and high data usability reported here demonstrate the feasibility of collecting and using brain imaging and biosamples in a large clinical cohort of youths seeking mental health care.

RELIEF: A structured multivariate approach for removal of latent inter-scanner effects.

Combining data collected from multiple study sites is becoming common and is advantageous to researchers to increase the generalizability and replicability of scientific discoveries. However, at the same time, unwanted inter-scanner biases are commonly observed across neuroimaging data collected from multiple study sites or scanners, rendering difficulties in integrating such data to obtain reliable findings. While several methods for handling such unwanted variations have been proposed, most of them use univariate approaches that could be too simple to capture all sources of scanner-specific variations. To address these challenges, we propose a novel multivariate harmonization method called RELIEF (REmoval of Latent Inter-scanner Effects through Factorization) for estimating and removing both explicit and latent scanner effects. Our method is the first approach to introduce the simultaneous dimension reduction and factorization of interlinked matrices to a data harmonization context, which provides a new direction in methodological research for correcting inter-scanner biases. Analyzing diffusion tensor imaging (DTI) data from the Social Processes Initiative in Neurobiology of the Schizophrenia (SPINS) study and conducting extensive simulation studies, we show that RELIEF outperforms existing harmonization methods in mitigating inter-scanner biases and retaining biological associations of interest to increase statistical power. RELIEF is publicly available as an R package.