Higher-order functional connectivity analysis of resting-state functional magnetic resonance imaging data using multivariate cumulants.

Blood-level oxygenation-dependent (BOLD) functional magnetic resonance imaging (fMRI) is the most common modality to study functional connectivity in the human brain. Most research to date has focused on connectivity between pairs of brain regions. However, attention has recently turned towards connectivity involving more than two regions, that is, higher-order connectivity. It is not yet clear how higher-order connectivity can best be quantified. The measures that are currently in use cannot distinguish between pairwise (i.e., second-order) and higher-order connectivity. We show that genuine higher-order connectivity can be quantified by using multivariate cumulants. We explore the use of multivariate cumulants for quantifying higher-order connectivity and the performance of block bootstrapping for statistical inference. In particular, we formulate a generative model for fMRI signals exhibiting higher-order connectivity and use it to assess bias, standard errors, and detection probabilities. Application to resting-state fMRI data from the Human Connectome Project demonstrates that spontaneous fMRI signals are organized into higher-order networks that are distinct from second-order resting-state networks. Application to a clinical cohort of patients with multiple sclerosis further demonstrates that cumulants can be used to classify disease groups and explain behavioral variability. Hence, we present a novel framework to reliably estimate genuine higher-order connectivity in fMRI data which can be used for constructing hyperedges, and finally, which can readily be applied to fMRI data from populations with neuropsychiatric disease or cognitive neuroscientific experiments.

Non-reversibility outperforms functional connectivity in characterisation of brain states in MEG data.

Characterising brain states during tasks is common practice for many neuroscientific experiments using electrophysiological modalities such as electroencephalography (EEG) and magnetoencephalography (MEG). Brain states are often described in terms of oscillatory power and correlated brain activity, i.e. functional connectivity. It is, however, not unusual to observe weak task induced functional connectivity alterations in the presence of strong task induced power modulations using classical time-frequency representation of the data. Here, we propose that non-reversibility, or the temporal asymmetry in functional interactions, may be more sensitive to characterise task induced brain states than functional connectivity. As a second step, we explore causal mechanisms of non-reversibility in MEG data using whole brain computational models. We include working memory, motor, language tasks and resting-state data from participants of the Human Connectome Project (HCP). Non-reversibility is derived from the lagged amplitude envelope correlation (LAEC), and is based on asymmetry of the forward and reversed cross-correlations of the amplitude envelopes. Using random forests, we find that non-reversibility outperforms functional connectivity in the identification of task induced brain states. Non-reversibility shows especially better sensitivity to capture bottom-up gamma induced brain states across all tasks, but also alpha band associated brain states. Using whole brain computational models we find that asymmetry in the effective connectivity and axonal conduction delays play a major role in shaping non-reversibility across the brain. Our work paves the way for better sensitivity in characterising brain states during both bottom-up as well as top-down modulation in future neuroscientific experiments.

Relationships between cortical myeloarchitecture and electrophysiological networks.

The human brain relies upon the dynamic formation and dissolution of a hierarchy of functional networks to support ongoing cognition. However, how functional connectivities underlying such networks are supported by cortical microstructure remains poorly understood. Recent animal work has demonstrated that electrical activity promotes myelination. Inspired by this, we test a hypothesis that gray-matter myelin is related to electrophysiological connectivity. Using ultra-high field MRI and the principle of structural covariance, we derive a structural network showing how myelin density differs across cortical regions and how separate regions can exhibit similar myeloarchitecture. Building upon recent evidence that neural oscillations mediate connectivity, we use magnetoencephalography to elucidate networks that represent the major electrophysiological pathways of communication in the brain. Finally, we show that a significant relationship exists between our functional and structural networks; this relationship differs as a function of neural oscillatory frequency and becomes stronger when integrating oscillations over frequency bands. Our study sheds light on the way in which cortical microstructure supports functional networks. Further, it paves the way for future investigations of the gray-matter structure/function relationship and its breakdown in pathology.

Optimising experimental design for MEG resting state functional connectivity measurement.

The study of functional connectivity using magnetoencephalography (MEG) is an expanding area of neuroimaging, and adds an extra dimension to the more common assessments made using fMRI. The importance of such metrics is growing, with recent demonstrations of their utility in clinical research, however previous reports suggest that whilst group level resting state connectivity is robust, single session recordings lack repeatability. Such robustness is critical if MEG measures in individual subjects are to prove clinically valuable. In the present paper, we test how practical aspects of experimental design affect the intra-subject repeatability of MEG findings; specifically we assess the effect of co-registration method and data recording duration. We show that the use of a foam head-cast, which is known to improve co-registration accuracy, increased significantly the between session repeatability of both beamformer reconstruction and connectivity estimation. We also show that recording duration is a critical parameter, with large improvements in repeatability apparent when using ten minute, compared to five minute recordings. Further analyses suggest that the origin of this latter effect is not underpinned by technical aspects of source reconstruction, but rather by a genuine effect of brain state; short recordings are simply inefficient at capturing the canonical MEG network in a single subject. Our results provide important insights on experimental design and will prove valuable for future MEG connectivity studies.

Measurement of dynamic task related functional networks using MEG.

The characterisation of dynamic electrophysiological brain networks, which form and dissolve in order to support ongoing cognitive function, is one of the most important goals in neuroscience. Here, we introduce a method for measuring such networks in the human brain using magnetoencephalography (MEG). Previous network analyses look for brain regions that share a common temporal profile of activity. Here distinctly, we exploit the high spatio-temporal resolution of MEG to measure the temporal evolution of connectivity between pairs of parcellated brain regions. We then use an ICA based procedure to identify networks of connections whose temporal dynamics covary. We validate our method using MEG data recorded during a finger movement task, identifying a transient network of connections linking somatosensory and primary motor regions, which modulates during the task. Next, we use our method to image the networks which support cognition during a Sternberg working memory task. We generate a novel neuroscientific picture of cognitive processing, showing the formation and dissolution of multiple networks which relate to semantic processing, pattern recognition and language as well as vision and movement. Our method tracks the dynamics of functional connectivity in the brain on a timescale commensurate to the task they are undertaking.

Abnormal task driven neural oscillations in multiple sclerosis: A visuomotor MEG study.

Multiple sclerosis (MS) is a debilitating disease commonly attributed to degradation of white matter myelin. Symptoms include fatigue, as well as problems associated with vision and movement. Although areas of demyelination in white matter are observed routinely in patients undergoing MRI scans, such measures are often a poor predictor of disease severity. For this reason, it is instructive to measure associated changes in brain function. Widespread white-matter demyelination may lead to delays of propagation of neuronal activity, and with its excellent temporal resolution, magnetoencephalography can be used to probe such delays in controlled conditions (e.g., during a task). In healthy subjects, responses to visuomotor tasks are well documented: in motor cortex, movement elicits a localised decrease in the power of beta band oscillations (event-related beta desynchronisation) followed by an increase above baseline on movement cessation (post-movement beta rebound (PMBR)). In visual cortex, visual stimulation generates increased gamma oscillations. In this study, we use a visuomotor paradigm to measure these responses in MS patients and compare them to age- and gender-matched healthy controls. We show a significant increase in the time-to-peak of the PMBR in patients which correlates significantly with the symbol digit modalities test: a measure of information processing speed. A significant decrease in the amplitude of visual gamma oscillations in patients is also seen. These findings highlight the potential value of electrophysiological imaging in generating a new understanding of visual disturbances and abnormal motor control in MS patients. Hum Brain Mapp 38:2441-2453, 2017. © 2017 Wiley Periodicals, Inc.

Cortical atrophy patterns in multiple sclerosis are non-random and clinically relevant.

Grey matter atrophy is common in multiple sclerosis. However, in contrast with other neurodegenerative diseases, it is unclear whether grey matter atrophy in multiple sclerosis is a diffuse 'global' process or develops, instead, according to distinct anatomical patterns. Using source-based morphometry we searched for anatomical patterns of co-varying cortical thickness and assessed their relationships with white matter pathology, physical disability and cognitive functioning. Magnetic resonance imaging was performed at 3 T in 208 patients with long-standing multiple sclerosis (141 females; age = 53.7 ± 9.6 years; disease duration = 20.2 ± 7.1 years) and 60 age- and sex-matched healthy controls. Spatial independent component analysis was performed on cortical thickness maps derived from 3D T1-weighted images across all subjects to identify co-varying patterns. The loadings, which reflect the presence of each cortical thickness pattern in a subject, were compared between patients with multiple sclerosis and healthy controls with generalized linear models. Stepwise linear regression analyses were used to assess whether white matter pathology was associated with these loadings and to identify the cortical thickness patterns that predict measures of physical and cognitive dysfunction. Ten cortical thickness patterns were identified, of which six had significantly lower loadings in patients with multiple sclerosis than in controls: the largest loading differences corresponded to the pattern predominantly involving the bilateral temporal pole and entorhinal cortex, and the pattern involving the bilateral posterior cingulate cortex. In patients with multiple sclerosis, overall white matter lesion load was negatively associated with the loadings of these two patterns. The final model for physical dysfunction as measured with Expanded Disability Status Scale score (adjusted R(2) = 0.297; P < 0.001) included the predictors age, overall white matter lesion load, the loadings of two cortical thickness patterns (bilateral sensorimotor cortex and bilateral insula), and global cortical thickness. The final model predicting average cognition (adjusted R(2) = 0.469; P < 0.001) consisted of age, the loadings of two cortical thickness patterns (bilateral posterior cingulate cortex and bilateral temporal pole), overall white matter lesion load and normal-appearing white matter integrity. Although white matter pathology measures were part of the final clinical regression models, they explained limited incremental variance (to a maximum of 4%). Several cortical atrophy patterns relevant for multiple sclerosis were found. This suggests that cortical atrophy in multiple sclerosis occurs largely in a non-random manner and develops (at least partly) according to distinct anatomical patterns. In addition, these cortical atrophy patterns showed stronger associations with clinical (especially cognitive) dysfunction than global cortical atrophy.

A multi-layer network approach to MEG connectivity analysis.

Recent years have shown the critical importance of inter-regional neural network connectivity in supporting healthy brain function. Such connectivity is measurable using neuroimaging techniques such as MEG, however the richness of the electrophysiological signal makes gaining a complete picture challenging. Specifically, connectivity can be calculated as statistical interdependencies between neural oscillations within a large range of different frequency bands. Further, connectivity can be computed between frequency bands. This pan-spectral network hierarchy likely helps to mediate simultaneous formation of multiple brain networks, which support ongoing task demand. However, to date it has been largely overlooked, with many electrophysiological functional connectivity studies treating individual frequency bands in isolation. Here, we combine oscillatory envelope based functional connectivity metrics with a multi-layer network framework in order to derive a more complete picture of connectivity within and between frequencies. We test this methodology using MEG data recorded during a visuomotor task, highlighting simultaneous and transient formation of motor networks in the beta band, visual networks in the gamma band and a beta to gamma interaction. Having tested our method, we use it to demonstrate differences in occipital alpha band connectivity in patients with schizophrenia compared to healthy controls. We further show that these connectivity differences are predictive of the severity of persistent symptoms of the disease, highlighting their clinical relevance. Our findings demonstrate the unique potential of MEG to characterise neural network formation and dissolution. Further, we add weight to the argument that dysconnectivity is a core feature of the neuropathology underlying schizophrenia.

Multi-parametric structural magnetic resonance imaging in relation to cognitive dysfunction in long-standing multiple sclerosis.

BACKGROUND: Cognitive deficits are common in multiple sclerosis. Most previous studies investigating the imaging substrate of cognitive deficits in multiple sclerosis included patients with relatively short disease durations and were limited to one modality/brain region. OBJECTIVE: To identify the strongest neuroimaging predictors for cognitive dysfunction in a large cohort of patients with long-standing multiple sclerosis. METHODS: Extensive neuropsychological testing and multimodal 3.0T MRI was performed in 202 patients with multiple sclerosis and 52 controls. Cognitive scores were compared between groups using Z-scores. Whole-brain, white matter, grey matter, deep grey matter and lesion volumes; cortical thickness, (juxta)cortical and cerebellar lesions; and extent and severity of diffuse white matter damage were measured. Stepwise linear regression was used to identify the strongest predictors for cognitive dysfunction. RESULTS: All cognitive domains were affected in patients. Patients showed extensive atrophy, focal pathology and damage in up to 75% of the investigated white matter. Associations between imaging markers and average cognition were two times stronger in cognitively impaired patients than in cognitively preserved patients. The final model for average cognition consisted of deep grey matter DGMV volume and fractional anisotropy severity (adjusted R²=0.490; p<0.001). CONCLUSION: From all imaging markers, deep grey matter atrophy and diffuse white matter damage emerged as the strongest predictors for cognitive dysfunction in long-standing multiple sclerosis.

Unraveling the relationship between regional gray matter atrophy and pathology in connected white matter tracts in long-standing multiple sclerosis.

INTRODUCTION: Gray matter (GM) atrophy is common in multiple sclerosis (MS), but the relationship with white matter (WM) pathology is largely unknown. Some studies found a co-occurrence in specific systems, but a regional analysis across the brain in different clinical phenotypes is necessary to further understand the disease mechanism underlying GM atrophy in MS. Therefore, we investigated the association between regional GM atrophy and pathology in anatomically connected WM tracts. METHODS: Conventional and diffusion tensor imaging was performed at 3T in 208 patients with long-standing MS and 60 healthy controls. Deep and cortical GM regions were segmented and quantified, and both lesion volumes and average normal appearing WM fractional anisotropy of their associated tracts were derived using an atlas obtained by probabilistic tractography in the controls. Linear regression was then performed to quantify the amount of regional GM atrophy that can be explained by WM pathology in the connected tract. RESULTS: MS patients showed extensive deep and cortical GM atrophy. Cortical atrophy was particularly present in frontal and temporal regions. Pathology in connected WM tracts statistically explained both regional deep and cortical GM atrophy in relapsing-remitting (RR) patients, but only deep GM atrophy in secondary-progressive (SP) patients. CONCLUSION: In RRMS patients, both deep and cortical GM atrophy were associated with pathology in connected WM tracts. In SPMS patients, only regional deep GM atrophy could be explained by pathology in connected WM tracts. This suggests that in SPMS patients cortical GM atrophy and WM damage are (at least partly) independent disease processes.

What explains gray matter atrophy in long-standing multiple sclerosis?

PURPOSE: To identify the measures of focal and diffuse white matter (WM) abnormalities that are related to whole-brain, deep, and cortical gray matter (GM) atrophy in long-standing multiple sclerosis (MS). MATERIALS AND METHODS: The institutional review board approved the study; all subjects gave written informed consent. Magnetic resonance (MR) imaging was performed at 3 T in 208 patients with MS of long-standing duration (disease duration ≥ 10 years) and in 60 healthy control subjects. Normalized GM volume (NGMV), normalized WM volume (NWMV), normalized deep GM volume (NDGMV), cortical thickness, and normalized lesion volume (NLV) were quantified. Tissue integrity of normal-appearing WM (NAWM) and lesions was measured by using diffusion-tensor MR imaging. Multivariate associations between measures of GM atrophy and WM abnormalities were assessed in the patient group by using multiple linear regression. RESULTS: NGMV, NDGMV, and cortical thickness were reduced in patients with MS (all P < .001). The final model for NGMV consisted of NWMV, NLV, and patient age and sex (adjusted R(2) = 0.58, P < .001). NWMV, NLV, and patient sex were the explanatory variables for NDGMV (adjusted R(2) = 0.75, P < .001). The model for cortical thickness consisted of fractional anisotropy of NAWM, NLV, and patient age and sex (adjusted R(2) = 0.32, P < .001). The relationship between GM atrophy and WM abnormalities was weaker in primary and secondary progressive disease than in relapsing-remitting disease. CONCLUSION: Whole-brain and deep GM atrophy were particularly explained by WM atrophy and lesion volume, while cortical atrophy was associated with NAWM integrity loss. The weaker relationship between GM atrophy and WM abnormalities in patients with progressive disease might indicate a more independent neurodegenerative disease process in these patients.

Mean upper cervical cord area (MUCCA) measurement in long-standing multiple sclerosis: relation to brain findings and clinical disability.

BACKGROUND: The majority of patients with multiple sclerosis (MS) present with spinal cord pathology. Spinal cord atrophy is thought to be a marker of disease severity, but in long-disease duration its relation to brain pathology and clinical disability is largely unknown. OBJECTIVE: Our aim was to investigate mean upper cervical cord area (MUCCA) in patients with long-standing MS and assess its relation to brain magnetic resonance imaging (MRI) measures and clinical disability. METHODS: MUCCA was measured in 196 MS patients and 55 healthy controls using 3DT1-weighted cervical images obtained at 3T MRI. Clinical disability was measured using the Expanded Disability Status Scale (EDSS), Nine-Hole-Peg test (9-HPT), and 25 feet Timed Walk Test (TWT). Stepwise linear regression was performed to assess the association between MUCCA and MRI measures, and between MUCCA and clinical disability. RESULTS: MUCCA was smaller (mean 11.7%) in MS patients compared with healthy controls (72.56±9.82 and 82.24±7.80 mm2 respectively; p<0.001), most prominently in male patients. MUCCA was associated with normalized brain volume, and number of cervical cord lesions. MUCCA was independently associated with EDSS, TWT, and 9-HPT. CONCLUSION: MUCCA was reduced in MS patients compared with healthy controls. It provides a relevant marker for clinical disability in long-standing disease, independent of other MRI measures.

Dissociation between phase and power correlation networks in the human brain is driven by co-occurrent bursts.

Well-known haemodynamic resting-state networks are better mirrored in power correlation networks than phase coupling networks in electrophysiological data. However, what do these power correlation networks reflect? We address this long-outstanding question in neuroscience using rigorous mathematical analysis, biophysical simulations with ground truth and application of these mathematical concepts to empirical magnetoencephalography (MEG) data. Our mathematical derivations show that for two non-Gaussian electrophysiological signals, their power correlation depends on their coherence, cokurtosis and conjugate-coherence. Only coherence and cokurtosis contribute to power correlation networks in MEG data, but cokurtosis is less affected by artefactual signal leakage and better mirrors haemodynamic resting-state networks. Simulations and MEG data show that cokurtosis may reflect co-occurrent bursting events. Our findings shed light on the origin of the complementary nature of power correlation networks to phase coupling networks and suggests that the origin of resting-state networks is partly reflected in co-occurent bursts in neuronal activity.

Preservation of thalamocortical circuitry is essential for good recovery after cardiac arrest.

Continuous electroencephalographam (EEG) monitoring contributes to prediction of neurological outcome in comatose cardiac arrest survivors. While the phenomenology of EEG abnormalities in postanoxic encephalopathy is well known, the pathophysiology, especially the presumed role of selective synaptic failure, is less understood. To further this understanding, we estimate biophysical model parameters from the EEG power spectra from individual patients with a good or poor recovery from a postanoxic encephalopathy. This biophysical model includes intracortical, intrathalamic, and corticothalamic synaptic strengths, as well as synaptic time constants and axonal conduction delays. We used continuous EEG measurements from hundred comatose patients recorded during the first 48 h postcardiac arrest, 50 with a poor neurological outcome [cerebral performance category ( CPC = 5 ) ] and 50 with a good neurological outcome ( CPC = 1 ). We only included patients that developed (dis-)continuous EEG activity within 48 h postcardiac arrest. For patients with a good outcome, we observed an initial relative excitation in the corticothalamic loop and corticothalamic propagation that subsequently evolved towards values observed in healthy controls. For patients with a poor outcome, we observed an initial increase in the cortical excitation-inhibition ratio, increased relative inhibition in the corticothalamic loop, delayed corticothalamic propagation of neuronal activity, and severely prolonged synaptic time constants that did not return to physiological values. We conclude that the abnormal EEG evolution in patients with a poor neurological recovery after cardiac arrest may result from persistent and selective synaptic failure that includes corticothalamic circuitry and also delayed corticothalamic propagation.

Early EEG monitoring predicts clinical outcome in patients with moderate to severe traumatic brain injury.

There is a need for reliable predictors in patients with moderate to severe traumatic brain injury to assist clinical decision making. We assess the ability of early continuous EEG monitoring at the intensive care unit (ICU) in patients with traumatic brain injury (TBI) to predict long term clinical outcome and evaluate its complementary value to current clinical standards. We performed continuous EEG measurements in patients with moderate to severe TBI during the first week of ICU admission. We assessed the Extended Glasgow Outcome Scale (GOSE) at 12 months, dichotomized into poor (GOSE 1-3) and good (GOSE 4-8) outcome. We extracted EEG spectral features, brain symmetry index, coherence, aperiodic exponent of the power spectrum, long range temporal correlations, and broken detailed balance. A random forest classifier using feature selection was trained to predict poor clinical outcome based on EEG features at 12, 24, 48, 72 and 96 h after trauma. We compared our predictor with the IMPACT score, the best available predictor, based on clinical, radiological and laboratory findings. In addition we created a combined model using EEG as well as the clinical, radiological and laboratory findings. We included hundred-seven patients. The best prediction model using EEG parameters was found at 72 h after trauma with an AUC of 0.82 (0.69-0.92), specificity of 0.83 (0.67-0.99) and sensitivity of 0.74 (0.63-0.93). The IMPACT score predicted poor outcome with an AUC of 0.81 (0.62-0.93), sensitivity of 0.86 (0.74-0.96) and specificity of 0.70 (0.43-0.83). A model using EEG and clinical, radiological and laboratory parameters resulted in a better prediction of poor outcome (p < 0.001) with an AUC of 0.89 (0.72-0.99), sensitivity of 0.83 (0.62-0.93) and specificity of 0.85 (0.75-1.00). EEG features have potential use for predicting clinical outcome and decision making in patients with moderate to severe TBI and provide complementary information to current clinical standards.

The Effect of Ketamine on Electrophysiological Connectivity in Major Depressive Disorder.

Major depressive disorder (MDD) is highly prevalent and frequently disabling. Only about 30% of patients respond to a first-line antidepressant treatment, and around 30% of patients are classified as "treatment-resistant" after failing to respond to multiple adequate trials. While most antidepressants target monoaminergic targets, ketamine is an N-methyl-D-aspartate (NMDA) antagonist that has shown rapid antidepressant effects when delivered intravenously or intranasally. While there is evidence that ketamine exerts its effects via enhanced α-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid (AMPA) throughput, its mechanism for relieving depressive symptoms is largely unknown. This study acquired resting-state magnetoencephalography (MEG) recordings after both ketamine and placebo infusions and investigated functional connectivity using a multilayer amplitude-amplitude correlation technique spanning the canonical frequency bands. Twenty-four healthy volunteers (HVs) and 27 unmedicated participants with MDD took part in a double-blind, placebo-controlled, crossover trial of 0.5 mg/kg IV ketamine. Order of infusion was randomized, and participants crossed over to receive the second infusion after two weeks. The results indicated widespread ketamine-induced reductions in connectivity in the alpha and beta bands that did not correlate with magnitude of antidepressant response. In contrast, the magnitude of ketamine's antidepressant effects in MDD participants was associated with cross-frequency connectivity for delta-alpha and delta-gamma bands, with HVs and ketamine non-responders showing connectivity decreases post-ketamine and ketamine responders demonstrating small increases in connectivity. These results may indicate functional subtypes of MDD and also suggest that neural responses to ketamine are fundamentally different between responders and non-responders.

Multilayer MEG functional connectivity as a potential marker for suicidal thoughts in major depressive disorder.

Major depressive disorder (MDD) is highly heterogeneous in its clinical presentation. The present exploratory study used magnetoencephalography (MEG) to investigate electrophysiological intrinsic connectivity differences between healthy volunteers and unmedicated participants with treatment-resistant MDD. The study examined canonical frequency bands from delta through gamma. In addition to group comparisons, correlational studies were conducted to determine whether connectivity was related to five symptom factors: depressed mood, tension, negative cognition, suicidal thoughts, and amotivation. The MDD and healthy volunteer groups did not differ significantly at baseline when corrected across all frequencies and clusters, although evidence of generalized slowing in MDD was observed. Notably, however, electrophysiological connectivity was strongly related to suicidal thoughts, particularly as coupling of low frequency power fluctuations (delta and theta) with alpha and beta power. This analysis revealed hub areas underlying this symptom cluster, including left hippocampus, left anterior insula, and bilateral dorsolateral prefrontal cortex. No other symptom cluster demonstrated a relationship with neurophysiological connectivity, suggesting a specificity to these results as markers of suicidal ideation.

Age-related differences in myeloarchitecture measured at 7 T.

We have used the magnetisation transfer (MT) MRI measure as a primary measure of myelination in both the gray matter (GM) of the 78 cortical automated anatomical labeling (AAL) regions of the brain, and the underlying white matter in each region, in a cohort of healthy adults (aged 19-62 year old). The results revealed a significant quadratic trend in myelination with age, with average global myelination peaking at 42.9 year old in gray matter, and at 41.7 year old in white matter. We also explored the possibility of using the Nuclear Overhauser Enhancement (NOE) effect, which is acquired in a similar method to MT, as an additional measure of myelination. We found that the MT and NOE signals were strongly correlated in the brain and that the NOE effects displayed similar (albeit weaker) parabolic trends with age. We also investigated differences in cortical thickness with age, and confirmed a previous result of a linear decline of 4.5 ± 1.2 µm/y.