Human brain state dynamics are highly reproducible and associated with neural and behavioral features.

Neural activity and behavior vary within an individual (states) and between individuals (traits). However, the mapping of state-trait neural variation to behavior is not well understood. To address this gap, we quantify moment-to-moment changes in brain-wide co-activation patterns derived from resting-state functional magnetic resonance imaging. In healthy young adults, we identify reproducible spatiotemporal features of co-activation patterns at the single-subject level. We demonstrate that a joint analysis of state-trait neural variations and feature reduction reveal general motifs of individual differences, encompassing state-specific and general neural features that exhibit day-to-day variability. The principal neural variations co-vary with the principal variations of behavioral phenotypes, highlighting cognitive function, emotion regulation, alcohol and substance use. Person-specific probability of occupying a particular co-activation pattern is reproducible and associated with neural and behavioral features. This combined analysis of state-trait variations holds promise for developing reproducible neuroimaging markers of individual life functional outcome.

Ketamine induces multiple individually distinct whole-brain functional connectivity signatures.

Background: Ketamine has emerged as one of the most promising therapies for treatment-resistant depression. However, inter-individual variability in response to ketamine is still not well understood and it is unclear how ketamine's molecular mechanisms connect to its neural and behavioral effects. Methods: We conducted a single-blind placebo-controlled study, with participants blinded to their treatment condition. 40 healthy participants received acute ketamine (initial bolus 0.23 mg/kg, continuous infusion 0.58 mg/kg/hr). We quantified resting-state functional connectivity via data-driven global brain connectivity and related it to individual ketamine-induced symptom variation and cortical gene expression targets. Results: We found that: (i) both the neural and behavioral effects of acute ketamine are multi-dimensional, reflecting robust inter-individual variability; (ii) ketamine's data-driven principal neural gradient effect matched somatostatin (SST) and parvalbumin (PVALB) cortical gene expression patterns in humans, while the mean effect did not; and (iii) behavioral data-driven individual symptom variation mapped onto distinct neural gradients of ketamine, which were resolvable at the single-subject level. Conclusions: These results highlight the importance of considering individual behavioral and neural variation in response to ketamine. They also have implications for the development of individually precise pharmacological biomarkers for treatment selection in psychiatry. Funding: This study was supported by NIH grants DP5OD012109-01 (A.A.), 1U01MH121766 (A.A.), R01MH112746 (J.D.M.), 5R01MH112189 (A.A.), 5R01MH108590 (A.A.), NIAAA grant 2P50AA012870-11 (A.A.); NSF NeuroNex grant 2015276 (J.D.M.); Brain and Behavior Research Foundation Young Investigator Award (A.A.); SFARI Pilot Award (J.D.M., A.A.); Heffter Research Institute (Grant No. 1-190420) (FXV, KHP); Swiss Neuromatrix Foundation (Grant No. 2016-0111) (FXV, KHP); Swiss National Science Foundation under the framework of Neuron Cofund (Grant No. 01EW1908) (KHP); Usona Institute (2015 - 2056) (FXV). Clinical trial number: NCT03842800.

Ketamine is a widely used anesthetic as well as a popular illegal recreational drug. Recently, it has also gained attention as a potential treatment for depression, particularly in cases that don't respond to conventional therapies. However, individuals can vary in their response to ketamine. For example, the drug can alter some people's perception, such as seeing objects as larger or small than they are, while other individuals are unaffected. Although a single dose of ketamine was shown to improve depression symptoms in approximately 65% of patients, the treatment does not work for a significant portion of patients. Understanding why ketamine does not work for everyone could help to identify which patients would benefit most from the treatment. Previous studies investigating ketamine as a treatment for depression have typically included a group of individuals given ketamine and a group given a placebo drug. Assuming people respond similarly to ketamine, the responses in each group were averaged and compared to one another. However, this averaging of results may have masked any individual differences in response to ketamine. As a result, Moujaes et al. set out to investigate whether individuals show differences in brain activity and behavior in response to ketamine. Moujaes et al. monitored the brain activity and behavior of 40 healthy individuals that were first given a placebo drug and then ketamine. The results showed that brain activity and behavior varied significantly between individuals after ketamine administration. Genetic analysis revealed that different gene expression patterns paired with differences in ketamine response in individuals ­ an effect that was hidden when the results were averaged. Ketamine also caused greater differences in brain activity and behavior between individuals than other drugs, such as psychedelics, suggesting ketamine generates a particularly complex response in people. In the future, extending these findings in healthy individuals to those with depression will be crucial for determining whether differences in response to ketamine align with how effective ketamine treatment is for an individual.

On the stability of canonical correlation analysis and partial least squares with application to brain-behavior associations.

Associations between datasets can be discovered through multivariate methods like Canonical Correlation Analysis (CCA) or Partial Least Squares (PLS). A requisite property for interpretability and generalizability of CCA/PLS associations is stability of their feature patterns. However, stability of CCA/PLS in high-dimensional datasets is questionable, as found in empirical characterizations. To study these issues systematically, we developed a generative modeling framework to simulate synthetic datasets. We found that when sample size is relatively small, but comparable to typical studies, CCA/PLS associations are highly unstable and inaccurate; both in their magnitude and importantly in the feature pattern underlying the association. We confirmed these trends across two neuroimaging modalities and in independent datasets with n ≈ 1000 and n = 20,000, and found that only the latter comprised sufficient observations for stable mappings between imaging-derived and behavioral features. We further developed a power calculator to provide sample sizes required for stability and reliability of multivariate analyses. Collectively, we characterize how to limit detrimental effects of overfitting on CCA/PLS stability, and provide recommendations for future studies.

Human brain state dynamics reflect individual neuro-phenotypes.

Neural activity and behavior manifest state and trait dynamics, as well as variation within and between individuals. However, the mapping of state-trait neural variation to behavior is not well understood. To address this gap, we quantify moment-to-moment changes in brain-wide co-activation patterns derived from resting-state functional magnetic resonance imaging. In healthy young adults, we identify reproducible spatio-temporal features of co-activation patterns at the single subject level. We demonstrate that a joint analysis of state-trait neural variations and feature reduction reveal general motifs of individual differences, encompassing state-specific and general neural features that exhibit day-to-day variability. The principal neural variations co-vary with the principal variations of behavioral phenotypes, highlighting cognitive function, emotion regulation, alcohol and substance use. Person-specific probability of occupying a particular co-activation pattern is reproducible and associated with neural and behavioral features. This combined analysis of state-trait variations holds promise for developing reproducible neuroimaging markers of individual life functional outcome.

Robust hierarchically organized whole-brain patterns of dysconnectivity in schizophrenia spectrum disorders observed after personalized intrinsic network topography.

BACKGROUND: Spatial patterns of brain functional connectivity can vary substantially at the individual level. Applying cortical surface-based approaches with individualized rather than group templates may accelerate the discovery of biological markers related to psychiatric disorders. We investigated cortico-subcortical networks from multi-cohort data in people with schizophrenia spectrum disorders (SSDs) and healthy controls (HC) using individualized connectivity profiles. METHODS: We utilized resting-state and anatomical MRI data from n = 406 participants (n = 203 SSD, n = 203 HC) from four cohorts. Functional timeseries were extracted from previously defined intrinsic network subregions of the striatum, thalamus, and cerebellum as well as 80 cortical regions of interest, representing six intrinsic networks using (1) volume-based approaches, (2) a surface-based group atlas approaches, and (3) Personalized Intrinsic Network Topography (PINT). RESULTS: The correlations between all cortical networks and the expected subregions of the striatum, cerebellum, and thalamus were increased using a surface-based approach (Cohen's D volume vs. surface 0.27-1.00, all p < 10-6 ) and further increased after PINT (Cohen's D surface vs. PINT 0.18-0.96, all p < 10-4 ). In SSD versus HC comparisons, we observed robust patterns of dysconnectivity that were strengthened using a surface-based approach and PINT (Number of differing pairwise-correlations: volume: 404, surface: 570, PINT: 628, FDR corrected). CONCLUSION: Surface-based and individualized approaches can more sensitively delineate cortical network dysconnectivity differences in people with SSDs. These robust patterns of dysconnectivity were visibly organized in accordance with the cortical hierarchy, as predicted by computational models.

Functional brain networks reflect spatial and temporal autocorrelation.

High-throughput experimental methods in neuroscience have led to an explosion of techniques for measuring complex interactions and multi-dimensional patterns. However, whether sophisticated measures of emergent phenomena can be traced back to simpler, low-dimensional statistics is largely unknown. To explore this question, we examined resting-state functional magnetic resonance imaging (rs-fMRI) data using complex topology measures from network neuroscience. Here we show that spatial and temporal autocorrelation are reliable statistics that explain numerous measures of network topology. Surrogate time series with subject-matched spatial and temporal autocorrelation capture nearly all reliable individual and regional variation in these topology measures. Network topology changes during aging are driven by spatial autocorrelation, and multiple serotonergic drugs causally induce the same topographic change in temporal autocorrelation. This reductionistic interpretation of widely used complexity measures may help link them to neurobiology.

QuNex-An integrative platform for reproducible neuroimaging analytics.

Introduction: Neuroimaging technology has experienced explosive growth and transformed the study of neural mechanisms across health and disease. However, given the diversity of sophisticated tools for handling neuroimaging data, the field faces challenges in method integration, particularly across multiple modalities and species. Specifically, researchers often have to rely on siloed approaches which limit reproducibility, with idiosyncratic data organization and limited software interoperability. Methods: To address these challenges, we have developed Quantitative Neuroimaging Environment & Toolbox (QuNex), a platform for consistent end-to-end processing and analytics. QuNex provides several novel functionalities for neuroimaging analyses, including a "turnkey" command for the reproducible deployment of custom workflows, from onboarding raw data to generating analytic features. Results: The platform enables interoperable integration of multi-modal, community-developed neuroimaging software through an extension framework with a software development kit (SDK) for seamless integration of community tools. Critically, it supports high-throughput, parallel processing in high-performance compute environments, either locally or in the cloud. Notably, QuNex has successfully processed over 10,000 scans across neuroimaging consortia, including multiple clinical datasets. Moreover, QuNex enables integration of human and non-human workflows via a cohesive translational platform. Discussion: Collectively, this effort stands to significantly impact neuroimaging method integration across acquisition approaches, pipelines, datasets, computational environments, and species. Building on this platform will enable more rapid, scalable, and reproducible impact of neuroimaging technology across health and disease.

Informing individualized multi-scale neural signatures of clozapine response in patients with treatment-refractory schizophrenia.

Clozapine is currently the only antipsychotic with demonstrated efficacy in treatment-refractory schizophrenia (TRS). However, response to clozapine differs widely between TRS patients, and there are no available clinical or neural predictive indicators that could be used to increase or accelerate the use of clozapine in patients who stand to benefit. Furthermore, it remains unclear how the neuropharmacology of clozapine contributes to its therapeutic effects. Identifying the mechanisms underlying clozapine's therapeutic effects across domains of symptomatology could be crucial for development of new optimized therapies for TRS. Here, we present results from a prospective neuroimaging study that quantitatively related heterogeneous patterns of clinical clozapine response to neural functional connectivity at baseline. We show that we can reliably capture specific dimensions of clozapine clinical response by quantifying the full variation across item-level clinical scales, and that these dimensions can be mapped to neural features that are sensitive to clozapine-induced symptom change. Thus, these features may act as "failure modes" that can provide an early indication of treatment (non-)responsiveness. Lastly, we related the response-relevant neural maps to spatial expression profiles of genes coding for receptors implicated in clozapine's pharmacology, demonstrating that distinct dimensions of clozapine symptom-informed neural features may be associated with specific receptor targets. Collectively, this study informs prognostic neuro-behavioral measures for clozapine as a more optimal treatment for selected patients with TRS. We provide support for the identification of neuro-behavioral targets linked to pharmacological efficacy that can be further developed to inform optimal early treatment decisions in schizophrenia.

Toward Mapping Neurobehavioral Heterogeneity of Psychedelic Neurobiology in Humans.

Precision psychiatry aims to identify markers of interindividual variability that allow for predicting the right treatment for each patient. However, bridging the gap between molecular-level manipulations and neural systems-level functional alterations remains an unsolved problem in psychiatry. After decades of low success rates in pharmaceutical research and development for psychiatric drugs, multiple studies now point to the potential of psychedelics as a promising, fast-acting, and long-lasting treatment for some psychiatric symptoms. Yet, given the highly psychoactive nature of these substances, a precision medicine approach is essential to map the neural signals related to clinical efficacy to identify patients who can maximally benefit from this treatment. Recent studies have shown that bridging the gap between pharmacology, systems-level neural response in humans, and individual experience is possible for psychedelic substances, therefore paving the way for a precision neuropsychiatric therapeutic development. Specifically, it has been shown that the integration of brain-wide positron emission tomography or transcriptomic data, i.e., receptor distribution for the serotonin 2A receptor, with computational neuroimaging methods can simulate the effect of psychedelics on the human brain. These novel computational psychiatry approaches allow for modeling interindividual differences in neural as well as subjective effects of psychedelic substances. Collectively, this review provides a deep dive into psychedelic pharmaconeuroimaging studies with a core focus on how recent computational psychiatry advances in biophysically based circuit modeling can be leveraged to predict individual responses. Finally, we emphasize the importance of human pharmacological neuroimaging for the continued precision therapeutic development of psychedelics.

Accelerated cortical thinning precedes and predicts conversion to psychosis: The NAPLS3 longitudinal study of youth at clinical high-risk.

Progressive grey matter loss has been demonstrated among clinical high-risk (CHR) individuals who convert to psychosis, but it is unknown whether these changes occur prior to psychosis onset. Identifying illness-related neurobiological mechanisms that occur prior to conversion is essential for targeted early intervention. Among participants in the third wave of the North American Prodrome Longitudinal Study (NAPLS3), this report investigated if steeper cortical thinning was observable prior to psychosis onset among CHR individuals who ultimately converted (CHR-C) and assessed the shortest possible time interval in which rates of cortical thinning differ between CHR-C, CHR non-converters (CHR-NC), and health controls (HC). 338 CHR-NC, 42 CHR-C, and 62 HC participants (age 19.3±4.2, 44.8% female, 52.5% racial/ethnic minority) completed up to 5 MRI scans across 8 months. Accelerated thinning among CHR-C compared to CHR-NC and HC was observed in multiple prefrontal, temporal, and parietal cortical regions. CHR-NC also exhibited accelerated cortical thinning compared to HC in several of these areas. Greater percent decrease in cortical thickness was observed among CHR-C compared to other groups across 2.9±1.8 months, on average, in several cortical areas. ROC analyses discriminating CHR-C from CHR-NC by percent thickness change in a left hemisphere region of interest, scanner, age, age2, and sex had an AUC of 0.74, with model predictive power driven primarily by percent thickness change. Findings indicate that accelerated cortical thinning precedes psychosis onset and differentiates CHR-C from CHR-NC and HC across short time intervals. Mechanisms underlying cortical thinning may provide novel treatment targets prior to psychosis onset.

A common symptom geometry of mood improvement under sertraline and placebo associated with distinct neural patterns.

Importance: Understanding the mechanisms of major depressive disorder (MDD) improvement is a key challenge to determine effective personalized treatments. Objective: To perform a secondary analysis quantifying neural-to-symptom relationships in MDD as a function of antidepressant treatment. Design: Double blind randomized controlled trial. Setting: Multicenter. Participants: Patients with early onset recurrent depression from the public Establishing Moderators and Biosignatures of Antidepressant Response in Clinical Care (EMBARC) study. Interventions: Either sertraline or placebo during 8 weeks (stage 1), and according to response a second line of treatment for 8 additional weeks (stage 2). Main Outcomes and Measures: To identify a data-driven pattern of symptom variations during these two stages, we performed a Principal Component Analysis (PCA) on the variations of individual items of four clinical scales measuring depression, anxiety, suicidal ideas and manic-like symptoms, resulting in a univariate measure of clinical improvement. We then investigated how initial clinical and neural factors predicted this measure during stage 1. To do so, we extracted resting-state global brain connectivity (GBC) at baseline at the individual level using a whole-brain functional network parcellation. In turn, we computed a linear model for each brain parcel with individual data-driven clinical improvement scores during stage 1 for each group. Results: 192 patients (127 women), age 37.7 years old (standard deviation: 13.5), were included. The first PC (PC1) capturing 20% of clinical variation was similar across treatment groups at stage 1 and stage 2, suggesting a reproducible pattern of symptom improvement. PC1 patients' scores significantly differed according to treatment during stage 1, whereas no difference of response was evidenced between groups with the Clinical Global Impressions (CGI). Baseline GBC correlated to stage 1 PC1 scores in the sertraline, but not in the placebo group. Conclusions and Relevance: Using data-driven reduction of symptoms scales, we identified a common profile of symptom improvement across placebo and sertraline. However, the neural patterns of baseline that mapped onto symptom improvement distinguished between treatment and placebo. Our results underscore that mapping from data-driven symptom improvement onto neural circuits is vital to detect treatment-responsive neural profiles that may aid in optimal patient selection for future trials.

Reward and loss incentives improve spatial working memory by shaping trial-by-trial posterior frontoparietal signals.

Integrating motivational signals with cognition is critical for goal-directed activities. The mechanisms that link neural changes with motivated working memory continue to be understood. Here, we tested how externally cued and non-cued (internally represented) reward and loss impact spatial working memory precision and neural circuits in human subjects using fMRI. We translated the classic delayed-response spatial working memory paradigm from non-human primate studies to take advantage of a continuous numeric measure of working memory precision, and the wealth of translational neuroscience yielded by these studies. Our results demonstrated that both cued and non-cued reward and loss improved spatial working memory precision. Visual association regions of the posterior prefrontal and parietal cortices, specifically the precentral sulcus (PCS) and intraparietal sulcus (IPS), had increased BOLD signal during incentivized spatial working memory. A subset of these regions had trial-by-trial increases in BOLD signal that were associated with better working memory precision, suggesting that these regions may be critical for linking neural signals with motivated working memory. In contrast, regions straddling executive networks, including areas in the dorsolateral prefrontal cortex, anterior parietal cortex and cerebellum displayed decreased BOLD signal during incentivized working memory. While reward and loss similarly impacted working memory processes, they dissociated during feedback when money won or avoided in loss was given based on working memory performance. During feedback, the trial-by-trial amount and valence of reward/loss received was dissociated amongst regions such as the ventral striatum, habenula and periaqueductal gray. Overall, this work suggests motivated spatial working memory is supported by complex sensory processes, and that the IPS and PCS in the posterior frontoparietal cortices may be key regions for integrating motivational signals with spatial working memory precision.

Computational Modeling of Electroencephalography and Functional Magnetic Resonance Imaging Paradigms Indicates a Consistent Loss of Pyramidal Cell Synaptic Gain in Schizophrenia.

BACKGROUND: Diminished synaptic gain-the sensitivity of postsynaptic responses to neural inputs-may be a fundamental synaptic pathology in schizophrenia. Evidence for this is indirect, however. Furthermore, it is unclear whether pyramidal cells or interneurons (or both) are affected, or how these deficits relate to symptoms. METHODS: People with schizophrenia diagnoses (PScz) (n = 108), their relatives (n = 57), and control subjects (n = 107) underwent 3 electroencephalography (EEG) paradigms-resting, mismatch negativity, and 40-Hz auditory steady-state response-and resting functional magnetic resonance imaging. Dynamic causal modeling was used to quantify synaptic connectivity in cortical microcircuits. RESULTS: Classic group differences in EEG features between PScz and control subjects were replicated, including increased theta and other spectral changes (resting EEG), reduced mismatch negativity, and reduced 40-Hz power. Across all 4 paradigms, characteristic PScz data features were all best explained by models with greater self-inhibition (decreased synaptic gain) in pyramidal cells. Furthermore, disinhibition in auditory areas predicted abnormal auditory perception (and positive symptoms) in PScz in 3 paradigms. CONCLUSIONS: First, characteristic EEG changes in PScz in 3 classic paradigms are all attributable to the same underlying parameter change: greater self-inhibition in pyramidal cells. Second, psychotic symptoms in PScz relate to disinhibition in neural circuits. These findings are more commensurate with the hypothesis that in PScz, a primary loss of synaptic gain on pyramidal cells is then compensated by interneuron downregulation (rather than the converse). They further suggest that psychotic symptoms relate to this secondary downregulation.

Transcriptomics-informed large-scale cortical model captures topography of pharmacological neuroimaging effects of LSD.

Psychoactive drugs can transiently perturb brain physiology while preserving brain structure. The role of physiological state in shaping neural function can therefore be investigated through neuroimaging of pharmacologically induced effects. Previously, using pharmacological neuroimaging, we found that neural and experiential effects of lysergic acid diethylamide (LSD) are attributable to agonism of the serotonin-2A receptor (Preller et al., 2018). Here, we integrate brain-wide transcriptomics with biophysically based circuit modeling to simulate acute neuromodulatory effects of LSD on human cortical large-scale spatiotemporal dynamics. Our model captures the inter-areal topography of LSD-induced changes in cortical blood oxygen level-dependent (BOLD) functional connectivity. These findings suggest that serotonin-2A-mediated modulation of pyramidal-neuronal gain is a circuit mechanism through which LSD alters cortical functional topography. Individual-subject model fitting captures patterns of individual neural differences in pharmacological response related to altered states of consciousness. This work establishes a framework for linking molecular-level manipulations to systems-level functional alterations, with implications for precision medicine.

Mapping brain-behavior space relationships along the psychosis spectrum.

Difficulties in advancing effective patient-specific therapies for psychiatric disorders highlight a need to develop a stable neurobiologically grounded mapping between neural and symptom variation. This gap is particularly acute for psychosis-spectrum disorders (PSD). Here, in a sample of 436 PSD patients spanning several diagnoses, we derived and replicated a dimensionality-reduced symptom space across hallmark psychopathology symptoms and cognitive deficits. In turn, these symptom axes mapped onto distinct, reproducible brain maps. Critically, we found that multivariate brain-behavior mapping techniques (e.g. canonical correlation analysis) do not produce stable results with current sample sizes. However, we show that a univariate brain-behavioral space (BBS) can resolve stable individualized prediction. Finally, we show a proof-of-principle framework for relating personalized BBS metrics with molecular targets via serotonin and glutamate receptor manipulations and neural gene expression maps derived from the Allen Human Brain Atlas. Collectively, these results highlight a stable and data-driven BBS mapping across PSD, which offers an actionable path that can be iteratively optimized for personalized clinical biomarker endpoints.

White matter changes in psychosis risk relate to development and are not impacted by the transition to psychosis.

Subtle alterations in white matter microstructure are observed in youth at clinical high risk (CHR) for psychosis. However, the timing of these changes and their relationships to the emergence of psychosis remain unclear. Here, we track the evolution of white matter abnormalities in a large, longitudinal cohort of CHR individuals comprising the North American Prodrome Longitudinal Study (NAPLS-3). Multi-shell diffusion magnetic resonance imaging data were collected across multiple timepoints (1-5 over 1 year) in 286 subjects (aged 12-32 years): 25 CHR individuals who transitioned to psychosis (CHR-P; 61 scans), 205 CHR subjects with unknown transition outcome after the 1-year follow-up period (CHR-U; 596 scans), and 56 healthy controls (195 scans). Linear mixed effects models were fitted to infer the impact of age and illness-onset on variation in the fractional anisotropy of cellular tissue (FAT) and the volume fraction of extracellular free water (FW). Baseline measures of white matter microstructure did not differentiate between HC, CHR-U and CHR-P individuals. However, age trajectories differed between the three groups in line with a developmental effect: CHR-P and CHR-U groups displayed higher FAT in adolescence, and 4% lower FAT by 30 years of age compared to controls. Furthermore, older CHR-P subjects (20+ years) displayed 4% higher FW in the forceps major (p < 0.05). Prospective analysis in CHR-P did not reveal a significant impact of illness onset on regional FAT or FW, suggesting that transition to psychosis is not marked by dramatic change in white matter microstructure. Instead, clinical high risk for psychosis-regardless of transition outcome-is characterized by subtle age-related white matter changes that occur in tandem with development.

A Whole-Brain and Cross-Diagnostic Perspective on Functional Brain Network Dysfunction.

A wide variety of mental disorders have been associated with resting-state functional network alterations, which are thought to contribute to the cognitive changes underlying mental illness. These observations appear to support theories postulating large-scale disruptions of brain systems in mental illness. However, existing approaches isolate differences in network organization without putting those differences in a broad, whole-brain perspective. Using a graph distance approach-connectome-wide similarity-we found that whole-brain resting-state functional network organization is highly similar across groups of individuals with and without a variety of mental diseases. This similarity was observed across autism spectrum disorder, attention-deficit hyperactivity disorder, and schizophrenia. Nonetheless, subtle differences in network graph distance were predictive of diagnosis, suggesting that while functional connectomes differ little across health and disease, those differences are informative. These results suggest a need to reevaluate neurocognitive theories of mental illness, with a role for subtle functional brain network changes in the production of an array of mental diseases. Such small network alterations suggest the possibility that small, well-targeted alterations to brain network organization may provide meaningful improvements for a variety of mental disorders.

Psilocybin Induces Time-Dependent Changes in Global Functional Connectivity.

BACKGROUND: The use of psilocybin in scientific and experimental clinical contexts has triggered renewed interest in the mechanism of action of psychedelics. However, its time-dependent systems-level neurobiology remains sparsely investigated in humans. METHODS: We conducted a double-blind, randomized, counterbalanced, crossover study comprising 23 healthy human participants who received placebo and 0.2 mg/kg of psilocybin orally on 2 different test days. Participants underwent magnetic resonance imaging at 3 time points between administration and peak effects: 20 minutes, 40 minutes, and 70 minutes after administration. Resting-state functional connectivity was quantified via a data-driven global brain connectivity method and compared with cortical gene expression maps. RESULTS: Psilocybin reduced associative, but concurrently increased sensory, brain-wide connectivity. This pattern emerged over time from administration to peak effects. Furthermore, we showed that baseline connectivity is associated with the extent of psilocybin-induced changes in functional connectivity. Lastly, psilocybin-induced changes correlated in a time-dependent manner with spatial gene expression patterns of the 5-HT2A (5-hydroxytryptamine 2A) and 5-HT1A (5-hydroxytryptamine 1A) receptors. CONCLUSIONS: These results suggest that the integration of functional connectivity in sensory regions and the disintegration in associative regions may underlie the psychedelic state and pinpoint the critical role of the serotonin 2A and 1A receptor systems. Furthermore, baseline connectivity may represent a predictive marker of the magnitude of changes induced by psilocybin and may therefore contribute to a personalized medicine approach within the potential framework of psychedelic treatment.