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.

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.

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.

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.

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.

Inferring pathobiology from structural MRI in schizophrenia and bipolar disorder: Modeling head motion and neuroanatomical specificity.

Despite over 400 peer-reviewed structural MRI publications documenting neuroanatomic abnormalities in bipolar disorder and schizophrenia, the confounding effects of head motion and the regional specificity of these defects are unclear. Using a large cohort of individuals scanned on the same research dedicated MRI with broadly similar protocols, we observe reduced cortical thickness indices in both illnesses, though less pronounced in bipolar disorder. While schizophrenia (n = 226) was associated with wide-spread surface area reductions, bipolar disorder (n = 227) and healthy comparison subjects (n = 370) did not differ. We replicate earlier reports that head motion (estimated from time-series data) influences surface area and cortical thickness measurements and demonstrate that motion influences a portion, but not all, of the observed between-group structural differences. Although the effect sizes for these differences were small to medium, when global indices were covaried during vertex-level analyses, between-group effects became nonsignificant. This analysis raises doubts about the regional specificity of structural brain changes, possible in contrast to functional changes, in affective and psychotic illnesses as measured with current imaging technology. Given that both schizophrenia and bipolar disorder showed cortical thickness reductions, but only schizophrenia showed surface area changes, and assuming these measures are influenced by at least partially unique sets of biological factors, then our results could indicate some degree of specificity between bipolar disorder and schizophrenia. Hum Brain Mapp 38:3757-3770, 2017. © 2017 Wiley Periodicals, Inc.

Optimal echo time for functional MRI of the infant brain identified in response to noxious stimulation.

PURPOSE: Blood oxygen level dependent (BOLD) brain activity, measured using functional MRI (fMRI), is dependent on the echo time (TE) and the reversible spin-spin relaxation time constant ( T2*) that describes the decay of transverse magnetization. Use of the optimal TE during fMRI experiments allows maximal sensitivity to BOLD to be achieved. Reports that T2* values are longer in infants (due to higher water concentrations and lower lipid content) have led to the use of longer TEs during infant fMRI experiments; however, the optimal TE has not been established. METHODS: In this study, acute experimental mildly noxious stimuli were applied to the heel in 12 term infants (mean gestational age = 40 weeks, mean postnatal age = 3 days); and the percentage change in BOLD activity was calculated across a range of TEs, from 30 to 70 ms, at 3 Tesla. In addition, T2* maps of the whole brain were collected in seven infants. RESULTS: The maximal change in BOLD occurred at a TE of 52 ms, and the average T2* across the whole brain was 99 ms. CONCLUSION: A TE of approximately 50 ms is recommended for use in 3T fMRI investigations in term infants. Magn Reson Med 78:625-631, 2017. © 2016 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine.

Hyperinsulinemia and elevated systolic blood pressure independently predict white matter hyperintensities with associated cognitive decrement in the middle-aged offspring of dementia patients.

Cerebrovascular disease is an independent risk factor for dementia that may also be synergistic with Alzheimer's disease. In recent years attention has switched from cerebral infarcts to microvascular disease as the primary cause of cerebrovascular cognitive decline, with damage to the white matter the primary mechanism. Uncertainties remain regarding the risks posed by different types vascular threat, the extent to which cerebrovascular damage occurs in middle age, and whether relatively "normal" amounts of white matter damage are accompanied by meaningful degrees of cognitive decline. We explored these issues via laboratory, cardiovascular, cognitive, and magnetic resonance imaging (MRI) data in 67 middle-aged cognitively normal offspring of dementia patients. The sample was enriched for vascular risk. Plasma insulin, 24-h systolic blood pressure, body mass index, age, and % small dense LDL cholesterol were the strongest correlates of MRI white matter hyperintensity (WMH) volume. With shared variance controlled for, 24 h systolic BP, plasma insulin, and age remained as significant predictors of WMH volume. An interaction variable (24 h BP * insulin) did not improve the prediction of WMH. WMH volume correlated negatively with cognition. No evidence for an ApoE ε4 effect emerged for either WMH or cognition. Hypertension and hyperinsulinemia appear to pose independent, consequential threats to the cerebral small vessel vasculature in middle age, reflected in the presence of areas of WMH on MRI scans. Our data show that even modest WMH volumes in middle age are associated with cognitive decrement, underscoring the importance of aggressive treatment and lifestyle modifications to address vascular risk throughout adulthood.

Faster permutation inference in brain imaging.

Permutation tests are increasingly being used as a reliable method for inference in neuroimaging analysis. However, they are computationally intensive. For small, non-imaging datasets, recomputing a model thousands of times is seldom a problem, but for large, complex models this can be prohibitively slow, even with the availability of inexpensive computing power. Here we exploit properties of statistics used with the general linear model (GLM) and their distributions to obtain accelerations irrespective of generic software or hardware improvements. We compare the following approaches: (i) performing a small number of permutations; (ii) estimating the p-value as a parameter of a negative binomial distribution; (iii) fitting a generalised Pareto distribution to the tail of the permutation distribution; (iv) computing p-values based on the expected moments of the permutation distribution, approximated from a gamma distribution; (v) direct fitting of a gamma distribution to the empirical permutation distribution; and (vi) permuting a reduced number of voxels, with completion of the remainder using low rank matrix theory. Using synthetic data we assessed the different methods in terms of their error rates, power, agreement with a reference result, and the risk of taking a different decision regarding the rejection of the null hypotheses (known as the resampling risk). We also conducted a re-analysis of a voxel-based morphometry study as a real-data example. All methods yielded exact error rates. Likewise, power was similar across methods. Resampling risk was higher for methods (i), (iii) and (v). For comparable resampling risks, the method in which no permutations are done (iv) was the absolute fastest. All methods produced visually similar maps for the real data, with stronger effects being detected in the family-wise error rate corrected maps by (iii) and (v), and generally similar to the results seen in the reference set. Overall, for uncorrected p-values, method (iv) was found the best as long as symmetric errors can be assumed. In all other settings, including for familywise error corrected p-values, we recommend the tail approximation (iii). The methods considered are freely available in the tool PALM - Permutation Analysis of Linear Models.