A comparison of in vivo MRI based cortical myelin mapping using T1w/T2w and R1 mapping at 3T

Air pollution and cortical myelin T1w/T2w ratio estimates in school-age children from the ABCD and NeuroSmog studies

Air pollution affects human health and may disrupt brain maturation, including axon myelination, critical for efficient neural signaling. Here, we assess the impact of prenatal and current long-term particulate matter (PM) and nitrogen dioxide (NO2) exposure on cortical T1w/T2w ratios - a proxy for myelin content - in school-age children from the Adolescent Brain Cognitive Development (ABCD) Study (United States; N = 2021) and NeuroSmog study (Poland; N = 577), using Siemens scanners. Across both samples, we found that NO2 and PM were not significantly associated with cortical T1w/T2w except for one association of PM10 with lower T1w/T2w in the precuneus in NeuroSmog. Superficially, ABCD Study analyses including data from all scanner types (Siemens, GE, Philips; N = 3089) revealed a negative association between NO₂ exposure and T1w/T2w ratios. However, this finding could be an artifact of between-site sociodemographic differences and large scanner-type-related measurement differences. While significant associations between air pollution and cortical myelin were largely absent, these findings do not rule out the possibility that air pollution affects cortical myelin during other exposure periods/stages of neurodevelopment. Future research should examine these relationships across diverse populations and developmental periods using unified analysis methods to better understand the potential neurotoxic effects of air pollution.

Approximating R1 and R2: A Quantitative Approach to Clinical Weighted MRI

Weighted MRI images are widely used in clinical as well as open-source neuroimaging databases. Weighted images such as T1-weighted, T2-weighted, and proton density-weighted (T1w, T2w, and PDw, respectively) are used for evaluating the brain's macrostructure; however, their values cannot be used for microstructural analysis, as they lack physical meaning. Quantitative MRI (qMRI) relaxation rate parameters (e.g., R1 and R2) do contain microstructural physical meaning. Nevertheless, qMRI is rarely done in large-scale clinical databases. Currently, the weighted images ratio T1w/T2w is used as a quantifier to approximate the brain's microstructure. In this paper, we test three additional quantifiers that approximate quantitative maps, which can help bring quantitative MRI to the clinic for easy use. Following the signal equations and using simple mathematical operations, we combine the T1w, T2w, and PDw images to estimate the R1 and R2 maps. We find that two of these quantifiers (T1w/PDw and T1w/ln(T2w)) can approximate R1, and that (ln(T2w/PDw)) can approximate R2, in 3 datasets that were tested. We find that this approach also can be applied to T2w scans taken from widely available DTI (Diffusion Tensor Imaging) datasets. We tested these quantifiers on both in vitro phantom and in vivo human datasets. We found that the quantifiers accurately represent the quantitative parameters across datasets. Finally, we tested the quantifiers within a clinical context, and found that they are robust across datasets. Our work provides a simple pipeline to enhance the usability and quantitative accuracy of MRI weighted images.

Radiomics score derived from T1-w/T2-w ratio image can predict motor symptom progression in Parkinson's disease

OBJECTIVES

To clarify the association between a radiomics score (Rad-score) derived from T1-weighted signal intensity to T2-weighted signal intensity (T1-w/T2-w) ratio images and the progression of motor symptoms in Parkinson's disease (PD).

MATERIALS AND METHODS

This retrospective study included patients with PD enrolled in the Parkinson's Progression Markers Initiative. The Movement Disorders Society-Unified Parkinson's Disease Rating Scale Part III score ≥ 33 and/or Hoehn and Yahr stage ≥ 3 indicated motor function decline. The Rad-score was constructed using radiomics features extracted from T1-w/T2-w ratio images. The Kaplan-Meier analysis and Cox regression analyses were used to assess the time differences in motor function decline between the high and low Rad-score groups.

RESULTS

A total of 171 patients with PD were divided into training (n = 101, mean age at baseline, 61.6 ± 9.3 years) and testing (n = 70, mean age at baseline, 61.6 ± 10 years). The patients in the high Rad-score group had a shorter time to motor function decline than those in the low Rad-score group in the training dataset (log-rank test, p < 0.001) and testing dataset (log-rank test, p < 0.001). The multivariate Cox regression using the Rad-score and clinical factors revealed a significant association between the Rad-score and motor function decline in the training dataset (HR = 2.368, 95%CI:1.423-3.943, p < 0.001) and testing dataset (HR = 2.931, 95%CI:1.472-5.837, p = 0.002).

CONCLUSION

Rad-scores based on radiomics features derived from T1-w/T2-w ratio images were associated with the progression of motor symptoms in PD.

CLINICAL RELEVANCE STATEMENT

The radiomics score derived from the T1-weighted/T2-weighted ratio images offers a predictive tool for assessing the progression of motor symptom in patients with PD.

KEY POINTS

Radiomics score derived from T1-weighted/T2-weighted ratio images is correlated with the motor symptoms of Parkinson's disease. A high radiomics score correlated with faster motor function decline in patients with Parkinson's disease. The proposed radiomics score offers predictive insight into the progression of motor symptoms of Parkinson's disease.

Microstructural asymmetry in the human cortex

The human cerebral cortex shows hemispheric asymmetry, yet the microstructural basis of this asymmetry remains incompletely understood. Here, we probe layer-specific microstructural asymmetry using one post-mortem male brain. Overall, anterior and posterior regions show leftward and rightward asymmetry respectively, but this pattern varies across cortical layers. A similar anterior-posterior pattern is observed using in vivo Human Connectome Project (N = 1101) T1w/T2w microstructural data, with average cortical asymmetry showing the strongest similarity with post-mortem-based asymmetry of layer III. Moreover, microstructural asymmetry is found to be heritable, varies as a function of age and sex, and corresponds to intrinsic functional asymmetry. We also observe a differential association of language and markers of mental health with microstructural asymmetry patterns at the individual level, illustrating a functional divergence between inferior-superior and anterior-posterior microstructural axes, possibly anchored in development. Last, we could show concordant evidence with alternative in vivo microstructural measures: magnetization transfer (N = 286) and quantitative T1 (N = 50). Together, our study highlights microstructural asymmetry in the human cortex and its functional and behavioral relevance.

Assessing regional intracortical myelination in schizophrenia spectrum and bipolar disorders using the optimized T1w/T2w-ratio

BACKGROUND

Dysmyelination could be part of the pathophysiology of schizophrenia spectrum (SCZ) and bipolar disorders (BPD), yet few studies have examined myelination of the cerebral cortex. The ratio of T1- and T2-weighted magnetic resonance images (MRI) correlates with intracortical myelin. We investigated the T1w/T2w-ratio and its age trajectories in patients and healthy controls (CTR) and explored associations with antipsychotic medication use and psychotic symptoms.

METHODS

Patients with SCZ (n = 64; mean age = 30.4 years, s.d. = 9.8), BPD (n = 91; mean age 31.0 years, s.d. = 10.2), and CTR (n = 155; mean age = 31.9 years, s.d. = 9.1) who participated in the TOP study (NORMENT, University of Oslo, Norway) were clinically assessed and scanned using a General Electric 3 T MRI system. T1w/T2w-ratio images were computed using an optimized pipeline with intensity normalization and field inhomogeneity correction. Vertex-wise regression models were used to compare groups and examine group × age interactions. In regions showing significant differences, we explored associations with antipsychotic medication use and psychotic symptoms.

RESULTS

No main effect of diagnosis was found. However, age slopes of the T1w/T2w-ratio differed significantly between SCZ and CTR, predominantly in frontal and temporal lobe regions: Lower T1w/T2w-ratio values with higher age were found in CTR, but not in SCZ. Follow-up analyses revealed a more positive age slope in patients who were using antipsychotics and patients using higher chlorpromazine-equivalent doses.

CONCLUSIONS

While we found no evidence of reduced intracortical myelin in SCZ or BPD relative to CTR, different regional age trajectories in SCZ may suggest a promyelinating effect of antipsychotic medication.

The Advanced BRain Imaging on ageing and Memory (ABRIM) data collection: Study design, data processing, and rationale

To understand the neurocognitive mechanisms that underlie heterogeneity in cognitive ageing, recent scientific efforts have led to a growing public availability of imaging cohort data. The Advanced BRain Imaging on ageing and Memory (ABRIM) project aims to add to these existing datasets by taking an adult lifespan approach to provide a cross-sectional, normative database with a particular focus on connectivity, myelinization and iron content of the brain in concurrence with cognitive functioning, mechanisms of reserve, and sleep-wake rhythms. ABRIM freely shares MRI and behavioural data from 295 participants between 18-80 years, stratified by age decade and sex (median age 52, IQR 36-66, 53.20% females). The ABRIM MRI collection consists of both the raw and pre-processed structural and functional MRI data to facilitate data usage among both expert and non-expert users. The ABRIM behavioural collection includes measures of cognitive functioning (i.e., global cognition, processing speed, executive functions, and memory), proxy measures of cognitive reserve (e.g., educational attainment, verbal intelligence, and occupational complexity), and various self-reported questionnaires (e.g., on depressive symptoms, pain, and the use of memory strategies in daily life and during a memory task). In a sub-sample (n = 120), we recorded sleep-wake rhythms using an actigraphy device (Actiwatch 2, Philips Respironics) for a period of 7 consecutive days. Here, we provide an in-depth description of our study protocol, pre-processing pipelines, and data availability. ABRIM provides a cross-sectional database on healthy participants throughout the adult lifespan, including numerous parameters relevant to improve our understanding of cognitive ageing. Therefore, ABRIM enables researchers to model the advanced imaging parameters and cognitive topologies as a function of age, identify the normal range of values of such parameters, and to further investigate the diverse mechanisms of reserve and resilience.

Intensity scaling of conventional brain magnetic resonance images avoiding cerebral reference regions: A systematic review

BACKGROUND

Conventional brain magnetic resonance imaging (MRI) produces image intensities that have an arbitrary scale, hampering quantification. Intensity scaling aims to overcome this shortfall. As neurodegenerative and inflammatory disorders may affect all brain compartments, reference regions within the brain may be misleading. Here we summarize approaches for intensity scaling of conventional T1-weighted (w) and T2w brain MRI avoiding reference regions within the brain.

METHODS

Literature was searched in the databases of Scopus, PubMed, and Web of Science. We included only studies that avoided reference regions within the brain for intensity scaling and provided validating evidence, which we divided into four categories: 1) comparative variance reduction, 2) comparative correlation with clinical parameters, 3) relation to quantitative imaging, or 4) relation to histology.

RESULTS

Of the 3825 studies screened, 24 fulfilled the inclusion criteria. Three studies used scaled T1w images, 2 scaled T2w images, and 21 T1w/T2w-ratio calculation (with double counts). A robust reduction in variance was reported. Twenty studies investigated the relation of scaled intensities to different types of quantitative imaging. Statistically significant correlations with clinical or demographic data were reported in 8 studies. Four studies reporting the relation to histology gave no clear picture of the main signal driver of conventional T1w and T2w MRI sequences.

CONCLUSIONS

T1w/T2w-ratio calculation was applied most often. Variance reduction and correlations with other measures suggest a biologically meaningful signal harmonization. However, there are open methodological questions and uncertainty on its biological underpinning. Validation evidence on other scaling methods is even sparser.

Changes in microstructural similarity of hippocampal subfield circuits in pathological cognitive aging

The hippocampal networks support multiple cognitive functions and may have biological roles and functions in pathological cognitive aging (PCA) and its associated diseases, which have not been explored. In the current study, a total of 116 older adults with 39 normal controls (NC) (mean age: 52.3 ± 13.64 years; 16 females), 39 mild cognitive impairment (MCI) (mean age: 68.15 ± 9.28 years, 14 females), and 38 dementia (mean age: 73.82 ± 8.06 years, 8 females) were included. The within-hippocampal subfields and the cortico-hippocampal circuits were assessed via a micro-structural similarity network approach using T1w/T2w ratio and regional gray matter tissue probability maps, respectively. An analysis of covariance was conducted to identify between-group differences in structural similarities among hippocampal subfields. The partial correlation analyses were performed to associate changes in micro-structural similarities with cognitive performance in the three groups, controlling the effect of age, sex, education, and cerebral small-vessel disease. Compared with the NC, an altered T1w/T2w ratio similarity between left CA3 and left subiculum was observed in the mild cognitive impairment (MCI) and dementia. The left CA3 was the most impaired region correlated with deteriorated cognitive performance. Using these regions as seeds for GM similarity comparisons between hippocampal subfields and cortical regions, group differences were observed primarily between the left subiculum and several cortical regions. By utilizing T1w/T2w ratio as a proxy measure for myelin content, our data suggest that the imbalanced synaptic weights within hippocampal CA3 provide a substrate to explain the abnormal firing characteristics of hippocampal neurons in PCA. Furthermore, our work depicts specific brain structural characteristics of normal and pathological cognitive aging and suggests a potential mechanism for cognitive aging heterogeneity.

An fMRI Dataset on Social Reward Processing and Decision Making in Younger and Older Adults

Behavioural and neuroimaging research has shown that older adults are less sensitive to financial losses compared to younger adults. Yet relatively less is known about age-related differences in social decisions and social reward processing. As part of a pilot study, we collected behavioural and functional magnetic resonance imaging (fMRI) data from 50 participants (Younger: N = 26, ages 18-34 years; Older: N = 24, ages 63-80 years) who completed three tasks in the scanner: an economic trust game as the investor with three partners (computer, stranger, friend) as the investee; a card-guessing task with monetary gains and losses shared with three partners (computer, stranger, friend); and an ultimatum game as responder to three anonymous proposers (computer, age-similar adults, age-dissimilar adults). We also collected B0 field maps and high-resolution structural images (T1-weighted and T2-weighted images). These data could be reused to answer questions about moment-to-moment variability in fMRI signal, representational similarity between tasks, and brain structure.

Regional cortical thinning, demyelination and iron loss in cerebral small vessel disease

The link between white matter hyperintensities (WMH) and cortical thinning is thought to be an important pathway by which WMH contributes to cognitive deficits in cerebral small vessel disease (SVD). However, the mechanism behind this association and the underlying tissue composition abnormalities are unclear. The objective of this study is to determine the association between WMH and cortical thickness, and the in vivo tissue composition abnormalities in the WMH-connected cortical regions. In this cross-sectional study, we included 213 participants with SVD who underwent standardized protocol including multimodal neuroimaging scans and cognitive assessment (i.e. processing speed, executive function and memory). We identified the cortex connected to WMH using probabilistic tractography starting from the WMH and defined the WMH-connected regions at three connectivity levels (low, medium and high connectivity level). We calculated the cortical thickness, myelin and iron of the cortex based on T1-weighted, quantitative R1, R2* and susceptibility maps. We used diffusion-weighted imaging to estimate the mean diffusivity of the connecting white matter tracts. We found that cortical thickness, R1, R2* and susceptibility values in the WMH-connected regions were significantly lower than in the WMH-unconnected regions (all Pcorrected < 0.001). Linear regression analyses showed that higher mean diffusivity of the connecting white matter tracts were related to lower thickness (β = -0.30, Pcorrected < 0.001), lower R1 (β = -0.26, Pcorrected = 0.001), lower R2* (β = -0.32, Pcorrected < 0.001) and lower susceptibility values (β = -0.39, Pcorrected < 0.001) of WMH-connected cortical regions at high connectivity level. In addition, lower scores on processing speed were significantly related to lower cortical thickness (β = 0.20, Pcorrected = 0.030), lower R1 values (β = 0.20, Pcorrected = 0.006), lower R2* values (β = 0.29, Pcorrected = 0.006) and lower susceptibility values (β = 0.19, Pcorrected = 0.024) of the WMH-connected regions at high connectivity level, independent of WMH volumes and the cortical measures of WMH-unconnected regions. Together, our study demonstrated that the microstructural integrity of white matter tracts passing through WMH is related to the regional cortical abnormalities as measured by thickness, R1, R2* and susceptibility values in the connected cortical regions. These findings are indicative of cortical thinning, demyelination and iron loss in the cortex, which is most likely through the disruption of the connecting white matter tracts and may contribute to processing speed impairment in SVD, a key clinical feature of SVD. These findings may have implications for finding intervention targets for the treatment of cognitive impairment in SVD by preventing secondary degeneration.

Neuroanatomical and Functional Dissociations between Variably Present Anterior Lateral Prefrontal Sulci

The lateral prefrontal cortex (LPFC) is an evolutionarily expanded region in humans that is critical for numerous complex functions, many of which are largely hominoid specific. Although recent work shows that the presence or absence of specific sulci in anterior LPFC is associated with cognitive performance across age groups, it is unknown whether the presence of these structures relates to individual differences in the functional organization of LPFC. To fill this gap in knowledge, we leveraged multimodal neuroimaging data from two samples encompassing 82 young adult humans (aged 22-36 years) and show that the dorsal and ventral components of the paraintermediate frontal sulcus, or pimfs, present distinct morphological (surface area), architectural (thickness and myelination), and functional (resting-state connectivity networks) properties. We further contextualize the pimfs components within classic and modern cortical parcellations. Taken together, the dorsal and ventral pimfs components mark transitions in LPFC anatomy and function, across metrics and parcellations. These results emphasize that the pimfs is a critical structure to consider when examining individual differences in the anatomical and functional organization of LPFC and suggest that future individual-level parcellations could benefit from incorporating sulcal anatomy when delineating LPFC cortical regions.

Family income is not significantly associated with T1w/T2w ratio in the Human Connectome Project in Development

Growing evidence indicates that brain development varies as a function of family socioeconomic status (SES). Numerous studies have demonstrated that children from low-SES backgrounds have thinner cortex than children from higher-SES backgrounds. A recent study in a large developmental sample found widespread associations between lower SES and greater cortical T1w/T2w ratio-thought to be an indirect proxy for cortical myelin. We evaluated the association of family income with cortical T1w/T2w ratio as a function of age in the Human Connectome Project in Development sample of 989 youth aged 8-21 years. We observed no associations between family income and T1w/T2w ratio that were significant after corrections for multiple comparisons at the region, network, or whole-brain level. Region of practical equivalence (ROPE) analyses were also consistent with the absence of an association between family income and T1w/T2w ratio. We discuss potential methodological sources of inconsistency between this and the previous study examining the same question. While the question of whether family income may influence cortical myelin development remains, these null results may indicate that the association between SES and cortical myelin development may not be as strong as with other aspects of brain structure.

Human white matter myelinates faster in utero than ex utero

The formation of myelin, the fatty sheath that insulates nerve fibers, is critical for healthy brain function. A fundamental open question is what impact being born has on myelin growth. To address this, we evaluated a large (n = 300) cross-sectional sample of newborns from the Developing Human Connectome Project (dHCP). First, we developed software for the automated identification of 20 white matter bundles in individual newborns that is well suited for large samples. Next, we fit linear models that quantify how T1w/T2w (a myelin-sensitive imaging contrast) changes over time at each point along the bundles. We found faster growth of T1w/T2w along the lengths of all bundles before birth than right after birth. Further, in a separate longitudinal sample of preterm infants (N = 34), we found lower T1w/T2w than in full-term peers measured at the same age. By applying the linear models fit on the cross-section sample to the longitudinal sample of preterm infants, we find that their delay in T1w/T2w growth is well explained by the amount of time they spent developing in utero and ex utero. These results suggest that white matter myelinates faster in utero than ex utero. The reduced rate of myelin growth after birth, in turn, explains lower myelin content in individuals born preterm and could account for long-term cognitive, neurological, and developmental consequences of preterm birth. We hypothesize that closely matching the environment of infants born preterm to what they would have experienced in the womb may reduce delays in myelin growth and hence improve developmental outcomes.

Resistance Training Maintains White Matter and Physical Function in Older Women with Cerebral Small Vessel Disease: An Exploratory Analysis of a Randomized Controlled Trial

Background

As the aging population grows, there is an increasing need to develop accessible interventions against risk factors for cognitive impairment and dementia, such as cerebral small vessel disease (CSVD). The progression of white matter hyperintensities (WMHs), a key hallmark of CSVD, can be slowed by resistance training (RT). We hypothesize RT preserves white matter integrity and that this preservation is associated with improved cognitive and physical function.

Objective

To determine if RT preserves regional white matter integrity and if any changes are associated with cognitive and physical outcomes.

Methods

Using magnetic resonance imaging data from a 12-month randomized controlled trial, we compared the effects of a twice-weekly 60-minute RT intervention versus active control on T1-weighted over T2-weighted ratio (T1w/T2w; a non-invasive proxy measure of white matter integrity) in a subset of study participants (N = 21 females, mean age = 69.7 years). We also examined the association between changes in T1w/T2w with two key outcomes of the parent study: (1) selective attention and conflict resolution, and (2) peak muscle power.

Results

Compared with an active control group, RT increased T1w/T2w in the external capsule (p = 0.024) and posterior thalamic radiations (p = 0.013) to a greater degree. Increased T1w/T2w in the external capsule was associated with an increase in peak muscle power (p = 0.043) in the RT group.

Conclusion

By maintaining white matter integrity, RT may be a promising intervention to counteract the pathological changes that accompany CSVD, while improving functional outcomes such as muscle power.

Neuroanatomical and functional dissociations between variably present anterior lateral prefrontal sulci

The lateral prefrontal cortex (LPFC) is an evolutionarily expanded region in humans that is critical for numerous complex functions, many of which are largely hominoid-specific. While recent work shows that the presence or absence of specific sulci in anterior LPFC is associated with cognitive performance across age groups, it is unknown whether the presence of these structures relates to individual differences in the functional organization of LPFC. To fill this gap in knowledge, we leveraged multimodal neuroimaging data from 72 young adult humans aged 22-36 and show that dorsal and ventral components of the paraintermediate frontal sulcus (pimfs) present distinct morphological (surface area), architectural (thickness and myelination), and functional (resting-state connectivity networks) properties. We further contextualize the pimfs components within classic and modern cortical parcellations. Taken together, the dorsal and ventral pimfs components mark transitions in anatomy and function in LPFC, across metrics and parcellations. These results emphasize that the pimfs is a critical structure to consider when examining individual differences in the anatomical and functional organization of LPFC and highlight the importance of considering individual anatomy when investigating structural and functional features of the cortex.

The intracortical myelin content of impulsive choices: results from T1- and T2-weighted MRI myelin mapping

Delay discounting (DD) refers to a phenomenon that humans tend to choose small-sooner over large-later rewards during intertemporal choices. Steep discounting of delayed outcome is related to a variety of maladaptive behaviors and is considered as a transdiagnostic process across psychiatric disorders. Previous studies have investigated the association between brain structure (e.g. gray matter volume) and DD; however, it is unclear whether the intracortical myelin (ICM) influences DD. Here, based on a sample of 951 healthy young adults drawn from the Human Connectome Project, we examined the relationship between ICM, which was measured by the contrast of T1w and T2w images, and DD and further tested whether the identified associations were mediated by the regional homogeneity (ReHo) of brain spontaneous activity. Vertex-wise regression analyses revealed that steeper DD was significantly associated with lower ICM in the left temporoparietal junction (TPJ) and right middle-posterior cingulate cortex. Region-of-interest analysis revealed that the ReHo values in the left TPJ partially mediated the association of its myelin content with DD. Our findings provide the first evidence that cortical myelination is linked with individual differences in decision impulsivity and suggest that the myelin content affects cognitive performances partially through altered local brain synchrony.

An attention-based context-informed deep framework for infant brain subcortical segmentation

Precise segmentation of subcortical structures from infant brain magnetic resonance (MR) images plays an essential role in studying early subcortical structural and functional developmental patterns and diagnosis of related brain disorders. However, due to the dynamic appearance changes, low tissue contrast, and tiny subcortical size in infant brain MR images, infant subcortical segmentation is a challenging task. In this paper, we propose a context-guided, attention-based, coarse-to-fine deep framework to precisely segment the infant subcortical structures. At the coarse stage, we aim to directly predict the signed distance maps (SDMs) from multi-modal intensity images, including T1w, T2w, and the ratio of T1w and T2w images, with an SDM-Unet, which can leverage the spatial context information, including the structural position information and the shape information of the target structure, to generate high-quality SDMs. At the fine stage, the predicted SDMs, which encode spatial-context information of each subcortical structure, are integrated with the multi-modal intensity images as the input to a multi-source and multi-path attention Unet (M2A-Unet) for achieving refined segmentation. Both the 3D spatial and channel attention blocks are added to guide the M2A-Unet to focus more on the important subregions and channels. We additionally incorporate the inner and outer subcortical boundaries as extra labels to help precisely estimate the ambiguous boundaries. We validate our method on an infant MR image dataset and on an unrelated neonatal MR image dataset. Compared to eleven state-of-the-art methods, the proposed framework consistently achieves higher segmentation accuracy in both qualitative and quantitative evaluations of infant MR images and also exhibits good generalizability in the neonatal dataset.

The Potential of Myelin-Sensitive Imaging: Redefining Spatiotemporal Patterns of Myeloarchitecture

Recent advances in magnetic resonance imaging (MRI) have paved the way for approximation of myelin content in vivo. In this review, our main goal was to determine how to best capitalize on myelin-sensitive imaging. First, we briefly overview the theoretical and empirical basis for the myelin sensitivity of different MRI markers and, in doing so, highlight how multimodal imaging approaches are important for enhancing specificity to myelin. Then, we discuss recent studies that have probed the nonuniform distribution of myelin across cortical layers and along white matter tracts. These approaches, collectively known as myelin profiling, have provided detailed depictions of myeloarchitecture in both the postmortem and living human brain. Notably, MRI-based profiling studies have recently focused on investigating whether it can capture interindividual variability in myelin characteristics as well as trajectories across the lifespan. Finally, another line of recent evidence emphasizes the contribution of region-specific myelination to large-scale organization, demonstrating the impact of myelination on global brain networks. In conclusion, we suggest that combining well-validated MRI markers with profiling techniques holds strong potential to elucidate individual differences in myeloarchitecture, which has important implications for understanding brain function and disease.

Macromolecular tissue volume mapping of lateral geniculate nucleus subdivisions in living human brains

The lateral geniculate nucleus (LGN) is a key thalamic nucleus in the visual system, which has an important function in relaying retinal visual input to the visual cortex. The human LGN is composed mainly of magnocellular (M) and parvocellular (P) subdivisions, each of which has different stimulus selectivity in neural response properties. Previous studies have discussed the potential relationship between LGN subdivisions and visual disorders based on psychophysical data on specific types of visual stimuli. However, these relationships remain speculative because non-invasive measurements of these subdivisions are difficult due to the small size of the LGN. Here we propose a method to identify these subdivisions by combining two structural MR measures: high-resolution proton-density weighted images and macromolecular tissue volume (MTV) maps. We defined the M and P subdivisions based on MTV fraction data and tested the validity of the definition by (1) comparing the data with that from human histological studies, (2) comparing the data with functional magnetic resonance imaging measurements on stimulus selectivity, and (3) analyzing the test-retest reliability. The findings demonstrated that the spatial organization of the M and P subdivisions was consistent across subjects and in line with LGN subdivisions observed in human histological data. Moreover, the difference in stimulus selectivity between the subdivisions identified using MTV was consistent with previous physiology literature. The definition of the subdivisions based on MTV was shown to be robust over measurements taken on different days. These results suggest that MTV mapping is a promising approach for evaluating the tissue properties of LGN subdivisions in living humans. This method potentially will enable neuroscientific and clinical hypotheses about the human LGN subdivisions to be tested.

Effects of non-modifiable risk factors of Alzheimer's disease on intracortical myelin content

BACKGROUND

Non-modifiable risk factors of Alzheimer's disease (AD) have lifelong effects on cortical integrity that could be mitigated if identified at early stages. However, it remains unknown whether cortical microstructure is affected in older individuals with non-modifiable AD risk factors and whether altered cortical tissue integrity produces abnormalities in brain functional networks in this AD-risk population.

METHODS

Using relative T1w/T2w (rT1w/T2w) ratio maps, we have compared tissue integrity of normal-appearing cortical GM between controls and cognitively normal older adults with either APOE4 (N = 50), with a first-degree family history (FH) of AD (N = 52), or with the co-occurrence of both AD risk factors (APOE4+FH) (N = 35). Additionally, individuals with only one risk factor (APOE4 or FH) were combined into one group (N = 102) and compared with controls. The same number of controls matched in age, sex, and years of education was employed for each of these comparisons. Group differences in resting state functional connectivity (rs-FC) patterns were also investigated, using as FC seeds those cortical regions showing significant changes in rT1w/T2w ratios.

RESULTS

Overall, individuals with non-modifiable AD risk factors exhibited significant variations in rT1w/T2w ratios compared to controls, being APOE4 and APOE4+FH at opposite ends of a continuum. The co-occurrence of APOE4 and FH was further accompanied by altered patterns of rs-FC.

CONCLUSIONS

These findings may have practical implications for early detection of cortical abnormalities in older populations with APOE4 and/or FH of AD and open new avenues to monitor changes in cortical tissue integrity associated with non-modifiable AD risk factors.

Spatial profiles provide sensitive MRI measures of the midbrain micro- and macrostructure

The midbrain is the rostral-most part of the brainstem. It contains numerous nuclei and white matter tracts, which are involved in motor, auditory and visual processing, and changes in their structure and function have been associated with aging, as well as neurodegenerative disorders. Current tools for estimating midbrain subregions and their structure with MRI require high resolution and multi-parametric quantitative MRI measures. We propose an approach that relies on morphology to calculate profiles along the midbrain and show these profiles are sensitive to the underlying macrostructure of the midbrain. First, we show that the midbrain structure can be sampled, within subject space, along three main axes of the left and right midbrain, producing profiles that are similar across subjects. We use two data sets with different field strengths, that contain R1, R2* and QSM maps and show that the profiles are highly correlated both across subjects and between datasets. Next, we compare profiles of the midbrain that sample ROIs, and show that the profiles along the first two axes sample the midbrain in a way that reliably separates the main structures, i.e., the substantia nigra, the red nucleus, and periaqueductal gray. We further show that age differences which are localized to specific nuclei, are reflected in the profiles. Finally, we generalize the same approach to calculate midbrain profiles on a third clinically relevant dataset using HCP subjects, with metrics such as the diffusion tensor and semi-quantitative data such as T1w/T2w maps. Our results suggest that midbrain profiles, both of quantitative and semi-quantitative estimates are sensitive to the underlying macrostructure of the midbrain. The midbrain profiles are calculated in native space, and rely on simple measurements. We show that it is robust and can be easily expanded to different datasets, and as such we hope that it will be of great use to the community and to the study of the midbrain in particular.

Cortical myelin profile variations in healthy aging brain: A T1w/T2w ratio study

Demyelination is observed in both healthy aging and age-related neurodegenerative disorders. While the significance of myelin within the cortex is well acknowledged, studies focused on intracortical demyelination and depth-specific structural alterations in normal aging are lacking. Using the recently available Human Connectome Project Aging dataset, we investigated intracortical myelin in a normal aging population using the T1w/T2w ratio. To capture the fine changes across cortical depths, we employed a surface-based approach by constructing cortical profiles traveling perpendicularly through the cortical ribbon and sampling T1w/T2w values. The curvatures of T1w/T2w cortical profiles may be influenced by differences in local myeloarchitecture and other tissue properties, which are known to vary across cortical regions. To quantify the shape of these profiles, we parametrized the level of curvature using a nonlinearity index (NLI) that measures the deviation of the profile from a straight line. We showed that NLI exhibited a steep decline in aging that was independent of local cortical thinning. Further examination of the profiles revealed that lower T1w/T2w near the gray-white matter boundary and superficial cortical depths were major contributors to the apparent NLI variations with age. These findings suggest that demyelination and changes in other T1w/T2w related tissue properties in normal aging may be depth-specific and highlight the potential of NLI as a unique marker of microstructural alterations within the cerebral cortex.

An Open MRI Dataset For Multiscale Neuroscience

Multimodal neuroimaging grants a powerful window into the structure and function of the human brain at multiple scales. Recent methodological and conceptual advances have enabled investigations of the interplay between large-scale spatial trends (also referred to as gradients) in brain microstructure and connectivity, offering an integrative framework to study multiscale brain organization. Here, we share a multimodal MRI dataset for Microstructure-Informed Connectomics (MICA-MICs) acquired in 50 healthy adults (23 women; 29.54 ± 5.62 years) who underwent high-resolution T1-weighted MRI, myelin-sensitive quantitative T1 relaxometry, diffusion-weighted MRI, and resting-state functional MRI at 3 Tesla. In addition to raw anonymized MRI data, this release includes brain-wide connectomes derived from (i) resting-state functional imaging, (ii) diffusion tractography, (iii) microstructure covariance analysis, and (iv) geodesic cortical distance, gathered across multiple parcellation scales. Alongside, we share large-scale gradients estimated from each modality and parcellation scale. Our dataset will facilitate future research examining the coupling between brain microstructure, connectivity, and function. MICA-MICs is available on the Canadian Open Neuroscience Platform data portal ( https://portal.conp.ca ) and the Open Science Framework ( https://osf.io/j532r/ ).

Empirical transmit field bias correction of T1w/T2w myelin maps

T1-weighted divided by T2-weighted (T1w/T2w) myelin maps were initially developed for neuroanatomical analyses such as identifying cortical areas, but they are increasingly used in statistical comparisons across individuals and groups with other variables of interest. Existing T1w/T2w myelin maps contain radiofrequency transmit field (B1+) biases, which may be correlated with these variables of interest, leading to potentially spurious results. Here we propose two empirical methods for correcting these transmit field biases using either explicit measures of the transmit field or alternatively a 'pseudo-transmit' approach that is highly correlated with the transmit field at 3T. We find that the resulting corrected T1w/T2w myelin maps are both better neuroanatomical measures (e.g., for use in cross-species comparisons), and more appropriate for statistical comparisons of relative T1w/T2w differences across individuals and groups (e.g., sex, age, or body-mass-index) within a consistently acquired study at 3T. We recommend that investigators who use the T1w/T2w approach for mapping cortical myelin use these B1+ transmit field corrected myelin maps going forward.

Graded Variation in T1w/T2w Ratio during Adolescence: Measurement, Caveats, and Implications for Development of Cortical Myelin

Adolescence is characterized by the maturation of cortical microstructure and connectivity supporting complex cognition and behavior. Axonal myelination influences brain connectivity during development by enhancing neural signaling speed and inhibiting plasticity. However, the maturational timing of cortical myelination during human adolescence remains poorly understood. Here, we take advantage of recent advances in high-resolution cortical T1w/T2w mapping methods, including principled correction of B1+ transmit field effects, using data from the Human Connectome Project in Development (HCP-D; N = 628, ages 8-21). We characterize microstructural changes relevant to myelination by estimating age-related differences in T1w/T2w throughout the cerebral neocortex from childhood to early adulthood. We apply Bayesian spline models and clustering analysis to demonstrate graded variation in age-dependent cortical T1w/T2w differences that are correlated with the sensorimotor-association (S-A) axis of cortical organization reported by others. In sensorimotor areas, T1w/T2w ratio measures start at high levels at early ages, increase at a fast pace, and decelerate at later ages (18-21). In intermediate multimodal areas along the S-A axis, T1w/T2w starts at intermediate levels and increases linearly at an intermediate pace. In transmodal/paralimbic association areas, T1w/T2w starts at low levels and increases linearly at the slowest pace. These data provide evidence for graded variation of the T1w/T2w ratio along the S-A axis that may reflect cortical myelination changes during adolescence underlying the development of complex information processing and psychological functioning. We discuss the implications of these results as well as caveats in interpreting magnetic resonance imaging (MRI)-based estimates of myelination.SIGNIFICANCE STATEMENT Myelin is a lipid membrane that is essential to healthy brain function. Myelin wraps axons to increase neural signaling speed, enabling complex neuronal functioning underlying learning and cognition. Here, we characterize the developmental timing of myelination across the cerebral cortex during adolescence using a noninvasive proxy measure, T1w/T2w mapping. Our results provide new evidence demonstrating graded variation across the cortex in the timing of T1w/T2w changes during adolescence, with rapid T1w/T2w increases in lower-order sensory areas and gradual T1w/T2w increases in higher-order association areas. This spatial pattern of microstructural brain development closely parallels the sensorimotor-to-association axis of cortical organization and plasticity during ontogeny.

Mapping microstructural gradients of the human striatum in normal aging and Parkinson's disease

Mapping structural spatial change (i.e., gradients) in the striatum is essential for understanding the function of the basal ganglia in both health and disease. We developed a method to identify and quantify gradients of microstructure in the single human brain in vivo. We found spatial gradients in the putamen and caudate nucleus of the striatum that were robust across individuals, clinical conditions, and datasets. By exploiting multiparametric quantitative MRI, we found distinct, spatially dependent, aging-related alterations in water content and iron concentration. Furthermore, we found cortico-striatal microstructural covariation, showing relations between striatal structural gradients and cortical hierarchy. In Parkinson's disease (PD) patients, we found abnormal gradients in the putamen, revealing changes in the posterior putamen that explain patients' dopaminergic loss and motor dysfunction. Our work provides a noninvasive approach for studying the spatially varying, structure-function relationship in the striatum in vivo, in normal aging and PD.

T1/T2-weighted ratio in multiple sclerosis: A longitudinal study with clinical associations

•

Alterations in T1-w/T2-w ratio precede lesion formation in CIS patients.

•

Longitudinal decreases in T1-w/T2-w were associated with disease activity in CIS.

•

Lower T1-w/T2-w was associated with longer disease duration and higher EDSS in MS.

Background

T1w/T2-w ratio has been proposed as a clinically feasible MRI biomarker to assess tissue integrity in multiple sclerosis. However, no data is available in the earliest stages of the disease and longitudinal studies analysing clinical associations are scarce.

Objective

To describe longitudinal changes in T1-w/T2-w in patients with clinically isolated syndrome (CIS) and multiple sclerosis, and to investigate their clinical associations.

Methods

T1-w/T2-w images were generated and the mean value obtained in the corresponding lesion, normal-appearing grey (NAGM) and white matter (NAWM) masks. By co-registering baseline to follow-up MRI, evolved lesions were assessed; and by placing the mask of new lesions to the baseline study, the pre-lesional tissue integrity was measured.

Results

We included 171 CIS patients and 22 established multiple sclerosis patients. In CIS, evolved lesions showed significant T1-w/T2-w increases compared to baseline (+7.6%, P < 0.001). T1-w/T2-w values in new lesions were lower than in pre-lesional tissue (-28.2%, P < 0.001), and pre-lesional tissue was already lower than baseline NAWM (-7.8%, P < 0.001). In CIS at baseline, higher NAGM T1-w/T2-w was associated with multiple sclerosis diagnosis, and longitudinal decreases in NAGM and NAWM T1-w/T2-w were associated with disease activity. In established multiple sclerosis, T1-w/T2-w was inversely correlated with clinical disability and disease duration.

Conclusion

A decrease in T1-w/T2-w ratio precedes lesion formation. In CIS, higher T1-w/T2-w was associated with multiple sclerosis diagnosis. In established multiple sclerosis, lower T1-w/T2-w values were associated with clinical disability. The possible differential impact of chronic inflammation, iron deposition and demyelination should be considered to interpret these findings.

T1w/T2w Ratio and Cognition in 9-to-11-Year-Old Children

Childhood is a period of extensive cortical and neural development. Among other things, axons in the brain gradually become more myelinated, promoting the propagation of electrical signals between different parts of the brain, which in turn may facilitate skill development. Myelin is difficult to assess in vivo, and measurement techniques are only just beginning to make their way into standard imaging protocols in human cognitive neuroscience. An approach that has been proposed as an indirect measure of cortical myelin is the T1w/T2w ratio, a contrast that is based on the intensities of two standard structural magnetic resonance images. Although not initially intended as such, researchers have recently started to use the T1w/T2w contrast for between-subject comparisons of cortical data with various behavioral and cognitive indices. As a complement to these earlier findings, we computed individual cortical T1w/T2w maps using data from the Adolescent Brain Cognitive Development study (N = 960; 449 females; aged 8.9 to 11.0 years) and related the T1w/T2w maps to indices of cognitive ability; in contrast to previous work, we did not find significant relationships between T1w/T2w values and cognitive performance after correcting for multiple testing. These findings reinforce existent skepticism about the applicability of T1w/T2w ratio for inter-individual comparisons.

The Myelin-Weighted Connectome in Parkinson's Disease

BACKGROUND

Even though Parkinson's disease (PD) is typically viewed as largely affecting gray matter, there is growing evidence that there are also structural changes in the white matter. Traditional connectomics methods that study PD may not be specific to underlying microstructural changes, such as myelin loss.

OBJECTIVE

The primary objective of this study is to investigate the PD-induced changes in myelin content in the connections emerging from the basal ganglia and the brainstem. For the weighting of the connectome, we used the longitudinal relaxation rate as a biologically grounded myelin-sensitive metric.

METHODS

We computed the myelin-weighted connectome in 35 healthy control subjects and 81 patients with PD. We used partial least squares to highlight the differences between patients with PD and healthy control subjects. Then, a ring analysis was performed on selected brainstem and subcortical regions to evaluate each node's potential role as an epicenter for disease propagation. Then, we used behavioral partial least squares to relate the myelin alterations with clinical scores.

RESULTS

Most connections (~80%) emerging from the basal ganglia showed a reduced myelin content. The connections emerging from potential epicentral nodes (substantia nigra, nucleus basalis of Meynert, amygdala, hippocampus, and midbrain) showed significant decrease in the longitudinal relaxation rate (P < 0.05). This effect was not seen for the medulla and the pons.

CONCLUSIONS

The myelin-weighted connectome was able to identify alteration of the myelin content in PD in basal ganglia connections. This could provide a different view on the importance of myelination in neurodegeneration and disease progression. © 2021 The Authors. Movement Disorders published by Wiley Periodicals LLC on behalf of International Parkinson and Movement Disorder Society.

Multisite reproducibility and test-retest reliability of the T1w/T2w-ratio: A comparison of processing methods

BACKGROUND

The ratio of T1-weighted (T1w) and T2-weighted (T2w) magnetic resonance imaging (MRI) images is often used as a proxy measure of cortical myelin. However, the T1w/T2w-ratio is based on signal intensities that are inherently non-quantitative and known to be affected by extrinsic factors. To account for this a variety of processing methods have been proposed, but a systematic evaluation of their efficacy is lacking. Given the dependence of the T1w/T2w-ratio on scanner hardware and T1w and T2w protocols, it is important to ensure that processing pipelines perform well also across different sites.

METHODS

We assessed a variety of processing methods for computing cortical T1w/T2w-ratio maps, including correction methods for nonlinear field inhomogeneities, local outliers, and partial volume effects as well as intensity normalisation. These were implemented in 33 processing pipelines which were applied to four test-retest datasets, with a total of 170 pairs of T1w and T2w images acquired on four different MRI scanners. We assessed processing pipelines across datasets in terms of their reproducibility of expected regional distributions of cortical myelin, lateral intensity biases, and test-retest reliability regionally and across the cortex. Regional distributions were compared both qualitatively with histology and quantitatively with two reference datasets, YA-BC and YA-B1+, from the Human Connectome Project.

RESULTS

Reproducibility of raw T1w/T2w-ratio distributions was overall high with the exception of one dataset. For this dataset, Spearman rank correlations increased from 0.27 to 0.70 after N3 bias correction relative to the YA-BC reference and from -0.04 to 0.66 after N4ITK bias correction relative to the YA-B1+ reference. Partial volume and outlier corrections had only marginal effects on the reproducibility of T1w/T2w-ratio maps and test-retest reliability. Before intensity normalisation, we found large coefficients of variation (CVs) and low intraclass correlation coefficients (ICCs), with total whole-cortex CV of 10.13% and whole-cortex ICC of 0.58 for the raw T1w/T2w-ratio. Intensity normalisation with WhiteStripe, RAVEL, and Z-Score improved total whole-cortex CVs to 5.91%, 5.68%, and 5.19% respectively, whereas Z-Score and Least Squares improved whole-cortex ICCs to 0.96 and 0.97 respectively.

CONCLUSIONS

In the presence of large intensity nonuniformities, bias field correction is necessary to achieve acceptable correspondence with known distributions of cortical myelin, but it can be detrimental in datasets with less intensity inhomogeneity. Intensity normalisation can improve test-retest reliability and inter-subject comparability. However, both bias field correction and intensity normalisation methods vary greatly in their efficacy and may affect the interpretation of results. The choice of T1w/T2w-ratio processing method must therefore be informed by both scanner and acquisition protocol as well as the given study objective. Our results highlight limitations of the T1w/T2w-ratio, but also suggest concrete ways to enhance its usefulness in future studies.

White matter plasticity in healthy older adults: The effects of aerobic exercise

White matter deterioration is associated with cognitive impairment in healthy aging and Alzheimer's disease. It is critical to identify interventions that can slow down white matter deterioration. So far, clinical trials have failed to demonstrate the benefits of aerobic exercise on the adult white matter using diffusion Magnetic Resonance Imaging. Here, we report the effects of a 6-month aerobic walking and dance interventions (clinical trial NCT01472744) on white matter integrity in healthy older adults (n = 180, 60-79 years) measured by changes in the ratio of calibrated T1- to T2-weighted images (T1w/T2w). Specifically, the aerobic walking and social dance interventions resulted in positive changes in the T1w/T2w signal in late-myelinating regions, as compared to widespread decreases in the T1w/T2w signal in the active control. Notably, in the aerobic walking group, positive change in the T1w/T2w signal correlated with improved episodic memory performance. Lastly, intervention-induced increases in cardiorespiratory fitness did not correlate with change in the T1w/T2w signal. Together, our findings suggest that white matter regions that are vulnerable to aging retain some degree of plasticity that can be induced by aerobic exercise training. In addition, we provided evidence that the T1w/T2w signal may be a useful and broadly accessible measure for studying short-term within-person plasticity and deterioration in the adult human white matter.

The Myelin Content of the Human Precentral Hand Knob Reflects Interindividual Differences in Manual Motor Control at the Physiological and Behavioral Level

The primary motor cortex hand area (M1_HAND) and adjacent dorsal premotor cortex (PMd) form the so-called motor hand knob in the precentral gyrus. M1_HAND and PMd are critical for dexterous hand use and are densely interconnected via corticocortical axons, lacking a sharp demarcating border. In 24 young right-handed volunteers, we performed multimodal mapping to delineate the relationship between structure and function in the right motor hand knob. Quantitative structural magnetic resonance imaging (MRI) at 3 tesla yielded regional R1 maps as a proxy of cortical myelin content. Participants also underwent functional MRI (fMRI). We mapped task-related activation and temporal precision, while they performed a visuomotor synchronization task requiring visually cued abduction movements with the left index or little finger. We also performed sulcus-aligned transcranial magnetic stimulation of the motor hand knob to localize the optimal site (hotspot) for evoking a motor evoked potential (MEP) in two intrinsic hand muscles. Individual motor hotspot locations varied along the rostrocaudal axis. The more rostral the motor hotspot location in the precentral crown, the longer were corticomotor MEP latencies. “Hotspot rostrality” was associated with the regional myelin content in the precentral hand knob. Cortical myelin content also correlated positively with task-related activation of the precentral crown and temporal precision during the visuomotor synchronization task. Together, our results suggest a link among cortical myelination, the spatial cortical representation, and temporal precision of finger movements. We hypothesize that the myelination of cortical axons facilitates neuronal integration in PMd and M1_HAND and, hereby, promotes the precise timing of movements.

SIGNIFICANCE STATEMENT Here we used magnetic resonance imaging and transcranial magnetic stimulation of the precentral motor hand knob to test for a link among cortical myelin content, functional corticomotor representations, and manual motor control. A higher myelin content of the precentral motor hand knob was associated with more rostral corticomotor presentations, with stronger task-related activation and a higher precision of movement timing during a visuomotor synchronization task. We propose that a high precentral myelin content enables fast and precise neuronal integration in M1 (primary motor cortex) and dorsal premotor cortex, resulting in higher temporal precision during dexterous hand use. Our results identify the degree of myelination as an important structural feature of the neocortex that is tightly linked to the function and behavior supported by the cortical area.

Cortical myelin content mediates differences in affective temperaments

BACKGROUND

Affective temperaments are regarded as subclinical forms and precursors of mental disorders. It may serve as candidates to facilitate the diagnosis and prediction of mental disorders. Cortical myelination likely characterizes the neurodevelopment and the evolution of cognitive functions and reflects brain functional demand. However, little is known about the relationship between affective temperaments and myelin plasticity. This study aims to analyze the association between the affective temperaments and cortical myelin content (CMC) in human brain.

METHODS

We measured affective temperaments using the Temperament Evaluation of Memphis, Pisa, Paris and San Diego Autoquestionnaire (TEMPS-A) on 106 healthy adults and used the ratio of T1- and T2-weighted images as the proxy for CMC. Using the unsupervised k-means clustering algorithm, we classified the cortical gray matter into heavily, intermediately, and lightly myelinated regions. The correlation between affective temperaments and CMC was calculated separately for different myelinated regions.

RESULTS

Hyperthymic temperament correlated negatively with CMC in the heavily myelinated (right postcentral gyrus and bilateral precentral gyrus) and lightly myelinated (bilateral frontal and lateral temporal) regions. Cyclothymic temperament showed a downward parabola-like correlation with CMC across the heavily, intermediately, and lightly myel0inated areas of the bilateral parietal-temporal regions.

LIMITATIONS

The analysis was constrained to cortical regions. The results were obtained from healthy subjects and we did not acquired data from patients of affective disorder, which may compromise the generalizability of the present findings.

CONCLUSION

The findings suggest that hyperthymic and cyclothymic temperaments have a CMC basis in extensive brain regions.

Overlooked Tertiary Sulci Serve as a Meso-Scale Link between Microstructural and Functional Properties of Human Lateral Prefrontal Cortex

Understanding the relationship between neuroanatomy and function in portions of cortex that perform functions largely specific to humans such as lateral prefrontal cortex (LPFC) is of major interest in systems and cognitive neuroscience. When considering neuroanatomical-functional relationships in LPFC, shallow indentations in cortex known as tertiary sulci have been largely unexplored. Here, by implementing a multimodal approach and manually defining 936 neuroanatomical structures in 72 hemispheres (in both males and females), we show that a subset of these overlooked tertiary sulci serve as a meso-scale link between microstructural (myelin content) and functional (network connectivity) properties of human LPFC in individual participants. For example, the posterior middle frontal sulcus (pmfs) is a tertiary sulcus with three components that differ in their myelin content, resting-state connectivity profiles, and engagement across meta-analyses of 83 cognitive tasks. Further, generating microstructural profiles of myelin content across cortical depths for each pmfs component and the surrounding middle frontal gyrus (MFG) shows that both gyral and sulcal components of the MFG have greater myelin content in deeper compared with superficial layers and that the myelin content in superficial layers of the gyral components is greater than sulcal components. These findings support a classic, yet largely unconsidered theory that tertiary sulci may serve as landmarks in association cortices, as well as a modern cognitive neuroscience theory proposing a functional hierarchy in LPFC. As there is a growing need for computational tools that automatically define tertiary sulci throughout cortex, we share pmfs probabilistic sulcal maps with the field.SIGNIFICANCE STATEMENT Lateral prefrontal cortex (LPFC) is critical for functions that are thought to be specific to humans compared with other mammals. However, relationships between fine-scale neuroanatomical structures largely specific to hominoid cortex and functional properties of LPFC remain elusive. Here, we show that these structures, which have been largely unexplored throughout history, surprisingly serve as markers for anatomical and functional organization in human LPFC. These findings have theoretical, methodological, developmental, and evolutionary implications for improved understanding of neuroanatomical-functional relationships not only in LPFC, but also in association cortices more broadly. Finally, these findings ignite new questions regarding how morphological features of these neglected neuroanatomical structures contribute to functions of association cortices that are critical for human-specific aspects of cognition.

Aberrant levels of cortical myelin distinguish individuals with depressive disorders from healthy controls

The association between depressive disorders and measures reflecting myelin content is underexplored, despite growing evidence of associations with white matter tract integrity. We characterized the T1w/T2w ratio using the Glasser atlas in 39 UD and 47 HC participants (ages = 19-44, 75% female). A logistic elastic net regularized regression with nested cross-validation and a subsequent linear discriminant analysis conducted on held-out samples were used to select brain regions and classify patients vs. healthy controls (HC). True-label model performance was compared against permuted-label model performance. The T1w/T2w ratio distinguished patients from HC with 68% accuracy (p < 0.001; sensitivity = 63.8%, specificity = 71.5%). Brain regions contributing to this classification performance were located in the orbitofrontal cortex, anterior cingulate, extended visual, and auditory cortices, and showed statistically significant differences in the T1w/T2w ratio for patients vs. HC. As the T1w/T2w ratio is thought to characterize cortical myelin, patterns of cortical myelin in these regions may be a biomarker distinguishing individuals with depressive disorders from HC.

Myelin quantification with MRI: A systematic review of accuracy and reproducibility

OBJECTIVES

Currently, multiple sclerosis is treated with anti-inflammatory therapies, but these treatments lack efficacy in progressive disease. New treatment strategies aim to repair myelin damage and efficacy evaluation of such new therapies would benefit from validated myelin imaging techniques. Several MRI methods for quantification of myelin density are available now. This systematic review aims to analyse the performance of these MRI methods.

METHODS

Studies comparing myelin quantification by MRI with histology, the current gold standard, or assessing reproducibility were retrieved from PubMed/MEDLINE and Embase (until December 2019). Included studies assessed both myelin histology and MRI quantitatively. Correlation or variance measurements were extracted from the studies. Non-parametric tests were used to analyse differences in study methodologies.

RESULTS

The search yielded 1348 unique articles. Twenty-two animal studies and 13 human studies correlated myelin MRI with histology. Eighteen clinical studies analysed the reproducibility. Overall bias risk was low or unclear. All MRI methods performed comparably, with a mean correlation between MRI and histology of R²=0.54 (SD=0.30) for animal studies, and R²=0.54 (SD=0.18) for human studies. Reproducibility for the MRI methods was good (ICC=0.75-0.93, R²=0.90-0.98, COV=1.3-27%), except for MTR (ICC=0.05-0.51).

CONCLUSIONS

Overall, MRI-based myelin imaging methods show a fairly good correlation with histology and a good reproducibility. However, the amount of validation data is too limited and the variability in performance between studies is too large to select the optimal MRI method for myelin quantification yet.

Effects of MP2RAGE B1+ sensitivity on inter-site T1 reproducibility and hippocampal morphometry at 7T

Most neuroanatomical studies are based on T₁-weighted MR images, whose intensity profiles are not solely determined by the tissue's longitudinal relaxation times (T₁), but also affected by varying non-T₁ contributions, hampering data reproducibility. In contrast, quantitative imaging using the MP2RAGE sequence, for example, allows direct characterization of the brain based on the tissue property of interest. Combined with 7 Tesla (7T) MRI, this offers unique opportunities to obtain robust high-resolution brain data characterized by a high reproducibility, sensitivity and specificity. However, specific MP2RAGE parameter choices - e.g., to emphasize intracortical myelin-dependent contrast variations - can substantially impact image quality and cortical analyses through remnants of B₁⁺-related intensity variations, as illustrated in our previous work. To follow up on this: we (1) validate this protocol effect using a dataset acquired with a particularly B₁⁺ insensitive set of MP2RAGE parameters combined with parallel transmission excitation; and (2) extend our analyses to evaluate the effects on hippocampal morphometry. The latter remained unexplored initially, but can provide important insights related to generalizability and reproducibility of neurodegenerative research using 7T MRI. We confirm that B₁⁺ inhomogeneities have a considerably variable effect on cortical T₁ estimates, as well as on hippocampal morphometry depending on the MP2RAGE setup. While T₁ differed substantially across datasets initially, we show the inter-site T₁ comparability improves after correcting for the spatially varying B₁⁺ field using a separately acquired Sa2RAGE B₁⁺ map. Finally, removal of B₁⁺ residuals affects hippocampal volumetry and boundary definitions, particularly near structures characterized by strong intensity changes (e.g. cerebral spinal fluid). Taken together, we show that the choice of MP2RAGE parameters can impact T₁ comparability across sites and present evidence that hippocampal segmentation results are modulated by B₁⁺ inhomogeneities. This calls for careful (1) consideration of sequence parameters when setting acquisition protocols, as well as (2) acquisition of a B₁⁺ map to correct MP2RAGE data for potential B₁⁺ variations to allow comparison across datasets.

The visual white matter connecting human area prostriata and the thalamus is retinotopically organized

The human visual system is capable of processing visual information from fovea to the far peripheral visual field. Recent fMRI studies have shown a full and detailed retinotopic map in area prostriata, located ventro-dorsally and anterior to the calcarine sulcus along the parieto-occipital sulcus with strong preference for peripheral and wide-field stimulation. Here, we report the anatomical pattern of white matter connections between area prostriata and the thalamus encompassing the lateral geniculate nucleus (LGN). To this end, we developed and utilized an automated pipeline comprising a series of Apps that run openly on the cloud computing platform brainlife.io to analyse 139 subjects of the Human Connectome Project (HCP). We observe a continuous and extended bundle of white matter fibers from which two subcomponents can be extracted: one passing ventrally parallel to the optic radiations (OR) and another passing dorsally circumventing the lateral ventricle. Interestingly, the loop travelling dorsally connects the thalamus with the central visual field representation of prostriata located anteriorly, while the other loop travelling more ventrally connects the LGN with the more peripheral visual field representation located posteriorly. We then analyse an additional cohort of 10 HCP subjects using a manual plane extraction method outside brainlife.io to study the relationship between the two extracted white matter subcomponents and eccentricity, myelin and cortical thickness gradients within prostriata. Our results are consistent with a retinotopic segregation recently demonstrated in the OR, connecting the LGN and V1 in humans and reveal for the first time a retinotopic segregation regarding the trajectory of a fiber bundle between the thalamus and an associative visual area.

B0-field dependence of MRI T1 relaxation in human brain

Tissue longitudinal relaxation characterized by recovery time T1 or rate R1 is a fundamental MRI contrast mechanism that is increasingly being used to study the brain's myelination patterns in both health and disease. Nevertheless, the quantitative relationship between T1 and myelination, and its dependence on B0 field strength, is still not well known. It has been theorized that in much of brain tissue, T1 field-dependence is driven by that of macromolecular protons (MP) through a mechanism called magnetization transfer (MT). Despite the explanatory power of this theory and substantial support from in-vitro experiments at low fields (<3 ​T), in-vivo evidence across clinically relevant field strengths is lacking. In this study, T1-weighted MRI was acquired in a group of eight healthy volunteers at four clinically relevant field strengths (0.55, 1.5, 3 and 7 ​T) using the same pulse sequence at a single site, and jointly analyzed based on the two-pool model of MT. MP fraction and free-water pool T1 were obtained in several brain structures at 3 and 7 ​T, which allowed distinguishing between contributions from macromolecular content and iron to tissue T1. Based on this, the T1 of MP in white matter, indirectly determined by assuming a field independent T1 of free water, was shown to increase approximately linearly with B0. This study advances our understanding of the T1 contrast mechanism and its relation to brain myelin content across the wide range of currently available MRI strengths, and it has the potential to inform design of T1 mapping methods for improved reproducibility in the human brain.