Individualized epidemic spreading models predict epilepsy surgery outcomes: A pseudo-prospective study.

Epilepsy surgery is the treatment of choice for drug-resistant epilepsy patients, but up to 50% of patients continue to have seizures one year after the resection. In order to aid presurgical planning and predict postsurgical outcome on a patient-by-patient basis, we developed a framework of individualized computational models that combines epidemic spreading with patient-specific connectivity and epileptogeneity maps: the Epidemic Spreading Seizure and Epilepsy Surgery framework (ESSES). ESSES parameters were fitted in a retrospective study (N = 15) to reproduce invasive electroencephalography (iEEG)-recorded seizures. ESSES reproduced the iEEG-recorded seizures, and significantly better so for patients with good (seizure-free, SF) than bad (nonseizure-free, NSF) outcome. We illustrate here the clinical applicability of ESSES with a pseudo-prospective study (N = 34) with a blind setting (to the resection strategy and surgical outcome) that emulated presurgical conditions. By setting the model parameters in the retrospective study, ESSES could be applied also to patients without iEEG data. ESSES could predict the chances of good outcome after any resection by finding patient-specific model-based optimal resection strategies, which we found to be smaller for SF than NSF patients, suggesting an intrinsic difference in the network organization or presurgical evaluation results of NSF patients. The actual surgical plan overlapped more with the model-based optimal resection, and had a larger effect in decreasing modeled seizure propagation, for SF patients than for NSF patients. Overall, ESSES could correctly predict 75% of NSF and 80.8% of SF cases pseudo-prospectively. Our results show that individualised computational models may inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection. This is the first time that such a model is validated with a fully independent cohort and without the need for iEEG recordings.

Individualized computational models of epilepsy surgery capture some of the key aspects of seizure propagation and the resective surgery. It is to be established whether this information can be integrated during the presurgical evaluation of the patient to improve surgical planning and the chances of a good surgical outcome. Here we address this question with a pseudo-prospective study that applies a computational framework of seizure propagation and epilepsy surgery­the ESSES framework­in a pseudo-prospective study mimicking the presurgical conditions. We found that within this pseudo-prospective setting, ESSES could correctly predict 75% of NSF and 80.8% of SF cases. This finding suggests the potential of individualised computational models to inform surgical planning by suggesting alternative resections and providing information on the likelihood of a good outcome after a proposed resection.

Inducing a meditative state by artificial perturbations: A mechanistic understanding of brain dynamics underlying meditation.

Contemplative neuroscience has increasingly explored meditation using neuroimaging. However, the brain mechanisms underlying meditation remain elusive. Here, we implemented a mechanistic framework to explore the spatiotemporal dynamics of expert meditators during meditation and rest, and controls during rest. We first applied a model-free approach by defining a probabilistic metastable substate (PMS) space for each condition, consisting of different probabilities of occurrence from a repertoire of dynamic patterns. Moreover, we implemented a model-based approach by adjusting the PMS of each condition to a whole-brain model, which enabled us to explore in silico perturbations to transition from resting-state to meditation and vice versa. Consequently, we assessed the sensitivity of different brain areas regarding their perturbability and their mechanistic local-global effects. Overall, our work reveals distinct whole-brain dynamics in meditation compared to rest, and how transitions can be induced with localized artificial perturbations. It motivates future work regarding meditation as a practice in health and as a potential therapy for brain disorders.

Our work explores brain dynamics in a group of expert meditators and controls. First, we characterized meditation and rest with a repertoire of brain patterns, each with its distinct probability of occurrence. Then, we generated whole-brain models of each condition, which enabled us to artificially perturb the systems to induce transitions between rest and meditation. Our results open new avenues in meditation research as a practice in health and disease.

In-vivo neuronal dysfunction by Aß and tau overlaps with brain-wide inflammatory mechanisms in Alzheimer's disease.

The molecular mechanisms underlying neuronal dysfunction in Alzheimer's disease (AD) remain uncharacterized. Here, we identify genes, molecular pathways and cellular components associated with whole-brain dysregulation caused by amyloid-beta (Aß) and tau deposits in the living human brain. We obtained in-vivo resting-state functional MRI (rs-fMRI), Aß- and tau-PET for 47 cognitively unimpaired and 16 AD participants from the Translational Biomarkers in Aging and Dementia cohort. Adverse neuronal activity impacts by Aß and tau were quantified with personalized dynamical models by fitting pathology-mediated computational signals to the participant's real rs-fMRIs. Then, we detected robust brain-wide associations between the spatial profiles of Aß-tau impacts and gene expression in the neurotypical transcriptome (Allen Human Brain Atlas). Within the obtained distinctive signature of in-vivo neuronal dysfunction, several genes have prominent roles in microglial activation and in interactions with Aß and tau. Moreover, cellular vulnerability estimations revealed strong association of microglial expression patterns with Aß and tau's synergistic impact on neuronal activity (q < 0.001). These results further support the central role of the immune system and neuroinflammatory pathways in AD pathogenesis. Neuronal dysregulation by AD pathologies also associated with neurotypical synaptic and developmental processes. In addition, we identified drug candidates from the vast LINCS library to halt or reduce the observed Aß-tau effects on neuronal activity. Top-ranked pharmacological interventions target inflammatory, cancer and cardiovascular pathways, including specific medications undergoing clinical evaluation in AD. Our findings, based on the examination of molecular-pathological-functional interactions in humans, may accelerate the process of bringing effective therapies into clinical practice.

Altered Structural Connectivity and Functional Brain Dynamics in Individuals With Heavy Alcohol Use Elucidated via Network Control Theory.

BACKGROUND: Heavy alcohol use and its associated conditions, such as alcohol use disorder, impact millions of individuals worldwide. While our understanding of the neurobiological correlates of alcohol use has evolved substantially, we still lack models that incorporate whole-brain neuroanatomical, functional, and pharmacological information under one framework. METHODS: Here, we utilized diffusion and functional magnetic resonance imaging to investigate alterations to brain dynamics in 130 individuals with a high amount of current alcohol use. We compared these alcohol-using individuals to 308 individuals with minimal use of any substances. RESULTS: We found that individuals with heavy alcohol use had less dynamic and complex brain activity, and through leveraging network control theory, had increased control energy to complete transitions between activation states. Furthermore, using separately acquired positron emission tomography data, we deployed an in silico evaluation demonstrating that decreased D2 receptor levels, as found previously in individuals with alcohol use disorder, may relate to our observed findings. CONCLUSIONS: This work demonstrates that whole-brain, multimodal imaging information can be combined under a network control framework to identify and evaluate neurobiological correlates and mechanisms of heavy alcohol use.

Gaming expertise induces meso­scale brain plasticity and efficiency mechanisms as revealed by whole-brain modeling.

Video games are a valuable tool for studying the effects of training and neural plasticity on the brain. However, the underlying mechanisms related to plasticity-associated brain structural changes and their impact on brain dynamics are unknown. Here, we used a semi-empirical whole-brain model to study structural neural plasticity mechanisms linked to video game expertise. We hypothesized that video game expertise is associated with neural plasticity-mediated changes in structural connectivity that manifest at the meso­scale level, resulting in a more segregated functional network topology. To test this hypothesis, we combined structural connectivity data of StarCraft II video game players (VGPs, n = 31) and non-players (NVGPs, n = 31), with generic fMRI data from the Human Connectome Project and computational models, to generate simulated fMRI recordings. Graph theory analysis on simulated data was performed during both resting-state conditions and external stimulation. VGPs' simulated functional connectivity was characterized by a meso­scale integration, with increased local connectivity in frontal, parietal, and occipital brain regions. The same analyses at the level of structural connectivity showed no differences between VGPs and NVGPs. Regions that increased their connectivity strength in VGPs are known to be involved in cognitive processes crucial for task performance such as attention, reasoning, and inference. In-silico stimulation suggested that differences in FC between VGPs and NVGPs emerge in noisy contexts, specifically when the noisy level of stimulation is increased. This indicates that the connectomes of VGPs may facilitate the filtering of noise from stimuli. These structural alterations drive the meso­scale functional changes observed in individuals with gaming expertise. Overall, our work sheds light on the mechanisms underlying structural neural plasticity triggered by video game experiences.

Turbulent dynamics and whole-brain modeling: toward new clinical applications for traumatic brain injury.

Traumatic Brain Injury (TBI) is a prevalent disorder mostly characterized by persistent impairments in cognitive function that poses a substantial burden on caregivers and the healthcare system worldwide. Crucially, severity classification is primarily based on clinical evaluations, which are non-specific and poorly predictive of long-term disability. In this Mini Review, we first provide a description of our model-free and model-based approaches within the turbulent dynamics framework as well as our vision on how they can potentially contribute to provide new neuroimaging biomarkers for TBI. In addition, we report the main findings of our recent study examining longitudinal changes in moderate-severe TBI (msTBI) patients during a one year spontaneous recovery by applying the turbulent dynamics framework (model-free approach) and the Hopf whole-brain computational model (model-based approach) combined with in silico perturbations. Given the neuroinflammatory response and heightened risk for neurodegeneration after TBI, we also offer future directions to explore the association with genomic information. Moreover, we discuss how whole-brain computational modeling may advance our understanding of the impact of structural disconnection on whole-brain dynamics after msTBI in light of our recent findings. Lastly, we suggest future avenues whereby whole-brain computational modeling may assist the identification of optimal brain targets for deep brain stimulation to promote TBI recovery.

The evolution of whole-brain turbulent dynamics during recovery from traumatic brain injury.

It has been previously shown that traumatic brain injury (TBI) is associated with reductions in metastability in large-scale networks in resting-state fMRI (rsfMRI). However, little is known about how TBI affects the local level of synchronization and how this evolves during the recovery trajectory. Here, we applied a novel turbulent dynamics framework to investigate whole-brain dynamics using an rsfMRI dataset from a cohort of moderate to severe TBI patients and healthy controls (HCs). We first examined how several measures related to turbulent dynamics differ between HCs and TBI patients at 3, 6, and 12 months post-injury. We found a significant reduction in these empirical measures after TBI, with the largest change at 6 months post-injury. Next, we built a Hopf whole-brain model with coupled oscillators and conducted in silico perturbations to investigate the mechanistic principles underlying the reduced turbulent dynamics found in the empirical data. A simulated attack was used to account for the effect of focal lesions. This revealed a shift to lower coupling parameters in the TBI dataset and, critically, decreased susceptibility and information-encoding capability. These findings confirm the potential of the turbulent framework to characterize longitudinal changes in whole-brain dynamics and in the reactivity to external perturbations after TBI.

Viscous dynamics associated with hypoexcitation and structural disintegration in neurodegeneration via generative whole-brain modeling.

INTRODUCTION: Alzheimer's disease (AD) and behavioral variant frontotemporal dementia (bvFTD) lack mechanistic biophysical modeling in diverse, underrepresented populations. Electroencephalography (EEG) is a high temporal resolution, cost-effective technique for studying dementia globally, but lacks mechanistic models and produces non-replicable results. METHODS: We developed a generative whole-brain model that combines EEG source-level metaconnectivity, anatomical priors, and a perturbational approach. This model was applied to Global South participants (AD, bvFTD, and healthy controls). RESULTS: Metaconnectivity outperformed pairwise connectivity and revealed more viscous dynamics in patients, with altered metaconnectivity patterns associated with multimodal disease presentation. The biophysical model showed that connectome disintegration and hypoexcitability triggered altered metaconnectivity dynamics and identified critical regions for brain stimulation. We replicated the main results in a second subset of participants for validation with unharmonized, heterogeneous recording settings. DISCUSSION: The results provide a novel agenda for developing mechanistic model-inspired characterization and therapies in clinical, translational, and computational neuroscience settings.

Metastability indexes global changes in the dynamic working point of the brain following brain stimulation.

Several studies have shown that coordination among neural ensembles is a key to understand human cognition. A well charted path is to identify coordination states associated with cognitive functions from spectral changes in the oscillations of EEG or MEG. A growing number of studies suggest that the tendency to switch between coordination states, sculpts the dynamic repertoire of the brain and can be indexed by a measure known as metastability. In this article, we characterize perturbations in the metastability of global brain network dynamics following Transcranial Magnetic Stimulation that could quantify the duration for which information processing is altered. Thus allowing researchers to understand the network effects of brain stimulation, standardize stimulation protocols and design experimental tasks. We demonstrate the effect empirically using publicly available datasets and use a digital twin (a whole brain connectome model) to understand the dynamic principles that generate such observations. We observed a significant reduction in metastability, concurrent with an increase in coherence following single-pulse TMS reflecting the existence of a window where neural coordination is altered. The reduction in complexity was validated by an additional measure based on the Lempel-Ziv complexity of microstate labeled EEG data. Interestingly, higher frequencies in the EEG signal showed faster recovery in metastability than lower frequencies. The digital twin shed light on how the phase resetting introduced by the single-pulse TMS in local cortical networks can propagate globally, giving rise to changes in metastability and coherence.

Low-dimensional organization of global brain states of reduced consciousness.

Brain states are frequently represented using a unidimensional scale measuring the richness of subjective experience (level of consciousness). This description assumes a mapping between the high-dimensional space of whole-brain configurations and the trajectories of brain states associated with changes in consciousness, yet this mapping and its properties remain unclear. We combine whole-brain modeling, data augmentation, and deep learning for dimensionality reduction to determine a mapping representing states of consciousness in a low-dimensional space, where distances parallel similarities between states. An orderly trajectory from wakefulness to patients with brain injury is revealed in a latent space whose coordinates represent metrics related to functional modularity and structure-function coupling, increasing alongside loss of consciousness. Finally, we investigate the effects of model perturbations, providing geometrical interpretation for the stability and reversibility of states. We conclude that conscious awareness depends on functional patterns encoded as a low-dimensional trajectory within the vast space of brain configurations.

Personalized whole brain modeling of status epilepticus.

How status epilepticus (SE) is generated and propagates in the brain is not known. As for seizures, a patient-specific approach is necessary, and the analysis should be performed at the whole brain level. Personalized brain models can be used to study seizure genesis and propagation at the whole brain scale in The Virtual Brain (TVB), using the Epileptor mathematical construct. Building on the fact that SE is part of the repertoire of activities that the Epileptor can generate, we present the first attempt to model SE at the whole brain scale in TVB, using data from a patient who experienced SE during presurgical evaluation. Simulations reproduced the patterns found with SEEG recordings. We find that if, as expected, the pattern of SE propagation correlates with the properties of the patient's structural connectome, SE propagation also depends upon the global state of the network, i.e., that SE propagation is an emergent property. We conclude that individual brain virtualization can be used to study SE genesis and propagation. This type of theoretical approach may be used to design novel interventional approaches to stop SE. This paper was presented at the 8th London-Innsbruck Colloquium on Status Epilepticus and Acute Seizures held in September 2022.

Improving the study of brain-behavior relationships by revisiting basic assumptions.

Neuroimaging research has been at the forefront of concerns regarding the failure of experimental findings to replicate. In the study of brain-behavior relationships, past failures to find replicable and robust effects have been attributed to methodological shortcomings. Methodological rigor is important, but there are other overlooked possibilities: most published studies share three foundational assumptions, often implicitly, that may be faulty. In this paper, we consider the empirical evidence from human brain imaging and the study of non-human animals that calls each foundational assumption into question. We then consider the opportunities for a robust science of brain-behavior relationships that await if scientists ground their research efforts in revised assumptions supported by current empirical evidence.

Dissection of a sensorimotor circuit underlying pathogen aversion in C. elegans.

BACKGROUND: Altering animal behavior to reduce pathogen exposure is a key line of defense against pathogen attack. In Caenorhabditis elegans, alterations in intestinal physiology caused by pathogen colonization and sensation of microbial metabolites may lead to activation of pathogen aversive behaviors ranging from aversive reflexes to learned avoidance. However, the neural circuitry between chemosensory neurons that sense pathogenic bacterial cues and the motor neurons responsible for avoidance-associated locomotion remains unknown. RESULTS: Using C. elegans, we found that backward locomotion was a component of learned pathogen avoidance, as animals pre-exposed to Pseudomonas aeruginosa or Enterococcus faecalis showed reflexive aversion to drops of the bacteria driven by chemosensory neurons, including the olfactory AWB neurons. This response also involved intestinal distention and, for E. faecalis, required expression of TRPM channels in the intestine and excretory system. Additionally, we uncovered a circuit composed of olfactory neurons, interneurons, and motor neurons that controls the backward locomotion crucial for learned reflexive aversion to pathogenic bacteria, learned avoidance, and the repulsive odor 2-nonanone. CONCLUSIONS: Using whole-brain simulation and functional assays, we uncovered a novel sensorimotor circuit governing learned reflexive aversion. The discovery of a complete sensorimotor circuit for reflexive aversion demonstrates the utility of using the C. elegans connectome and computational modeling in uncovering new neuronal regulators of behavior.

Meditation-induced effects on whole-brain structural and effective connectivity.

In the past decades, there has been a growing scientific interest in characterizing neural correlates of meditation training. Nonetheless, the mechanisms underlying meditation remain elusive. In the present work, we investigated meditation-related changes in functional dynamics and structural connectivity (SC). For this purpose, we scanned experienced meditators and control (naive) subjects using magnetic resonance imaging (MRI) to acquire structural and functional data during two conditions, resting-state and meditation (focused attention on breathing). In this way, we aimed to characterize and distinguish both short-term and long-term modifications in the brain's structure and function. First, to analyze the fMRI data, we calculated whole-brain effective connectivity (EC) estimates, relying on a dynamical network model to replicate BOLD signals' spatio-temporal structure, akin to functional connectivity (FC) with lagged correlations. We compared the estimated EC, FC, and SC links as features to train classifiers to predict behavioral conditions and group identity. Then, we performed a network-based analysis of anatomical connectivity. We demonstrated through a machine-learning approach that EC features were more informative than FC and SC solely. We showed that the most informative EC links that discriminated between meditators and controls involved several large-scale networks mainly within the left hemisphere. Moreover, we found that differences in the functional domain were reflected to a smaller extent in changes at the anatomical level as well. The network-based analysis of anatomical pathways revealed strengthened connectivity for meditators compared to controls between four areas in the left hemisphere belonging to the somatomotor, dorsal attention, subcortical and visual networks. Overall, the results of our whole-brain model-based approach revealed a mechanism underlying meditation by providing causal relationships at the structure-function level.

Rare long-range cortical connections enhance human information processing.

What are the key topological features of connectivity critically relevant for generating the dynamics underlying efficient cortical function? A candidate feature that has recently emerged is that the connectivity of the mammalian cortex follows an exponential distance rule, which includes a small proportion of long-range high-weight anatomical exceptions to this rule. Whole-brain modeling of large-scale human neuroimaging data in 1,003 participants offers the unique opportunity to create two models, with and without long-range exceptions, and explicitly study their functional consequences. We found that rare long-range exceptions are crucial for significantly improving information processing. Furthermore, modeling in a simplified ring architecture shows that this improvement is greatly enhanced by the turbulent regime found in empirical neuroimaging data. Overall, the results provide strong empirical evidence for the immense functional benefits of long-range exceptions combined with turbulence for information processing.

Whole-brain modelling of resting state fMRI differentiates ADHD subtypes and facilitates stratified neuro-stimulation therapy.

Recent advances in non-linear computational and dynamical modelling have opened up the possibility to parametrize dynamic neural mechanisms that drive complex behavior. Importantly, building models of neuronal processes is of key importance to fully understand disorders of the brain as it may provide a quantitative platform that is capable of binding multiple neurophysiological processes to phenotype profiles. In this study, we apply a newly developed adaptive frequency-based model of whole-brain oscillations to resting-state fMRI data acquired from healthy controls and a cohort of attention deficit hyperactivity disorder (ADHD) subjects. As expected, we found that healthy control subjects differed from ADHD in terms of attractor dynamics. However, we also found a marked dichotomy in neural dynamics within the ADHD cohort. Next, we classified the ADHD group according to the level of distance of each individual's empirical network from the two model-based simulated networks. Critically, the model was mirrored in the empirical behavior data with the two ADHD subgroups displaying distinct behavioral phenotypes related to emotional instability (i.e., depression and hypomanic personality traits). Finally, we investigated the applicability and feasibility of our whole-brain model in a therapeutic setting by conducting in silico excitatory stimulations to parsimoniously mimic clinical neuro-stimulation paradigms in ADHD. We tested the effect of stimulating any individual brain region on the key network measures derived from the simulated brain network and its contribution in rectifying the brain dynamics to that of the healthy brain, separately for each ADHD subgroup. This showed that this was indeed possible for both subgroups. However, the current effect sizes were small suggesting that the stimulation protocol needs to be tailored at the individual level. These findings demonstrate the potential of this new modelling framework to unveil hidden neurophysiological profiles and establish tailored clinical interventions.

Dynamic coupling of whole-brain neuronal and neurotransmitter systems.

Remarkable progress has come from whole-brain models linking anatomy and function. Paradoxically, it is not clear how a neuronal dynamical system running in the fixed human anatomical connectome can give rise to the rich changes in the functional repertoire associated with human brain function, which is impossible to explain through long-term plasticity. Neuromodulation evolved to allow for such flexibility by dynamically updating the effectivity of the fixed anatomical connectivity. Here, we introduce a theoretical framework modeling the dynamical mutual coupling between the neuronal and neurotransmitter systems. We demonstrate that this framework is crucial to advance our understanding of whole-brain dynamics by bidirectional coupling of the two systems through combining multimodal neuroimaging data (diffusion magnetic resonance imaging [dMRI], functional magnetic resonance imaging [fMRI], and positron electron tomography [PET]) to explain the functional effects of specific serotoninergic receptor (5-HT2AR) stimulation with psilocybin in healthy humans. This advance provides an understanding of why psilocybin is showing considerable promise as a therapeutic intervention for neuropsychiatric disorders including depression, anxiety, and addiction. Overall, these insights demonstrate that the whole-brain mutual coupling between the neuronal and the neurotransmission systems is essential for understanding the remarkable flexibility of human brain function despite having to rely on fixed anatomical connectivity.

Imaging Connectomics and the Understanding of Brain Diseases.

Neuroimaging-based personalized medicine is emerging to characterize brain disorders and their evolution at the patient level. In this chapter, we present the most classic methods used to infer large-scale brain connectivity based on functional MRI. We adopt a modeling perspective where every connectivity measure is linked to a specific model that allows to interpret the connectivity estimate. This perspective allows to analyze the quality of retrieved connectivity profiles in terms of modeling error and estimation error. In the first part of the chapter, we present undirected functional connectivity (Pearson's correlation and MI) and effective connectivity (partial correlation), as well as directed effective connectivity (VAR, MOU, Granger causality, DCM). In addition, some of these measures correspond to fully connected graphs (Pearson's correlation) while others to sparse ones (MOU, DCM), where the sparsity can come from the integration of functional and structural data. In the second part, we claim that machine learning tools are better suited than null-hypothesis testing to link the estimated connectomes with diagnosis and prognosis of neuropsychiatric diseases. Finally, we propose that linear models and features selection are preferable to more complex and nonlinear tools (when prediction performance is on a par) for building interpretable algorithms to predict clinical variables.

Reliable local dynamics in the brain across sessions are revealed by whole-brain modeling of resting state activity.

Resting state fMRI is a tool for studying the functional organization of the human brain. Ongoing brain activity at "rest" is highly dynamic, but procedures such as correlation or independent component analysis treat functional connectivity (FC) as if, theoretically, it is stationary and therefore the fluctuations observed in FC are thought as noise. Consequently, FC is not usually used as a single-subject level marker and it is limited to group studies. Here we develop an imaging-based technique capable of reliably portraying information of local dynamics at a single-subject level by using a whole-brain model of ongoing dynamics that estimates a local parameter, which reflects if each brain region presents stable, asynchronous or transitory oscillations. Using 50 longitudinal resting-state sessions of one single subject and single resting-state sessions from a group of 50 participants we demonstrate that brain dynamics can be quantified consistently with respect to group dynamics using a scanning time of 20 min. We show that brain hubs are closer to a transition point between synchronous and asynchronous oscillatory dynamics and that dynamics in frontal areas have larger heterogeneity in its values compared to other lobules. Nevertheless, frontal regions and hubs showed higher consistency within the same subject while the inter-session variability found in primary visual and motor areas was only as high as the one found across subjects. The framework presented here can be used to study functional brain dynamics at group and, more importantly, at individual level, opening new avenues for possible clinical applications.

Poor Sleep Quality Associates With Decreased Functional and Structural Brain Connectivity in Normative Aging: A MRI Multimodal Approach.

Sleep is a ubiquitous phenomenon, essential to the organism homeostasis. Notwithstanding, there has been an increasing concern with its disruption, not only within the context of pathological conditions, such as neurologic and psychiatric diseases, but also in health. In fact, sleep complaints are becoming particularly common, especially in middle-aged and older adults, which may suggest an underlying susceptibility to sleep quality loss and/or its consequences. Thus, a whole-brain modeling approach to study the shifts in the system can cast broader light on sleep quality mechanisms and its associated morbidities. Following this line, we sought to determine the association between the standard self-reported measure of sleep quality, the Pittsburgh Sleep Quality Index (PSQI) and brain correlates in a normative aging cohort. To this purpose, 86 participants (age range 52-87 years) provided information regarding sociodemographic parameters, subjective sleep quality and associated psychological variables. A multimodal magnetic resonance imaging (MRI) approach was used, with whole-brain functional and structural connectomes being derived from resting-state functional connectivity (FC) and probabilistic white matter tractography (structural connectivity, SC). Brain regional volumes and white matter properties associations were also explored. Results show that poor sleep quality was associated with a decrease in FC and SC of distinct networks, overlapping in right superior temporal pole, left middle temporal and left inferior occipital regions. Age displayed important associations with volumetric changes in the cerebellum cortex and white matter, thalamus, hippocampus, right putamen, left supramarginal and left lingual regions. Overall, results suggest that not only the PSQI global score may act as a proxy of changes in FC/SC in middle-aged and older individuals, but also that the age-related regional volumetric changes may be associated to an adjustment of brain connectivity. These findings may also represent a step further in the comprehension of the role of sleep disturbance in disease, since the networks found share regions that have been shown to be affected in pathologies, such as depression and Alzheimer's disease.