SARS-CoV-2 is associated with changes in brain structure in UK Biobank.

There is strong evidence of brain-related abnormalities in COVID-191-13. However, it remains unknown whether the impact of SARS-CoV-2 infection can be detected in milder cases, and whether this can reveal possible mechanisms contributing to brain pathology. Here we investigated brain changes in 785 participants of UK Biobank (aged 51-81 years) who were imaged twice using magnetic resonance imaging, including 401 cases who tested positive for infection with SARS-CoV-2 between their two scans-with 141 days on average separating their diagnosis and the second scan-as well as 384 controls. The availability of pre-infection imaging data reduces the likelihood of pre-existing risk factors being misinterpreted as disease effects. We identified significant longitudinal effects when comparing the two groups, including (1) a greater reduction in grey matter thickness and tissue contrast in the orbitofrontal cortex and parahippocampal gyrus; (2) greater changes in markers of tissue damage in regions that are functionally connected to the primary olfactory cortex; and (3) a greater reduction in global brain size in the SARS-CoV-2 cases. The participants who were infected with SARS-CoV-2 also showed on average a greater cognitive decline between the two time points. Importantly, these imaging and cognitive longitudinal effects were still observed after excluding the 15 patients who had been hospitalised. These mainly limbic brain imaging results may be the in vivo hallmarks of a degenerative spread of the disease through olfactory pathways, of neuroinflammatory events, or of the loss of sensory input due to anosmia. Whether this deleterious effect can be partially reversed, or whether these effects will persist in the long term, remains to be investigated with additional follow-up.

Widespread, depth-dependent cortical microstructure alterations in pediatric focal epilepsy.

OBJECTIVE: Tissue abnormalities in focal epilepsy may extend beyond the presumed focus. The underlying pathophysiology of these broader changes is unclear, and it is not known whether they result from ongoing disease processes or treatment-related side effects, or whether they emerge earlier. Few studies have focused on the period of onset for most focal epilepsies, childhood. Fewer still have utilized quantitative magnetic resonance imaging (MRI), which may provide a more sensitive and interpretable measure of tissue microstructural change. Here, we aimed to determine common spatial modes of changes in cortical architecture in children with heterogeneous drug-resistant focal epilepsy and, secondarily, whether changes were related to disease severity. METHODS: To assess cortical microstructure, quantitative T1 and T2 relaxometry (qT1 and qT2) was measured in 43 children with drug-resistant focal epilepsy (age range = 4-18 years) and 46 typically developing children (age range = 2-18 years). We assessed depth-dependent qT1 and qT2 values across the neocortex, as well as their gradient of change across cortical depths. We also determined whether global changes seen in group analyses were driven by focal pathologies in individual patients. Finally, as a proof-of-concept, we trained a classifier using qT1 and qT2 gradient maps from patients with radiologically defined abnormalities (MRI positive) and healthy controls, and tested whether this could classify patients without reported radiological abnormalities (MRI negative). RESULTS: We uncovered depth-dependent qT1 and qT2 increases in widespread cortical areas in patients, likely representing microstructural alterations in myelin or gliosis. Changes did not correlate with disease severity measures, suggesting they may represent antecedent neurobiological alterations. Using a classifier trained with MRI-positive patients and controls, sensitivity was 71.4% at 89.4% specificity on held-out MRI-negative patients. SIGNIFICANCE: These findings suggest the presence of a potential imaging endophenotype of focal epilepsy, detectable irrespective of radiologically identified abnormalities.

Multimodal study of the neural sources of error monitoring in adolescents and adults.

The ability to monitor performance during a goal-directed behavior differs among children and adults in ways that can be measured with several tasks and techniques. As well, recent work has shown that individual differences in error monitoring moderate temperamental risk for anxiety and that this moderation changes with age. We investigated age differences in neural responses linked to performance monitoring using a multimodal approach. The approach combined functional MRI and source localization of event-related potentials (ERPs) in 12-year-old, 15-year-old, and adult participants. Neural generators of two components related to performance and error monitoring, the N2 and ERN, lay within specific areas of fMRI clusters. Whereas correlates of the N2 component appeared similar across age groups, age-related differences manifested in the location of the generators of the ERN component. The dorsal anterior cingulate cortex (dACC) was the predominant source location for the 12-year-old group; this area manifested posteriorly for the 15-year-old and adult groups. A fMRI-based ROI analysis confirmed this pattern of activity. These results suggest that changes in the underlying neural mechanisms are related to developmental changes in performance monitoring.

Multimodal Imaging Brain Markers in Early Adolescence Are Linked with a Physically Active Lifestyle.

The World Health Organization promotes physical exercise and a healthy lifestyle as means to improve youth development. However, relationships between physical lifestyle and human brain development are not fully understood. Here, we asked whether a human brain-physical latent mode of covariation underpins the relationship between physical activity, fitness, and physical health measures with multimodal neuroimaging markers. In 50 12-year old school pupils (26 females), we acquired multimodal whole-brain MRI, characterizing brain structure, microstructure, function, myelin content, and blood perfusion. We also acquired physical variables measuring objective fitness levels, 7 d physical activity, body mass index, heart rate, and blood pressure. Using canonical correlation analysis, we unravel a latent mode of brain-physical covariation, independent of demographics, school, or socioeconomic status. We show that MRI metrics with greater involvement in this mode also showed spatially extended patterns across the brain. Specifically, global patterns of greater gray matter perfusion, volume, cortical surface area, greater white matter extra-neurite density, and resting state networks activity covaried positively with measures reflecting a physically active phenotype (high fit, low sedentary individuals). Showing that a physically active lifestyle is linked with systems-level brain MRI metrics, these results suggest widespread associations relating to several biological processes. These results support the notion of close brain-body relationships and underline the importance of investigating modifiable lifestyle factors not only for physical health but also for brain health early in adolescence.SIGNIFICANCE STATEMENT An active lifestyle is key for healthy development. In this work, we answer the following question: How do brain neuroimaging markers relate with young adolescents' level of physical activity, fitness, and physical health? Combining advanced whole-brain multimodal MRI metrics with computational approaches, we show a robust relationship between physically active lifestyles and spatially extended, multimodal brain imaging-derived phenotypes. Suggesting a wider effect on brain neuroimaging metrics than previously thought, this work underlies the importance of studying physical lifestyle, as well as other brain-body relationships in an effort to foster brain health at this crucial stage in development.

The Enhancing NeuroImaging Genetics through Meta-Analysis Consortium: 10 Years of Global Collaborations in Human Brain Mapping.

This Special Issue of Human Brain Mapping is dedicated to a 10-year anniversary of the Enhancing NeuroImaging Genetics through Meta-Analysis (ENIGMA) Consortium. It reports updates from a broad range of international neuroimaging projects that pool data from around the world to answer fundamental questions in neuroscience. Since ENIGMA was formed in December 2009, the initiative grew into a worldwide effort with over 2,000 participating scientists from 45 countries, and over 50 working groups leading large-scale studies of human brain disorders. Over the last decade, many lessons were learned on how best to pool brain data from diverse sources. Working groups were created to develop methods to analyze worldwide data from anatomical and diffusion magnetic resonance imaging (MRI), resting state and task-based functional MRI, electroencephalography (EEG), magnetoencephalography (MEG), and magnetic resonance spectroscopy (MRS). The quest to understand genetic effects on human brain development and disease also led to analyses of brain scans on an unprecedented scale. Genetic roadmaps of the human cortex were created by researchers worldwide who collaborated to perform statistically well-powered analyses of common and rare genetic variants on brain measures and rates of brain development and aging. Here, we summarize the 31 papers in this Special Issue, covering: (a) technical approaches to harmonize analysis of different types of brain imaging data, (b) reviews of the last decade of work by several of ENIGMA's clinical and technical working groups, and (c) new empirical papers reporting large-scale international brain mapping analyses in patients with substance use disorders, schizophrenia, bipolar disorders, major depression, posttraumatic stress disorder, obsessive compulsive disorder, epilepsy, and stroke.

Patients with chronic pain exhibit individually unique cortical signatures of pain encoding.

Chronic pain is characterised by an ongoing and fluctuating intensity over time. Here, we investigated how the trajectory of the patients' endogenous pain is encoded in the brain. In repeated functional MRI (fMRI) sessions, 20 patients with chronic back pain and 20 patients with chronic migraine were asked to continuously rate the intensity of their endogenous pain. Linear mixed effects models were used to disentangle cortical processes related to pain intensity and to pain intensity changes. At group level, we found that the intensity of pain in patients with chronic back pain is encoded in the anterior insular cortex, the frontal operculum, and the pons; the change of pain in chronic back pain and chronic migraine patients is mainly encoded in the anterior insular cortex. At the individual level, we identified a more complex picture where each patient exhibited their own signature of endogenous pain encoding. The diversity of the individual cortical signatures of chronic pain encoding results bridge between clinical observations and neuroimaging; they add to the understanding of chronic pain as a complex and multifaceted disease.

ENIGMA-anxiety working group: Rationale for and organization of large-scale neuroimaging studies of anxiety disorders.

Anxiety disorders are highly prevalent and disabling but seem particularly tractable to investigation with translational neuroscience methodologies. Neuroimaging has informed our understanding of the neurobiology of anxiety disorders, but research has been limited by small sample sizes and low statistical power, as well as heterogenous imaging methodology. The ENIGMA-Anxiety Working Group has brought together researchers from around the world, in a harmonized and coordinated effort to address these challenges and generate more robust and reproducible findings. This paper elaborates on the concepts and methods informing the work of the working group to date, and describes the initial approach of the four subgroups studying generalized anxiety disorder, panic disorder, social anxiety disorder, and specific phobia. At present, the ENIGMA-Anxiety database contains information about more than 100 unique samples, from 16 countries and 59 institutes. Future directions include examining additional imaging modalities, integrating imaging and genetic data, and collaborating with other ENIGMA working groups. The ENIGMA consortium creates synergy at the intersection of global mental health and clinical neuroscience, and the ENIGMA-Anxiety Working Group extends the promise of this approach to neuroimaging research on anxiety disorders.

Translating Big Data to Clinical Outcomes in Anxiety: Potential for Multimodal Integration.

PURPOSE OF THE REVIEW: This review describes approaches to research on anxiety that attempt to link neural correlates to treatment response and novel therapies. The review emphasizes pediatric anxiety disorders since most anxiety disorders begin before adulthood. RECENT FINDINGS: Recent literature illustrates how current treatments for anxiety manifest diverse relations with a range of neural markers. While some studies demonstrate post-treatment normalization of markers in anxious individuals, others find persistence of group differences. For other markers, which show no pretreatment association with anxiety, the markers nevertheless distinguish treatment-responders from non-responders. Heightened error related negativity represents the risk marker discussed in the most depth; however, limitations in measures related to error responding necessitate multimodal and big-data approaches. Single risk markers show limits as correlates of treatment response. Large-scale, multimodal data analyzed with predictive models may illuminate additional risk markers related to anxiety disorder treatment outcomes. Such work may identify novel targets and eventually guide improvements in treatment response/outcomes.

An In-vivo 1H-MRS short-echo time technique at 7T: Quantification of metabolites in chronic multiple sclerosis and neuromyelitis optica brain lesions and normal appearing brain tissue.

Magnetic Resonance Spectroscopy (MRS) allows for the non-invasive quantification of neurochemicals and has the potential to differentiate between the pathologically distinct diseases, multiple sclerosis (MS) and AQP4Ab-positive neuromyelitis optica spectrum disorder (AQP4Ab-NMOSD). In this study we characterised the metabolite profiles of brain lesions in 11 MS and 4 AQP4Ab-NMOSD patients using an optimised MRS methodology at ultra-high field strength (7T) incorporating correction for T2 water relaxation differences between lesioned and normal tissue. MS metabolite results were in keeping with the existing literature: total N-acetylaspartate (NAA) was lower in lesions compared to normal appearing brain white matter (NAWM) with reciprocal findings for myo-Inositol. An unexpected subtlety revealed by our technique was that total NAA differences were likely driven by NAA-glutamate (NAAG), a ubiquitous CNS molecule with functions quite distinct from NAA though commonly quantified together with NAA in MRS studies as total NAA. Surprisingly, AQP4Ab-NMOSD showed no significant differences for total NAA, NAA, NAAG or myo-Inositol between lesion and NAWM sites, nor were there any differences between MS and AQP4Ab-NMOSD for a priori hypotheses. Post-hoc testing revealed a significant correlation between NAWM Ins:NAA and disability (as measured by EDSS) for disease groups combined, driven by the AP4Ab-NMOSD group. Utilising an optimised MRS methodology, our study highlights some under-explored subtleties in MRS profiles, such as the absence of myo-Inositol concentration differences in AQP4Ab-NMOSD brain lesions versus NAWM and the potential influence of NAAG differences between lesions and normal appearing white matter in MS.

Minimal Relationship between Local Gyrification and General Cognitive Ability in Humans.

Previous studies suggest that gyrification is associated with superior cognitive abilities in humans, but the strength of this relationship remains unclear. Here, in two samples of related individuals (total N = 2882), we calculated an index of local gyrification (LGI) at thousands of cortical surface points using structural brain images and an index of general cognitive ability (g) using performance on cognitive tests. Replicating previous studies, we found that phenotypic and genetic LGI-g correlations were positive and statistically significant in many cortical regions. However, all LGI-g correlations in both samples were extremely weak, regardless of whether they were significant or nonsignificant. For example, the median phenotypic LGI-g correlation was 0.05 in one sample and 0.10 in the other. These correlations were even weaker after adjusting for confounding neuroanatomical variables (intracranial volume and local cortical surface area). Furthermore, when all LGIs were considered together, at least 89% of the phenotypic variance of g remained unaccounted for. We conclude that the association between LGI and g is too weak to have profound implications for our understanding of the neurobiology of intelligence. This study highlights potential issues when focusing heavily on statistical significance rather than effect sizes in large-scale observational neuroimaging studies.

Neuroimaging Association Scores: reliability and validity of aggregate measures of brain structural features linked to mental disorders in youth.

In genetics, aggregation of many loci with small effect sizes into a single score improved prediction. Nevertheless, studies applying easily replicable weighted scores to neuroimaging data are lacking. Our aim was to assess the reliability and validity of the Neuroimaging Association Score (NAS), which combines information from structural brain features previously linked to mental disorders. Participants were 726 youth (aged 6-14) from two cities in Brazil who underwent MRI and psychopathology assessment at baseline and 387 at 3-year follow-up. Results were replicated in two samples: IMAGEN (n = 1627) and the Healthy Brain Network (n = 843). NAS were derived by summing the product of each standardized brain feature by the effect size of the association of that brain feature with seven psychiatric disorders documented by previous meta-analyses. NAS were calculated for surface area, cortical thickness and subcortical volumes using T1-weighted scans. NAS reliability, temporal stability and psychopathology and cognition prediction were analyzed. NAS for surface area showed high internal consistency and 3-year stability and predicted general psychopathology and cognition with higher replicability than specific symptomatic domains for all samples. They also predicted general psychopathology with higher replicability than single structures alone, accounting for 1-3% of the variance, but without directionality. The NAS for cortical thickness and subcortical volumes showed lower internal consistency and less replicable associations with behavioural phenotypes. These findings indicate the NAS based on surface area might be replicable markers of general psychopathology, but these links are unlikely to be causal or clinically useful yet.

Permutation inference for canonical correlation analysis.

Canonical correlation analysis (CCA) has become a key tool for population neuroimaging, allowing investigation of associations between many imaging and non-imaging measurements. As age, sex and other variables are often a source of variability not of direct interest, previous work has used CCA on residuals from a model that removes these effects, then proceeded directly to permutation inference. We show that a simple permutation test, as typically used to identify significant modes of shared variation on such data adjusted for nuisance variables, produces inflated error rates. The reason is that residualisation introduces dependencies among the observations that violate the exchangeability assumption. Even in the absence of nuisance variables, however, a simple permutation test for CCA also leads to excess error rates for all canonical correlations other than the first. The reason is that a simple permutation scheme does not ignore the variability already explained by previous canonical variables. Here we propose solutions for both problems: in the case of nuisance variables, we show that transforming the residuals to a lower dimensional basis where exchangeability holds results in a valid permutation test; for more general cases, with or without nuisance variables, we propose estimating the canonical correlations in a stepwise manner, removing at each iteration the variance already explained, while dealing with different number of variables in both sides. We also discuss how to address the multiplicity of tests, proposing an admissible test that is not conservative, and provide a complete algorithm for permutation inference for CCA.

Multiple testing correction over contrasts for brain imaging.

The multiple testing problem arises not only when there are many voxels or vertices in an image representation of the brain, but also when multiple contrasts of parameter estimates (that represent hypotheses) are tested in the same general linear model. We argue that a correction for this multiplicity must be performed to avoid excess of false positives. Various methods for correction have been proposed in the literature, but few have been applied to brain imaging. Here we discuss and compare different methods to make such correction in different scenarios, showing that one classical and well known method is invalid, and argue that permutation is the best option to perform such correction due to its exactness and flexibility to handle a variety of common imaging situations.

Combining fMRI during resting state and an attention bias task in children.

Neuroimaging studies typically focus on either resting state or task-based fMRI data. Prior research has shown that similarity in functional connectivity between rest and cognitive tasks, interpreted as reconfiguration efficiency, is related to task performance and IQ. Here, we extend this approach from adults to children, and from cognitive tasks to a threat-based attention task. The goal of the current study was to examine whether similarity in functional connectivity during rest and an attention bias task relates to threat bias, IQ, anxiety symptoms, and social reticence. fMRI was measured during resting state and during the dot-probe task in 41 children (M = 13.44, SD = 0.70). Functional connectivity during rest and dot-probe was positively correlated, suggesting that functional hierarchies in the brain are stable. Similarity in functional connectivity between rest and the dot-probe task only related to threat bias (puncorr < .03). This effect did not survive correction for multiple testing. Overall, children who allocate more attention towards threat also may possess greater reconfiguration efficiency in switching from intrinsic to threat-related attention states. Finally, functional connectivity correlated negatively across the two conditions of the dot-probe task. Opposing patterns of modulation of functional connectivity by threat-congruent and threat-incongruent trials may reflect task-specific network changes during two different attentional processes.

Schizophrenia Exhibits Bi-directional Brain-Wide Alterations in Cortico-Striato-Cerebellar Circuits.

Distributed neural dysconnectivity is considered a hallmark feature of schizophrenia (SCZ), yet a tension exists between studies pinpointing focal disruptions versus those implicating brain-wide disturbances. The cerebellum and the striatum communicate reciprocally with the thalamus and cortex through monosynaptic and polysynaptic connections, forming cortico-striatal-thalamic-cerebellar (CSTC) functional pathways that may be sensitive to brain-wide dysconnectivity in SCZ. It remains unknown if the same pattern of alterations persists across CSTC systems, or if specific alterations exist along key functional elements of these networks. We characterized connectivity along major functional CSTC subdivisions using resting-state functional magnetic resonance imaging in 159 chronic patients and 162 matched controls. Associative CSTC subdivisions revealed consistent brain-wide bi-directional alterations in patients, marked by hyper-connectivity with sensory-motor cortices and hypo-connectivity with association cortex. Focusing on the cerebellar and striatal components, we validate the effects using data-driven k-means clustering of voxel-wise dysconnectivity and support vector machine classifiers. We replicate these results in an independent sample of 202 controls and 145 patients, additionally demonstrating that these neural effects relate to cognitive performance across subjects. Taken together, these results from complementary approaches implicate a consistent motif of brain-wide alterations in CSTC systems in SCZ, calling into question accounts of exclusively focal functional disturbances.

Stable between-subject statistical inference from unstable within-subject functional connectivity estimates.

Spatial or temporal aspects of neural organization are known to be important indices of how cognition is organized. However, measurements and estimations are often noisy and many of the algorithms used are probabilistic, which in combination have been argued to limit studies exploring the neural basis of specific aspects of cognition. Focusing on static and dynamic functional connectivity estimations, we propose to leverage this variability to improve statistical efficiency in relating these estimations to behavior. To achieve this goal, we use a procedure based on permutation testing that provides a way of combining the results from many individual tests that refer to the same hypothesis. This is needed when testing a measure whose value is obtained from a noisy process, which can be repeated multiple times, referred to as replications. Focusing on functional connectivity, this noisy process can be: (a) computational, for example, when using an approximate inference algorithm for which different runs can produce different results or (b) observational, if we have the capacity to acquire data multiple times, and the different acquired data sets can be considered noisy examples of some underlying truth. In both cases, we are not interested in the individual replications but on the unobserved process generating each replication. In this note, we show how results can be combined instead of choosing just one of the estimated models. Using both simulations and real data, we show the benefits of this approach in practice.

Accelerated estimation and permutation inference for ACE modeling.

There are a wealth of tools for fitting linear models at each location in the brain in neuroimaging analysis, and a wealth of genetic tools for estimating heritability for a small number of phenotypes. But there remains a need for computationally efficient neuroimaging genetic tools that can conduct analyses at the brain-wide scale. Here we present a simple method for heritability estimation on twins that replaces a variance component model-which requires iterative optimisation-with a (noniterative) linear regression model, by transforming data to squared twin-pair differences. We demonstrate that the method has comparable bias, mean squared error, false positive risk, and power to best practice maximum-likelihood-based methods, while requiring a small fraction of the computation time. Combined with permutation, we call this approach "Accelerated Permutation Inference for the ACE Model (APACE)" where ACE refers to the additive genetic (A) effects, and common (C), and unique (E) environmental influences on the trait. We show how the use of spatial statistics like cluster size can dramatically improve power, and illustrate the method on a heritability analysis of an fMRI working memory dataset.

Joint Analysis of Cortical Area and Thickness as a Replacement for the Analysis of the Volume of the Cerebral Cortex.

Cortical surface area is an increasingly used brain morphology metric that is ontogenetically and phylogenetically distinct from cortical thickness and offers a separate index of neurodevelopment and disease. However, the various existing methods for assessment of cortical surface area from magnetic resonance images have never been systematically compared. We show that the surface area method implemented in FreeSurfer corresponds closely to the exact, but computationally more demanding, mass-conservative (pycnophylactic) method, provided that images are smoothed. Thus, the data produced by this method can be interpreted as estimates of cortical surface area, as opposed to areal expansion. In addition, focusing on the joint analysis of thickness and area, we compare an improved, analytic method for measuring cortical volume to a permutation-based nonparametric combination (NPC) method. We use the methods to analyze area, thickness and volume in young adults born preterm with very low birth weight, and show that NPC analysis is a more sensitive option for studying joint effects on area and thickness, giving equal weight to variation in both of these 2 morphological features.

Investigating resting-state functional connectivity in the cervical spinal cord at 3T.

The study of spontaneous fluctuations in the blood-oxygen-level-dependent (BOLD) signal has recently been extended from the brain to the spinal cord. Two ultra-high field functional magnetic resonance imaging (fMRI) studies in humans have provided evidence for reproducible resting-state connectivity between the dorsal horns as well as between the ventral horns, and a study in non-human primates has shown that these resting-state signals are impacted by spinal cord injury. As these studies were carried out at ultra-high field strengths using region-of-interest (ROI) based analyses, we investigated whether such resting-state signals could also be observed at the clinically more prevalent field strength of 3T. In a reanalysis of a sample of 20 healthy human participants who underwent a resting-state fMRI acquisition of the cervical spinal cord, we were able to observe significant dorsal horn connectivity as well as ventral horn connectivity, but no consistent effects for connectivity between dorsal and ventral horns, thus replicating the human 7T results. These effects were not only observable when averaging along the acquired length of the spinal cord, but also when we examined each of the acquired spinal segments separately, which showed similar patterns of connectivity. Finally, we investigated the robustness of these resting-state signals against variations in the analysis pipeline by varying the type of ROI creation, temporal filtering, nuisance regression and connectivity metric. We observed that - apart from the effects of band-pass filtering - ventral horn connectivity showed excellent robustness, whereas dorsal horn connectivity showed moderate robustness. Together, our results provide evidence that spinal cord resting-state connectivity is a robust and spatially consistent phenomenon that could be a valuable tool for investigating the effects of pathology, disease progression, and treatment response in neurological conditions with a spinal component, such as spinal cord injury.

Time related effects on functional brain connectivity after serotonergic and cholinergic neuromodulation.

Psychopharmacological research, if properly designed, may offer insight into both timing and area of effect, increasing our understanding of the brain's neurotransmitter systems. For that purpose, the acute influence of the selective serotonin reuptake inhibitor citalopram (30 mg) and the acetylcholinesterase inhibitor galantamine (8 mg) was repeatedly measured in 12 healthy young volunteers with resting state functional magnetic resonance imaging (RS-fMRI). Eighteen RS-fMRI scans were acquired per subject during this randomized, double blind, placebo-controlled, crossover study. Within-group comparisons of voxelwise functional connectivity with 10 functional networks were examined (P < 0.05, FWE-corrected) using a non-parametric multivariate approach with cerebrospinal fluid, white matter, heart rate, and baseline measurements as covariates. Although both compounds did not change cognitive performance on several tests, significant effects were found on connectivity with multiple resting state networks. Serotonergic stimulation primarily reduced connectivity with the sensorimotor network and structures that are related to self-referential mechanisms, whereas galantamine affected networks and regions that are more involved in learning, memory, and visual perception and processing. These results are consistent with the serotonergic and cholinergic trajectories and their functional relevance. In addition, this study demonstrates the power of using repeated measures after drug administration, which offers the chance to explore both combined and time specific effects. Hum Brain Mapp 38:308-325, 2017. © 2016 Wiley Periodicals, Inc.