Detection of pathological contrast enhancement with synthetic brain imaging from quantitative multiparametric MRI

BACKGROUND AND PURPOSE

We aimed to test whether synthetic T1-weighted imaging derived from a post-contrast Quantitative Transient-state Imaging (QTI) acquisition enabled revealing pathological contrast enhancement in intracranial lesions.

METHODS

The analysis included 141 patients who underwent a 3 Tesla-MRI brain exam with intravenous contrast media administration, with the post-contrast acquisition protocol comprising a three-dimensional fast spoiled gradient echo (FSPGR) sequence and a QTI acquisition. Synthetic T1-weighted images were generated from QTI-derived quantitative maps of relaxation times and proton density. Two neuroradiologists assessed synthetic and conventional post-contrast T1-weighted images for the presence and pattern of pathological contrast enhancement in intracranial lesions. Enhancement volumes were quantitatively compared.

RESULTS

Using conventional imaging as a reference, synthetic T1-weighted imaging was 93% sensitive in revealing the presence of contrast enhancing lesions. The agreement for the presence/absence of contrast enhancement was almost perfect both between readers (k = 1 for both conventional and synthetic imaging) and between sequences (k = 0.98 for both readers). In 91% of lesions, synthetic T1-weighted imaging showed the same pattern of contrast enhancement visible in conventional imaging. Differences in enhancement pattern in the remaining lesions can be due to the lower spatial resolution and the longer acquisition delay from contrast media administration of QTI compared to FSPGR. Overall, enhancement volumes appeared larger in synthetic imaging.

CONCLUSIONS

QTI-derived post-contrast synthetic T1-weighted imaging captures pathological contrast enhancement in most intracranial enhancing lesions. Further comparative studies employing quantitative imaging with higher spatial resolution is needed to support our data and explore possible future applications in clinical trials.

Magnetic resonance fingerprinting-based myelin water fraction mapping for the assessment of white matter maturation and integrity in typical development and leukodystrophies

A quantitative biomarker for myelination, such as myelin water fraction (MWF), would boost the understanding of normative and pathological neurodevelopment, improving patients' diagnosis and follow-up. We quantified the fraction of a rapidly relaxing pool identified as MW using multicomponent three-dimensional (3D) magnetic resonance fingerprinting (MRF) to evaluate white matter (WM) maturation in typically developing (TD) children and alterations in leukodystrophies (LDs). We acquired DTI and 3D MRF-based R1, R2 and MWF data of 15 TD children and 17 LD patients (9 months-12.5 years old) at 1.5 T. We computed normative maturation curves in corpus callosum and corona radiata and performed WM tract profile analysis, comparing MWF with R1, R2 and fractional anisotropy (FA). Normative maturation curves demonstrated a steep increase for all tissue parameters in the first 3 years of age, followed by slower growth for MWF while R1, R2R2 and FA reached a plateau. Unlike FA, MWF values were similar for regions of interest (ROIs) with different degrees of axonal packing, suggesting independence from fiber bundle macro-organization and higher myelin specificity. Tract profile analysis indicated a specific spatial pattern of myelination in the major fiber bundles, consistent across subjects. LD were better distinguished from TD by MWF rather than FA, showing reduced MWF with respect to age-matched controls in both ROI-based and tract analysis. In conclusion, MRF-based MWF provides myelin-specific WM maturation curves and is sensitive to alteration due to LDs, suggesting its potential as a biomarker for WM disorders. As MRF allows fast simultaneous acquisition of relaxometry and MWF, it can represent a valuable diagnostic tool to study and follow up developmental WM disorders in children.

Case report: Exploring chemoradiotherapy-induced leukoencephalopathy with 7T imaging and quantitative susceptibility mapping

Chemotherapy and radiotherapy are widely used in the treatment of central nervous system tumors and acute lymphocytic leukemia even in the pediatric population. However, such treatments run the risk of a broad spectrum of cognitive and neurological deficits. Even though the correlation with cognitive decline is still not clear, neuroradiological defects linked to white matter injury and vasculopathies may be identified. Thanks to the use of 7T MRI it is possible to better define the vascular pattern of the brain lesions with the added advantage of identifying their characteristics and anatomical localization, which, however, are not evident with a conventional brain scan. Moreover, the use of Quantitative Susceptibility Mapping (QSM) makes it possible to discriminate between calcium deposits on vessels (chemo-radiation-induced) and hemoglobin deposition in radio-induced cavernomas, speculating, as a result, about the pathophysiology of iatrogenic brain damage. We describe the case of a 9 year-old boy with a T-type acute lymphoid leukemia who had previously been treated with polychemotherapy and high-dose RT. To better define the child's neuroradiological pattern, 7T MRI and QSM were performed in addition to conventional imaging examinations. Our case report suggests the potential usefulness of a QSM study to distinguish radio-induced vascular malformations from mineralizing microangiopathy.

Fast high-resolution electric properties tomography using three-dimensional quantitative transient-state imaging-based water fraction estimation

In this study, we aimed to develop a fast and robust high-resolution technique for clinically feasible electrical properties tomography based on water content maps (wEPT) using Quantitative Transient-state Imaging (QTI), a multiparametric transient state-based method that is similar to MR fingerprinting. Compared with the original wEPT implementation based on standard spin-echo acquisition, QTI provides robust electrical properties quantification towards B1 + inhomogeneities and full quantitative relaxometry data. To validate the proposed approach, 3D QTI data of 12 healthy volunteers were acquired on a 1.5 T scanner. QTI-provided T1 maps were used to compute water content maps of the tissues using an empirical relationship based on literature ex-vivo measurements. Assuming that electrical properties are modulated mainly by tissue water content, the water content maps were used to derive electrical conductivity and relative permittivity maps. The proposed technique was compared with a conventional phase-only Helmholtz EPT (HH-EPT) acquisition both within whole white matter, gray matter, and cerebrospinal fluid masks, and within different white and gray matter subregions. In addition, QTI-based wEPT was retrospectively applied to four multiple sclerosis adolescent and adult patients, compared with conventional contrast-weighted imaging in terms of lesion delineation, and quantitatively assessed by measuring the variation of electrical properties in lesions. Results obtained with the proposed approach agreed well with theoretical predictions and previous in vivo findings in both white and gray matter. The reconstructed maps showed greater anatomical detail and lower variability compared with standard phase-only HH-EPT. The technique can potentially improve delineation of pathology when compared with conventional contrast-weighted imaging and was able to detect significant variations in lesions with respect to normal-appearing tissues. In conclusion, QTI can reliably measure conductivity and relative permittivity of brain tissues within a short scan time, opening the way to the study of electric properties in clinical settings.

Brainstem anatomy with 7-T MRI: in vivo assessment and ex vivo comparison

BACKGROUND

The brainstem contains grey matter nuclei and white matter tracts to be identified in clinical practice. The small size and the low contrast among them make their in vivo visualisation challenging using conventional magnetic resonance imaging (MRI) sequences at high magnetic field strengths. Combining higher spatial resolution, signal- and contrast-to-noise ratio and sensitivity to magnetic susceptibility (χ), susceptibility-weighted 7-T imaging could improve the assessment of brainstem anatomy.

METHODS

We acquired high-resolution 7-T MRI of the brainstem in a 46-year-old female healthy volunteer (using a three-dimensional multi-echo gradient-recalled-echo sequence; spatial resolution 0.3 × 0.3 × 1.2 mm3) and in a brainstem sample from a 48-year-old female body donor that was sectioned and stained. Images were visually assessed; nuclei and tracts were labelled and named according to the official nomenclature.

RESULTS

This in vivo imaging revealed structures usually evaluated through light microscopy, such as the accessory olivary nuclei, oculomotor nucleus and the medial longitudinal fasciculus. Some fibre tracts, such as the medial lemniscus, were visible for most of their course. Overall, in in vivo acquisitions, χ and frequency maps performed better than T2*-weighted imaging and allowed for the evaluation of a greater number of anatomical structures. All the structures identified in vivo were confirmed by the ex vivo imaging and histology.

CONCLUSIONS

The use of multi-echo GRE sequences at 7 T allowed the visualisation of brainstem structures that are not visible in detail at conventional magnetic field and opens new perspectives in the diagnostic and therapeutical approach to brain disorders.

RELEVANCE STATEMENT

In vivo MR imaging at UHF provides detailed anatomy of CNS substructures comparable to that obtained with histology. Anatomical details are fundamentals for diagnostic purposes but also to plan a direct targeting for a minimally invasive brain stimulation or ablation.

KEY POINTS

• The in vivo brainstem anatomy was explored with ultrahigh field MRI (7 T). • In vivo T2*-weighted magnitude, χ, and frequency images revealed many brainstem structures. • Ex vivo imaging and histology confirmed all the structures identified in vivo. • χ and frequency imaging revealed more brainstem structures than magnitude imaging.

Quantitative T1 mapping detects blood-brain barrier breakdown in apparently non-enhancing multiple sclerosis lesions

OBJECTIVES

The disruption of the blood-brain barrier (BBB) is a key and early feature in the pathogenesis of demyelinating multiple sclerosis (MS) lesions and has been neuropathologically demonstrated in both active and chronic plaques. The local overt BBB disruption in acute demyelinating lesions is captured as signal hyperintensity in post-contrast T1-weighted images because of the contrast-related shortening of the T1 relaxation time. On the contrary, the subtle BBB disruption in chronic lesions is not visible at conventional radiological evaluation but it might be of clinical relevance. Indeed, persistent, subtle BBB leakage might be linked to low-grade inflammation and plaque evolution. Here we hypothesised that 3D Quantitative Transient-state Imaging (QTI) was able to reveal and measure T1 shortening (ΔT1) reflecting small amounts of contrast media leakage in apparently non-enhancing lesions (ANELs).

MATERIALS AND METHODS

Thirty-four patients with relapsing remitting MS were included in the study. All patients underwent a 3 T MRI exam of the brain including conventional sequences and QTI acquisitions (1.1 mm isotropic voxel) performed both before and after contrast media administration. For each patient, a ΔT1 map was obtained via voxel-wise subtraction of pre- and post- contrast QTI-derived T1 maps. ΔT1 values measured in ANELs were compared with those recorded in enhancing lesions and in the normal appearing white matter. A reference distribution of ΔT1 in the white matter was obtained from datasets acquired in 10 non-MS patients with unrevealing MR imaging.

RESULTS

Mean ΔT1 in ANELs (57.45 ± 48.27 ms) was significantly lower than in enhancing lesions (297.71 ± 177.52 ms; p < 0. 0001) and higher than in the normal appearing white matter (36.57 ± 10.53 ms; p < 0.005). Fifty-two percent of ANELs exhibited ΔT1 higher than those observed in the white matter of non-MS patients.

CONCLUSIONS

QTI-derived quantitative ΔT1 mapping enabled to measure contrast-related T1 shortening in ANELs. ANELs exhibiting ΔT1 values that deviate from the reference distribution in non-MS patients may indicate persistent, subtle, BBB disruption. Access to this information may be proved useful to better characterise pathology and objectively monitor disease activity and response to therapy.

Generating Synthetic Radiological Images with PySynthMRI: An Open-Source Cross-Platform Tool

Synthetic MR Imaging allows for the reconstruction of different image contrasts from a single acquisition, reducing scan times. Commercial products that implement synthetic MRI are used in research. They rely on vendor-specific acquisitions and do not include the possibility of using custom multiparametric imaging techniques. We introduce PySynthMRI, an open-source tool with a user-friendly interface that uses a set of input images to generate synthetic images with diverse radiological contrasts by varying representative parameters of the desired target sequence, including the echo time, repetition time and inversion time(s). PySynthMRI is written in Python 3.6, and it can be executed under Linux, Windows, or MacOS as a python script or an executable. The tool is free and open source and is developed while taking into consideration the possibility of software customization by the end user. PySynthMRI generates synthetic images by calculating the pixelwise signal intensity as a function of a set of input images (e.g., T1 and T2 maps) and simulated scanner parameters chosen by the user via a graphical interface. The distribution provides a set of default synthetic contrasts, including T1w gradient echo, T2w spin echo, FLAIR and Double Inversion Recovery. The synthetic images can be exported in DICOM or NiFTI format. PySynthMRI allows for the fast synthetization of differently weighted MR images based on quantitative maps. Specialists can use the provided signal models to retrospectively generate contrasts and add custom ones. The modular architecture of the tool can be exploited to add new features without impacting the codebase.

Normative values of the topological metrics of the structural connectome: A multi-site reproducibility study across the Italian Neuroscience network

PURPOSE

The use of topological metrics to derive quantitative descriptors from structural connectomes is receiving increasing attention but deserves specific studies to investigate their reproducibility and variability in the clinical context. This work exploits the harmonization of diffusion-weighted acquisition for neuroimaging data performed by the Italian Neuroscience and Neurorehabilitation Network initiative to obtain normative values of topological metrics and to investigate their reproducibility and variability across centers.

METHODS

Different topological metrics, at global and local level, were calculated on multishell diffusion-weighted data acquired at high-field (e.g. 3 T) Magnetic Resonance Imaging scanners in 13 different centers, following the harmonization of the acquisition protocol, on young and healthy adults. A "traveling brains" dataset acquired on a subgroup of subjects at 3 different centers was also analyzed as reference data. All data were processed following a common processing pipeline that includes data pre-processing, tractography, generation of structural connectomes and calculation of graph-based metrics. The results were evaluated both with statistical analysis of variability and consistency among sites with the traveling brains range. In addition, inter-site reproducibility was assessed in terms of intra-class correlation variability.

RESULTS

The results show an inter-center and inter-subject variability of <10%, except for "clustering coefficient" (variability of 30%). Statistical analysis identifies significant differences among sites, as expected given the wide range of scanners' hardware.

CONCLUSIONS

The results show low variability of connectivity topological metrics across sites running a harmonised protocol.

Quality assessment, variability and reproducibility of anatomical measurements derived from T1-weighted brain imaging: The RIN-Neuroimaging Network case study

Initiatives for the collection of harmonized MRI datasets are growing continuously, opening questions on the reliability of results obtained in multi-site contexts. Here we present the assessment of the brain anatomical variability of MRI-derived measurements obtained from T1-weighted images, acquired according to the Standard Operating Procedures, promoted by the RIN-Neuroimaging Network. A multicentric dataset composed of 77 brain T1w acquisitions of young healthy volunteers (mean age = 29.7 ± 5.0 years), collected in 15 sites with MRI scanners of three different vendors, was considered. Parallelly, a dataset of 7 "traveling" subjects, each undergoing three acquisitions with scanners from different vendors, was also used. Intra-site, intra-vendor, and inter-site variabilities were evaluated in terms of the percentage standard deviation of volumetric and cortical thickness measures. Image quality metrics such as contrast-to-noise and signal-to-noise ratio in gray and white matter were also assessed for all sites and vendors. The results showed a measured global variability that ranges from 11% to 19% for subcortical volumes and from 3% to 10% for cortical thicknesses. Univariate distributions of the normalized volumes of subcortical regions, as well as the distributions of the thickness of cortical parcels appeared to be significantly different among sites in 8 subcortical (out of 17) and 21 cortical (out of 68) regions of i nterest in the multicentric study. The Bland-Altman analysis on "traveling" brain measurements did not detect systematic scanner biases even though a multivariate classification approach was able to classify the scanner vendor from brain measures with an accuracy of 0.60 ± 0.14 (chance level 0.33).

Simulated and experimental phantom data for multi-center quality assurance of quantitative susceptibility maps at 3 T, 7 T and 9.4 T

PURPOSE

To develop methods for quality assurance of quantitative susceptibility mapping (QSM) using MRI at different magnetic field strengths, and scanners, using different MR-sequence protocols, and post-processing pipelines.

METHODS

We built a custom phantom based on iron in two forms: homogeneous susceptibility ('free iron') and with fine-scaled variations in susceptibility ('clustered iron') at different iron concentrations. The phantom was measured at 3.0 T (two scanners), 7.0 T and 9.4 T using multi-echo, gradient echo acquisition sequences. A digital phantom analogue to the iron-phantom, tailored to obtain similar results as in experimentation was developed, with similar geometry and susceptibility values. Morphology enabled dipole inversion was applied to the phase images to obtain QSM for experimental and simulated data using the MEDI + 0 approach for background regularization.

RESULTS

Across all scanners, QSM-values showed a linear increase with iron concentrations. The QSM-relaxivity was 0.231 ± 0.047 ppm/mM for free and 0.054 ± 0.013 ppm/mM for clustered iron, with adjusted determination coefficients (DoC) ≥ 0.87. Similarly, the simulations yielded linear increases (DoC ≥ 0.99). In both the experimental and digital phantoms, the estimated molar susceptibility was lower with clustered iron, because clustering led to highly localized field effects.

CONCLUSION

Our iron phantom can be used to evaluate the capability of QSM to detect local variations in susceptibility across different field strengths, when using different MR-sequence protocols. The devised simulation method captures the effect of iron clustering in QSM as seen experimentally and could be used in the future to optimize QSM processing pipelines and achieve higher accuracy for local field effects, as also seen in Alzheimer's beta-amyloid plaques.

DEVELOPMENT OF BOLD RESPONSE TO MOTION IN HUMAN INFANTS

Behavioral studies suggest that motion perception is rudimentary at birth and matures steadily over the first few years. We demonstrated previously that the major cortical associative areas serving motion processing, like MT+, V6 and PVICS in adults, show selective responses to coherent flow in 8-week-old infants. Here we study the BOLD response to the same motion stimuli in 5-week-old infants (four females and four males) and compare the maturation between these two ages. The results show that MT+ and PVIC areas show a similar motion response at 5 and 8 weeks, while response in the V6 area shows a reduced BOLD response to motion at 5 weeks and Cuneus associative areas are not identifiable at this young age. In infants and in adults, V1 does not show a selectivity for coherent motion, but shows a very fast development between 5 and 8 weeks of age in response to appearance of motion stimuli. Resting-state correlations demonstrate adult-like functional connectivity between the motion-selective associative areas, but not between primary cortex and temporo-occipital and posterior-insular cortices. The results are consistent with a differential developmental trajectory of motion area respect to other occipital regions, probably reflecting also a different development trajectory of the central and peripheral visual field.SIGNIFICANCE STATEMENT:How the cortical visual areas attain the specialization that we observed in human adults in the first few months of life is unknown. However, this knowledge is crucial to understanding the consequence of peri-natal brain damage and its outcome. Here we show that motion selective areas are already functioning well in 5-week-old infants with greater responses for detecting coherent motion over random motion, suggesting that very little experience is needed to attain a motion selectivity.

Virtual Ontogeny of Cortical Growth Preceding Mental Illness

BACKGROUND

Morphology of the human cerebral cortex differs across psychiatric disorders, with neurobiology and developmental origins mostly undetermined. Deviations in the tangential growth of the cerebral cortex during pre/perinatal periods may be reflected in individual variations in cortical surface area later in life.

METHODS

Interregional profiles of group differences in surface area between cases and controls were generated using T1-weighted magnetic resonance imaging from 27,359 individuals including those with attention-deficit/hyperactivity disorder, autism spectrum disorder, bipolar disorder, major depressive disorder, schizophrenia, and high general psychopathology (through the Child Behavior Checklist). Similarity of interregional profiles of group differences in surface area and prenatal cell-specific gene expression was assessed.

RESULTS

Across the 11 cortical regions, group differences in cortical area for attention-deficit/hyperactivity disorder, schizophrenia, and Child Behavior Checklist were dominant in multimodal association cortices. The same interregional profiles were also associated with interregional profiles of (prenatal) gene expression specific to proliferative cells, namely radial glia and intermediate progenitor cells (greater expression, larger difference), as well as differentiated cells, namely excitatory neurons and endothelial and mural cells (greater expression, smaller difference). Finally, these cell types were implicated in known pre/perinatal risk factors for psychosis. Genes coexpressed with radial glia were enriched with genes implicated in congenital abnormalities, birth weight, hypoxia, and starvation. Genes coexpressed with endothelial and mural genes were enriched with genes associated with maternal hypertension and preterm birth.

CONCLUSIONS

Our findings support a neurodevelopmental model of vulnerability to mental illness whereby prenatal risk factors acting through cell-specific processes lead to deviations from typical brain development during pregnancy.

Evaluation of iron overload in nigrosome 1 via quantitative susceptibility mapping as a progression biomarker in prodromal stages of synucleinopathies

Idiopathic rapid eye movement (REM) sleep behavior disorder (iRBD) is a prodromal stage of α-synucleinopathies, such as Parkinson's disease (PD), which are characterized by the loss of dopaminergic neurons in substantia nigra, associated with abnormal iron load. The assessment of presymptomatic biomarkers predicting the onset of neurodegenerative disorders is critical for monitoring early signs, screening patients for neuroprotective clinical trials and understanding the causal relationship between iron accumulation processes and disease development. Here, we used Quantitative Susceptibility Mapping (QSM) and 7T MRI to quantify iron deposition in Nigrosome 1 (N1) in early PD (ePD) patients, iRBD patients and healthy controls and investigated group differences and correlation with disease progression. We evaluated the radiological appearance of N1 and analyzed its iron content in 35 ePD, 30 iRBD patients and 14 healthy controls via T2*-weighted sequences and susceptibility (χ) maps. N1 regions of interest (ROIs) were manually drawn on control subjects and warped onto a study-specific template to obtain probabilistic N1 ROIs. For each subject the N1 with the highest mean χ was considered for statistical analysis. The appearance of N1 was rated pathological in 45% of iRBD patients. ePD patients showed increased N1 χ compared to iRBD patients and HC but no correlation with disease duration, indicating that iron load remains stable during the early stages of disease progression. Although no difference was reported in iron content between iRBD and HC, N1 χ in the iRBD group increases as the disease evolves. QSM can reveal temporal changes in N1 iron content and its quantification may represent a valuable presymptomatic biomarker to assess neurodegeneration in the prodromal stages of PD.

Diagnostic accuracy of quantitative susceptibility mapping in multiple system atrophy: The impact of echo time and the potential of histogram analysis

•

We performed histogram analysis on χ maps at different TEs on MSA patients and HC.

•

We found altered χ distribution in Pu, SN, GP, CN for MSAp and in SN, DN for MSAc.

•

QSM diagnostic accuracy is TE-dependent and is enhanced at short TEs.

•

Short TEs capture rapidly-decaying contributions of high χ sources.

•

Histogram features detect χ spatial heterogeneities improving diagnostic accuracy.

The non-invasive quantification of iron stores via Quantitative Susceptibility Mapping (QSM) could play an important role in the diagnosis and the differential diagnosis of atypical Parkinsonisms. However, the susceptibility (χ) values measured via QSM depend on echo time (TE). This effect relates to the microstructural organization within the voxel, whose composition can be altered by the disease. Moreover, pathological iron deposition in a brain area may not be spatially uniform, and conventional Region of Interest (ROI)-based analysis may fail in detecting alterations. Therefore, in this work we evaluated the impact of echo time on the diagnostic accuracy of QSM on a population of patients with Multiple System Atrophy (MSA) of either Parkinsonian (MSAp) or cerebellar (MSAc) phenotypes. In addition, we tested the potential of histogram analysis to improve QSM classification accuracy.

We enrolled 32 patients (19 MSAp and 13 MSAc) and 16 healthy controls, who underwent a 7T MRI session including a gradient-recalled multi-echo sequence for χ mapping. Nine histogram features were extracted from the χ maps computed for each TE in atlas-based ROIs covering deep brain nuclei, and compared among groups.

Alterations of susceptibility distribution were found in the Putamen, Substantia Nigra, Globus Pallidus and Caudate Nucleus for MSAp and in the Substantia Nigra and Dentate Nucleus for MSAc. Increased iron deposition was observed in a larger number of ROIs for the two shortest TEs and the standard deviation, the 75^th and the 90^th percentile were the most informative features yielding excellent diagnostic accuracy with area under the ROC curve > 0.9.

In conclusion, short TEs may enhance QSM diagnostic performances, as they can capture variations in rapidly-decaying contributions of high χ sources. The analysis of histogram features allowed to reveal fine heterogeneities in the spatial distribution of susceptibility alteration, otherwise undetected by a simple evaluation of ROI χ mean values.

Quantitative MRI Harmonization to Maximize Clinical Impact: The RIN-Neuroimaging Network

Neuroimaging studies often lack reproducibility, one of the cardinal features of the scientific method. Multisite collaboration initiatives increase sample size and limit methodological flexibility, therefore providing the foundation for increased statistical power and generalizable results. However, multisite collaborative initiatives are inherently limited by hardware, software, and pulse and sequence design heterogeneities of both clinical and preclinical MRI scanners and the lack of benchmark for acquisition protocols, data analysis, and data sharing. We present the overarching vision that yielded to the constitution of RIN-Neuroimaging Network, a national consortium dedicated to identifying disease and subject-specific in-vivo neuroimaging biomarkers of diverse neurological and neuropsychiatric conditions. This ambitious goal needs efforts toward increasing the diagnostic and prognostic power of advanced MRI data. To this aim, 23 Italian Scientific Institutes of Hospitalization and Care (IRCCS), with technological and clinical specialization in the neurological and neuroimaging field, have gathered together. Each IRCCS is equipped with high- or ultra-high field MRI scanners (i.e., ≥3T) for clinical or preclinical research or has established expertise in MRI data analysis and infrastructure. The actions of this Network were defined across several work packages (WP). A clinical work package (WP1) defined the guidelines for a minimum standard clinical qualitative MRI assessment for the main neurological diseases. Two neuroimaging technical work packages (WP2 and WP3, for clinical and preclinical scanners) established Standard Operative Procedures for quality controls on phantoms as well as advanced harmonized quantitative MRI protocols for studying the brain of healthy human participants and wild type mice. Under FAIR principles, a web-based e-infrastructure to store and share data across sites was also implemented (WP4). Finally, the RIN translated all these efforts into a large-scale multimodal data collection in patients and animal models with dementia (i.e., case study). The RIN-Neuroimaging Network can maximize the impact of public investments in research and clinical practice acquiring data across institutes and pathologies with high-quality and highly-consistent acquisition protocols, optimizing the analysis pipeline and data sharing procedures.

Short-term plasticity in the human visual thalamus

While there is evidence that the visual cortex retains a potential for plasticity in adulthood, less is known about the subcortical stages of visual processing. Here we asked whether short-term ocular dominance plasticity affects the human visual thalamus. We addressed this question in normally sighted adult humans, using ultra-high field (7T) magnetic resonance imaging combined with the paradigm of short-term monocular deprivation. With this approach, we previously demonstrated transient shifts of perceptual eye dominance and ocular dominance in visual cortex (Binda et al., 2018). Here we report evidence for short-term plasticity in the ventral division of the pulvinar (vPulv), where the deprived eye representation was enhanced over the non-deprived eye. This ventral-pulvinar plasticity was similar as previously seen in visual cortex and it was correlated with the ocular dominance shift measured behaviorally. In contrast, there was no effect of monocular deprivation in two adjacent thalamic regions: dorsal pulvinar (dPulv), and Lateral Geniculate Nucleus (LGN). We conclude that the visual thalamus retains potential for short-term plasticity in adulthood; the plasticity effect differs across thalamic subregions, possibly reflecting differences in their cortico-fugal connectivity.

Subtly altered topological asymmetry of brain structural covariance networks in autism spectrum disorder across 43 datasets from the ENIGMA consortium

Small average differences in the left-right asymmetry of cerebral cortical thickness have been reported in individuals with autism spectrum disorder (ASD) compared to typically developing controls, affecting widespread cortical regions. The possible impacts of these regional alterations in terms of structural network effects have not previously been characterized. Inter-regional morphological covariance analysis can capture network connectivity between different cortical areas at the macroscale level. Here, we used cortical thickness data from 1455 individuals with ASD and 1560 controls, across 43 independent datasets of the ENIGMA consortium's ASD Working Group, to assess hemispheric asymmetries of intra-individual structural covariance networks, using graph theory-based topological metrics. Compared with typical features of small-world architecture in controls, the ASD sample showed significantly altered average asymmetry of networks involving the fusiform, rostral middle frontal, and medial orbitofrontal cortex, involving higher randomization of the corresponding right-hemispheric networks in ASD. A network involving the superior frontal cortex showed decreased right-hemisphere randomization. Based on comparisons with meta-analyzed functional neuroimaging data, the altered connectivity asymmetry particularly affected networks that subserve executive functions, language-related and sensorimotor processes. These findings provide a network-level characterization of altered left-right brain asymmetry in ASD, based on a large combined sample. Altered asymmetrical brain development in ASD may be partly propagated among spatially distant regions through structural connectivity.

Simultaneous relaxometry and morphometry of human brain structures with 3D magnetic resonance fingerprinting: a multicenter, multiplatform, multifield-strength study

Relaxation times and morphological information are fundamental magnetic resonance imaging-derived metrics of the human brain that reflect the status of the underlying tissue. Magnetic resonance fingerprinting (MRF) enables simultaneous acquisition of T1 and T2 maps inherently aligned to the anatomy, allowing whole-brain relaxometry and morphometry in a single scan. In this study, we revealed the feasibility of 3D MRF for simultaneous brain structure-wise morphometry and relaxometry. Comprehensive test-retest scan analyses using five 1.5-T and three 3.0-T systems from a single vendor including different scanner types across 3 institutions demonstrated that 3D MRF-derived morphological information and relaxation times are highly repeatable at both 1.5 T and 3.0 T. Regional cortical thickness and subcortical volume values showed high agreement and low bias across different field strengths. The ability to acquire a set of regional T1, T2, thickness, and volume measurements of neuroanatomical structures with high repeatability and reproducibility facilitates the ability of longitudinal multicenter imaging studies to quantitatively monitor changes associated with underlying pathologies, disease progression, and treatments.

Complementing canonical fMRI with functional Quantitative Susceptibility Mapping (fQSM) in modern neuroimaging research

Functional Quantitative Susceptibility Mapping (fQSM) allows for the quantitative measurement of time-varying magnetic susceptibility across cortical and subcortical brain structures with a potentially higher spatial specificity than conventional fMRI. While the usefulness of fQSM with General Linear Model and "On/Off" paradigms has been assessed, little is known about the potential applications and limitations of this technique in more sophisticated experimental paradigms and analyses, such as those currently used in modern neuroimaging. To thoroughly characterize fQSM activations, here we used 7T MRI, tonotopic mapping, as well as univariate (i.e., GLM and population Receptive Field) and multivariate (Representational Similarity Analysis; RSA) analyses. Although fQSM detected less tone-responsive voxels than fMRI, they were more consistently localized in gray matter. Also, the majority of active gray matter voxels exhibited negative fQSM response, signaling the expected oxyhemoglobin increase, whereas positive fQSM activations were mainly in white matter. Though fMRI- and fQSM-based tonotopic maps were overall comparable, the representation of frequency tunings in tone-sensitive regions was significantly more balanced for fQSM. Lastly, RSA revealed that frequency information from the auditory cortex could be successfully retrieved by using either methods. Overall, fQSM produces complementary results to conventional fMRI, as it captures small-scale variations in the activation pattern which inform multivariate measures. Although positive fQSM responses deserve further investigation, they do not impair the interpretation of contrasts of interest. The quantitative nature of fQSM, its spatial specificity and the possibility to simultaneously acquire canonical fMRI support the use of this technique for longitudinal and multicentric studies and pre-surgical mapping.

Predictive visuo-motor communication through neural oscillations

The mechanisms coordinating action and perception over time are poorly understood. The sensory cortex needs to prepare for upcoming changes contingent on action, and this requires temporally precise communication that takes into account the variable delays between sensory and motor processing. Several theorists^1,2 have proposed synchronization of the endogenous oscillatory activity observed in most regions of the brain³ as the basis for an efficient and flexible communication protocol between distal brain areas,^2,4 a concept known as "communication through coherence." Synchronization of endogenous oscillations^5,6 occurs after a salient sensory stimulus, such as a flash or a sound,^7-11 and after a voluntary action,^12-18 and this directly impacts perception, causing performance to oscillate rhythmically over time. Here we introduce a novel fMRI paradigm to probe the neural sources of oscillations, based on the concept of perturbative signals, which overcomes the low temporal resolution of BOLD signals. The assumption is that a synchronized endogenous rhythm will modulate cortical excitability rhythmically, which should be reflected in the BOLD responses to brief stimuli presented at different phases of the oscillation cycle. We record rhythmic oscillations of V1 BOLD synchronized by a simple voluntary action, in phase with behaviorally measured oscillations in visual sensitivity in the theta range. The functional connectivity between V1 and M1 also oscillates at the same rhythm. By demonstrating oscillatory temporal coupling between primary motor and sensory cortices, our results strongly implicate communication through coherence to achieve precise coordination and to encode sensory-motor timing.

Neural substrates of neuropsychological profiles in dystrophynopathies: A pilot study of diffusion tractography imaging

Introduction

Cognitive difficulties and neuropsychological alterations in Duchenne and Becker muscular dystrophy (DMD, BMD) boys are not yet sufficiently explored, although this topic could have a relevant impact, finding novel biomarkers of disease both at genetics and neuroimaging point of view. The current study aims to: 1) analyze the neuropsychological profile of a group of DMD and BMD boys without cognitive impairment with an assessment of their executive functions; 2) explore the structural connectivity in DMD, BMD, and age-matched controls focusing on cortico-subcortical tracts that connect frontal cortex, basal ganglia, and cerebellum via the thalamus; 3) explore possible correlations between altered structural connectivity and clinical neuropsychological measures.

Materials and methods

This pilot study included 15 boys (5 DMD subjects, 5 BMD subjects, and 5 age-matched typically developing, TD). They were assessed using a neuropsychological assessment protocol including cognitive and executive functioning assessment and performed a 1.5T MRI brain exam including advance Diffusion Weighted Imaging (DWI) method for tractography. Structural connectivity measurements were extracted along three specific tracts: Cortico-Ponto-Cerebellar Tract (CPCT), Cerebellar-Thalamic Tract (CTT), and Superior Longitudinal Fasciculus (SLF). Cortical-Spinal Tract (CST) was selected for reference, as control tract.

Results

Regarding intellectual functioning, a major impairment in executive functions compared to the general intellectual functioning was observed both for DMD (mean score = 86.20; SD = 11.54) and for BMD children (mean score = 88; SD = 3.67). Mean FA resulted tendentially always lower in DMD compared to both BMD and TD groups for all the examined tracts. The differences in FA were statistically significant for the right CTT (DMD vs BMD, p = 0.002, and DMD vs TD, p = 0.0015) and the right CPCT (DMD vs TD, p = 0.008). Concerning DMD, significant correlations emerged between FA-R-CTT and intellectual quotients (FIQ, p = 0.044; ρs = 0.821), and executive functions (Denomination Total, p = 0.044, ρs = 0.821; Inhibition Total, p = 0.019, ρs = 0.900). BMD showed a significant correlation between FA-R-CPCT and working memory index (p = 0.007; ρs = 0.949).

Discussion and conclusion

In this pilot study, despite the limitation of sample size, the findings support the hypothesis of the involvement of a cerebellar-thalamo-cortical loop for the neuropsychological profile of DMD, as the CTT and the CPCT are involved in the network and the related brain structures are known to be implied in executive functions. Our results suggest that altered WM connectivity and reduced fibre organization in cerebellar tracts, probably due to the lack of dystrophin in the brain, may render less efficient some neuropsychological functions in children affected by dystrophinopathies. The wider multicentric study could help to better establish the role of cerebellar connectivity in neuropsychological profile for dystrophinopathies, identifying possible novel diagnostic and prognostic biomarkers.

Quantitative imaging metrics derived from magnetic resonance fingerprinting using ISMRM/NIST MRI system phantom: An international multicenter repeatability and reproducibility study

PURPOSE

To compare the bias and inherent reliability of the quantitative (T₁ and T₂ ) imaging metrics generated from the magnetic resonance fingerprinting (MRF) technique using the ISMRM/NIST system phantom in an international multicenter setting.

METHOD

ISMRM/NIST MRI system phantom provides standard reference T₁ and T₂ relaxation values (vendor-provided) for each of the 14 vials in T₁ and T₂ arrays. MRF-SSFP scans repeated over 30 days on GE 1.5 and 3.0 T scanners at three collaborative centers. MRF estimated T_1, and T₂ values averaged over 30 days were compared with the phantom vendor-provided and spin-echo (SE) based convention gold standard (GS) method. Repeatability and reproducibility were characterized by the within-case coefficient of variation (wCV) of the MRF data acquired over 30 days, along with the biases.

RESULT

For the wide ranges of MRF estimated T₁ values, vials #1-8 (T₁ relaxation time between 2033 and 184 ms) exhibited a wCV less than 3% and 4%, respectively, on 3.0 and 1.5 T scanners. T₂ values in vials #1-8 (T₂ relaxation, 1044-45 ms) have shown wCV to be <7% on both 3.0 and 1.5 T scanners. A stronger linear correlation overall for T₁ (R²  = 0.9960 and 0.9963 at center-1 and center-2 on 3.0 T scanner, and R²  = 0.9951 and 0.9988 at center-1 and center-3 on 1.5 T scanner) compared to T₂ (R²  = 0.9971 and 0.9972 at center-1 and center-2 on 3.0 T scanner, and R²  = 0.9815 and 0.9754 at center-1 and center-3 on 1.5 T scanner). Bland-Altman (BA) analysis showed MRF based T₁ and T₂ values were within the limit of agreement (LOA) except for one data point. The mean difference or bias and 95% lower bound (LB) and upper bound (UB) LOA are reported in the format; mean bias: 95% LB LOA: 95% UB LOA. The biases for T₁ values were 21.34: -50.00: 92.69, 21.32: -47.29: 89.94 ms, and for T₂ values were -19.88: -42.37: 2.61, -19.06: -43.58: 5.45 ms on 3.0 T scanner at center-1 and center-2, respectively. Similarly, on 1.5 T scanner biases for T₁ values were 26.54: -53.41: 106.50, 9.997: -51.94: 71.94 ms, and for T₂ values were -23.84: -135.40: 87.76, -37.30: 134.30: 59.73 ms at center-1 and center-3, respectively. Additionally, the correlation between the SE based GS and MRF estimated T₁ and T₂ values (R²  = 0.9969 and 0.9977) showed a similar trend as we observed between vendor-provided and MRF estimated T₁ and T₂ values (R²  = 0.9963 and 0.9972). In addition to correlation, BA analysis showed that all the vials are within the LOA between the GS and vendor-provided for the T₁ values and except one vial for T₂ . All the vials are within the LOA between GS and MRF except one vial in T₁ and T₂ array. The wCV for reproducibility was <3% for both T₁ and T₂ values in vials #1-8, for all the 14 vials, wCV calculated for reproducibility was <4% for T₁ values and <5% for T₂ .

CONCLUSION

This study shows that MRF is highly repeatable (wCV <4% for T₁ and <7% for T₂ ) and reproducible (wCV < 3% for both T₁ and T₂ ) in certain vials (vials #1-8).

Three dimensional MRF obtains highly repeatable and reproducible multi-parametric estimations in the healthy human brain at 1.5T and 3T

Magnetic resonance fingerprinting (MRF) is highly promising as a quantitative MRI technique due to its accuracy, robustness, and efficiency. Previous studies have found high repeatability and reproducibility of 2D MRF acquisitions in the brain. Here, we have extended our investigations to 3D MRF acquisitions covering the whole brain using spiral projection k-space trajectories. Our travelling head study acquired test/retest data from the brains of 12 healthy volunteers and 8 MRI systems (3 systems at 3 T and 5 at 1.5 T, all from a single vendor), using a study design not requiring all subjects to be scanned at all sites. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived PD T1 and T2 maps to an anatomical atlas, coefficients of variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated a high repeatability (CVs 0.7-1.3% for T1, 2.0-7.8% for T2, 1.4-2.5% for normalized PD) and reproducibility (CVs of 2.0-5.8% for T1, 7.4-10.2% for T2, 5.2-9.2% for normalized PD) in gray and white matter. Both repeatability and reproducibility improved when compared to similar experiments using 2D acquisitions. Three-dimensional MRF obtains highly repeatable and reproducible estimations of T1 and T2, supporting the translation of MRF-based fast quantitative imaging into clinical applications.

Virtual Histology of Cortical Thickness and Shared Neurobiology in 6 Psychiatric Disorders

IMPORTANCE

Large-scale neuroimaging studies have revealed group differences in cortical thickness across many psychiatric disorders. The underlying neurobiology behind these differences is not well understood.

OBJECTIVE

To determine neurobiologic correlates of group differences in cortical thickness between cases and controls in 6 disorders: attention-deficit/hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BD), major depressive disorder (MDD), obsessive-compulsive disorder (OCD), and schizophrenia.

DESIGN, SETTING, AND PARTICIPANTS

Profiles of group differences in cortical thickness between cases and controls were generated using T1-weighted magnetic resonance images. Similarity between interregional profiles of cell-specific gene expression and those in the group differences in cortical thickness were investigated in each disorder. Next, principal component analysis was used to reveal a shared profile of group difference in thickness across the disorders. Analysis for gene coexpression, clustering, and enrichment for genes associated with these disorders were conducted. Data analysis was conducted between June and December 2019. The analysis included 145 cohorts across 6 psychiatric disorders drawn from the ENIGMA consortium. The numbers of cases and controls in each of the 6 disorders were as follows: ADHD: 1814 and 1602; ASD: 1748 and 1770; BD: 1547 and 3405; MDD: 2658 and 3572; OCD: 2266 and 2007; and schizophrenia: 2688 and 3244.

MAIN OUTCOMES AND MEASURES

Interregional profiles of group difference in cortical thickness between cases and controls.

RESULTS

A total of 12 721 cases and 15 600 controls, ranging from ages 2 to 89 years, were included in this study. Interregional profiles of group differences in cortical thickness for each of the 6 psychiatric disorders were associated with profiles of gene expression specific to pyramidal (CA1) cells, astrocytes (except for BD), and microglia (except for OCD); collectively, gene-expression profiles of the 3 cell types explain between 25% and 54% of variance in interregional profiles of group differences in cortical thickness. Principal component analysis revealed a shared profile of difference in cortical thickness across the 6 disorders (48% variance explained); interregional profile of this principal component 1 was associated with that of the pyramidal-cell gene expression (explaining 56% of interregional variation). Coexpression analyses of these genes revealed 2 clusters: (1) a prenatal cluster enriched with genes involved in neurodevelopmental (axon guidance) processes and (2) a postnatal cluster enriched with genes involved in synaptic activity and plasticity-related processes. These clusters were enriched with genes associated with all 6 psychiatric disorders.

CONCLUSIONS AND RELEVANCE

In this study, shared neurobiologic processes were associated with differences in cortical thickness across multiple psychiatric disorders. These processes implicate a common role of prenatal development and postnatal functioning of the cerebral cortex in these disorders.

Autism Spectrum Disorder and Childhood Apraxia of Speech: Early Language-Related Hallmarks across Structural MRI Study

Autism Spectrum Disorder (ASD) and Childhood Apraxia of Speech (CAS) are developmental disorders with distinct diagnostic criteria and different epidemiology. However, a common genetic background as well as overlapping clinical features between ASD and CAS have been recently reported. To date, brain structural language-related abnormalities have been detected in both the conditions, but no study directly compared young children with ASD, CAS and typical development (TD). In the current work, we aim: (i) to test the hypothesis that ASD and CAS display neurostructural differences in comparison with TD through morphometric Magnetic Resonance Imaging (MRI)-based measures (ASD vs. TD and CAS vs. TD); (ii) to investigate early possible disease-specific brain structural patterns in the two clinical groups (ASD vs. CAS); (iii) to evaluate predictive power of machine-learning (ML) techniques in differentiating the three samples (ASD, CAS, TD). We retrospectively analyzed the T1-weighted brain MRI scans of 68 children (age range: 34-74 months) grouped into three cohorts: (1) 26 children with ASD (mean age ± standard deviation: 56 ± 11 months); (2) 24 children with CAS (57 ± 10 months); (3) 18 children with TD (55 ± 13 months). Furthermore, a ML analysis based on a linear-kernel Support Vector Machine (SVM) was performed. All but one brain structures displayed significant higher volumes in both ASD and CAS children than TD peers. Specifically, ASD alterations involved fronto-temporal regions together with basal ganglia and cerebellum, while CAS alterations are more focused and shifted to frontal regions, suggesting a possible speech-related anomalies distribution. Caudate, superior temporal and hippocampus volumes directly distinguished the two conditions in terms of greater values in ASD compared to CAS. The ML analysis identified significant differences in brain features between ASD and TD children, whereas only some trends in the ML classification capability were detected in CAS as compared to TD peers. Similarly, the MRI structural underpinnings of two clinical groups were not significantly different when evaluated with linear-kernel SVM. Our results may represent the first step towards understanding shared and specific neural substrate in ASD and CAS conditions, which subsequently may contribute to early differential diagnosis and tailoring specific early intervention.

7-T MRI of explanted liver and ex-vivo pancreatic specimens: prospective study protocol of radiological-pathological correlation feasibility (the EXLIPSE project)

The study focuses on radiological-pathological correlation between imaging of ex vivo samples obtained by a 7-T scanner and histological examination. The specimens will be derived from native explanted cirrhotic livers, liver grafts excluded from donation because of severe steatosis, and primary pancreatic tumours. Magnetic resonance imaging (MRI) examinations will be performed within 24 h from liver or pancreatic lesion surgical removal. The MRI protocol will include morphological sequences, quantitative T1, T2, and fat-, water-fraction maps with Cartesian k-space acquisition, and multiparametric methods based on a transient-state “MRI fingerprinting”. Finally, the specimen will be fixed by formalin. Qualitative imaging analysis will be performed by two independent blinded radiologists to assess image consistency score. Quantitative analysis will be performed by drawing regions of interest on different tissue zones to measure T1 and T2 relaxation times as well as fat- and water-fraction. The same tissue areas will be analysed by the pathologists. This study will provide the possibility to improve our knowledge about qualitative and quantitative abdominal imaging assessment at 7 T, by correlating imaging characteristics and the corresponding histological composition of ex vivo specimens, in order to identify imaging biomarkers. Trial registration: ClinicalTrials.gov: 13646. Registered 9 July 2019—retrospectively registered

Reorganization of action observation and sensory-motor networks after action observation therapy in children with congenital hemiplegia: A pilot study

New rehabilitation programs based on action observation therapy (AOT) are effective in improving motor function in children with congenital hemiplegia. In this pilot study we tested the potential effects of AOT on the reorganization of the motor system by functional magnetic resonance imaging (fMRI). As part of a randomized trial, eight subjects (age range: 6.2-14.5 years) with congenital hemiplegia were randomly assigned to an experimental (EG) or control (CG) group. All children underwent a clinical and neurophysiological assessment with Assisting Hand Assessment (AHA), MRI, and fMRI at baseline (T0), 1(T1), and 8(T2) weeks after the end of 3-week treatment. For the EG, AOT consisted in the observation of uni/bimanual goal-directed actions followed by their execution. CG watched same-duration computer games and then performed the same actions in the same order used in the EG. fMRI study was carried out using two different paradigms, for exploring sensory-motor network (SMN) localization and action observation network (AON). The pattern of brain activation was generally similar between T0 and T1 for both groups, while it was more widespread at T2, compared to T0 and T1, in the EG. This enlargement was coupled with functional improvement at AHA. Single-subject analysis shows a reduction of lateralization indexes both for the AON and the SMN. This pilot study, despite the small sample, showed the fMRI feasibility for providing relevant biomarkers of brain plasticity for monitoring the AOT response in children with congenital hemiplegia. The study was registered at http://www.clinicaltrials.gov (identifier NCT01016496).

Subcortical Brain Volume, Regional Cortical Thickness, and Cortical Surface Area Across Disorders: Findings From the ENIGMA ADHD, ASD, and OCD Working Groups

OBJECTIVE

Attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), and obsessive-compulsive disorder (OCD) are common neurodevelopmental disorders that frequently co-occur. The authors sought to directly compare these disorders using structural brain imaging data from ENIGMA consortium data.

METHODS

Structural T1-weighted whole-brain MRI data from healthy control subjects (N=5,827) and from patients with ADHD (N=2,271), ASD (N=1,777), and OCD (N=2,323) from 151 cohorts worldwide were analyzed using standardized processing protocols. The authors examined subcortical volume, cortical thickness, and cortical surface area differences within a mega-analytical framework, pooling measures extracted from each cohort. Analyses were performed separately for children, adolescents, and adults, using linear mixed-effects models adjusting for age, sex, and site (and intracranial volume for subcortical and surface area measures).

RESULTS

No shared differences were found among all three disorders, and shared differences between any two disorders did not survive correction for multiple comparisons. Children with ADHD compared with those with OCD had smaller hippocampal volumes, possibly influenced by IQ. Children and adolescents with ADHD also had smaller intracranial volume than control subjects and those with OCD or ASD. Adults with ASD showed thicker frontal cortices compared with adult control subjects and other clinical groups. No OCD-specific differences were observed across different age groups and surface area differences among all disorders in childhood and adulthood.

CONCLUSIONS

The study findings suggest robust but subtle differences across different age groups among ADHD, ASD, and OCD. ADHD-specific intracranial volume and hippocampal differences in children and adolescents, and ASD-specific cortical thickness differences in the frontal cortex in adults, support previous work emphasizing structural brain differences in these disorders.

Rapid three-dimensional multiparametric MRI with quantitative transient-state imaging

Novel methods for quantitative, transient-state multiparametric imaging are increasingly being demonstrated for assessment of disease and treatment efficacy. Here, we build on these by assessing the most common Non-Cartesian readout trajectories (2D/3D radials and spirals), demonstrating efficient anti-aliasing with a k-space view-sharing technique, and proposing novel methods for parameter inference with neural networks that incorporate the estimation of proton density. Our results show good agreement with gold standard and phantom references for all readout trajectories at 1.5 T and 3 T. Parameters inferred with the neural network were within 6.58% difference from the parameters inferred with a high-resolution dictionary. Concordance correlation coefficients were above 0.92 and the normalized root mean squared error ranged between 4.2 and 12.7% with respect to gold-standard phantom references for T1 and T2. In vivo acquisitions demonstrate sub-millimetric isotropic resolution in under five minutes with reconstruction and inference times < 7 min. Our 3D quantitative transient-state imaging approach could enable high-resolution multiparametric tissue quantification within clinically acceptable acquisition and reconstruction times.

Musculoskeletal MRI at 7 T: do we need more or is it more than enough?

Ultra-high field magnetic resonance imaging (UHF-MRI) provides important diagnostic improvements in musculoskeletal imaging. The higher signal-to-noise ratio leads to higher spatial and temporal resolution which results in improved anatomic detail and higher diagnostic confidence. Several methods, such as T2, T2*, T1rho mapping, delayed gadolinium-enhanced, diffusion, chemical exchange saturation transfer, and magnetisation transfer techniques, permit a better tissue characterisation. Furthermore, UHF-MRI enables in vivo measurements by low-γ nuclei (²³Na, ³¹P, ¹³C, and ³⁹K) and the evaluation of different tissue metabolic pathways. European Union and Food and Drug Administration approvals for clinical imaging at UHF have been the first step towards a more routinely use of this technology, but some drawbacks are still present limiting its widespread clinical application. This review aims to provide a clinically oriented overview about the application of UHF-MRI in the different anatomical districts and tissues of musculoskeletal system and its pros and cons. Further studies are needed to consolidate the added value of the use of UHF-MRI in the routine clinical practice and promising efforts in technology development are already in progress.

Retrospective rigid motion correction of three-dimensional magnetic resonance fingerprinting of the human brain

PURPOSE

To obtain three-dimensional (3D), quantitative and motion-robust imaging with magnetic resonance fingerprinting (MRF).

METHODS

Our acquisition is based on a 3D spiral projection k-space scheme. We compared different orderings of trajectory interleaves in terms of rigid motion-correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7-s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole-brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k-space data, and the aligned data were matched with the dictionary to obtain motion-corrected maps.

RESULTS

A significant improvement on the motion-affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T₁ and T₂ estimations after motion correction. In addition, the average motion-induced quantification bias of 70 ms for T₁ and 18 ms for T₂ values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations.

CONCLUSION

We established a method that allows correcting 3D rigid motion on a 7-s timescale during the reconstruction of MRF data using self-navigators, improving the image quality and the quantification robustness.

European Ultrahigh-Field Imaging Network for Neurodegenerative Diseases (EUFIND)

Introduction

The goal of European Ultrahigh-Field Imaging Network in Neurodegenerative Diseases (EUFIND) is to identify opportunities and challenges of 7 Tesla (7T) MRI for clinical and research applications in neurodegeneration. EUFIND comprises 22 European and one US site, including over 50 MRI and dementia experts as well as neuroscientists.

Methods

EUFIND combined consensus workshops and data sharing for multisite analysis, focusing on 7 core topics: clinical applications/clinical research, highest resolution anatomy, functional imaging, vascular systems/vascular pathology, iron mapping and neuropathology detection, spectroscopy, and quality assurance. Across these topics, EUFIND considered standard operating procedures, safety, and multivendor harmonization.

Results

The clinical and research opportunities and challenges of 7T MRI in each subtopic are set out as a roadmap. Specific MRI sequences for each subtopic were implemented in a pilot study presented in this report. Results show that a large multisite 7T imaging network with highly advanced and harmonized imaging sequences is feasible and may enable future multicentre ultrahigh-field MRI studies and clinical trials.

Discussion

The EUFIND network can be a major driver for advancing clinical neuroimaging research using 7T and for identifying use-cases for clinical applications in neurodegeneration.

Multi-site repeatability and reproducibility of MR fingerprinting of the healthy brain at 1.5 and 3.0 T

Fully-quantitative MR imaging methods are useful for longitudinal characterization of disease and assessment of treatment efficacy. However, current quantitative MRI protocols have not been widely adopted in the clinic, mostly due to lengthy scan times. Magnetic Resonance Fingerprinting (MRF) is a new technique that can reconstruct multiple parametric maps from a single fast acquisition in the transient state of the MR signal. Due to the relative novelty of this technique, the repeatability and reproducibility of quantitative measurements obtained using MRF has not been extensively studied. Our study acquired test/retest data from the brains of nine healthy volunteers, each scanned on five MRI systems (two at 3.0 T and three at 1.5 T, all from a single vendor) located at two different centers. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived M0, T1 and T2 maps to an anatomical atlas, coefficients-of-variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated an excellent repeatability (CVs of 2-3% for T1, 5-8% for T2, 3% for normalized-M0) and a good reproducibility (CVs of 3-8% for T1, 8-14% for T2, 5% for normalized-M0) in grey and white matter.

Echo-time dependency of quantitative susceptibility mapping reproducibility at different magnetic field strengths

Quantitative Susceptibility Mapping (QSM) provides a way of measuring iron concentration and myelination non-invasively and has the potential of becoming a tool of paramount importance in the study of a host of different pathologies. However, several experimental factors and the physical properties of magnetic susceptibility (χ) can impair the reliability of QSM, and it is therefore essential to assess QSM reproducibility for repeated acquisitions and different field strength. In particular, it has recently been demonstrated that QSM measurements strongly depend on echo time (TE): the same tissue, measured on the same scanner, exhibits different apparent frequency shifts depending on the TE used. This study aims to assess the influence of TE on intra-scanner and inter-scanner reproducibility of QSM, by using MRI systems operating at 3T and 7T. To maximize intra-scanner reproducibility it is necessary to match the TEs of the acquisition protocol, but the application of this rule leads to inconsistent QSM values across scanners operating at different static magnetic field. This study however demonstrates that, provided a careful choice of acquisition parameters, and in particular by using TEs at 3T that are approximately 2.6 times longer than those at 7T, highly reproducible whole-brain χ maps can be achieved also across different scanners, which renders QSM a suitable technique for longitudinal follow-up in clinical settings and in multi-center studies.

Magnetic resonance fingerprinting with dictionary-based fat and water separation (DBFW MRF): A multi-component approach

Purpose

To obtain a fast and robust fat‐water separation with simultaneous estimation of water T₁, fat T₁, and fat fraction maps.

Methods

We modified an MR fingerprinting (MRF) framework to use a single dictionary combination of a water and fat dictionary. A variable TE acquisition pattern with maximum TE = 4.8 ms was used to increase the fat–water separability. Radiofrequency (RF) spoiling was used to reduce the size of the dictionary by reducing T₂ sensitivity. The technique was compared both in vitro and in vivo to an MRF method that incorporated 3‐point Dixon (DIXON MRF), as well as Cartesian IDEAL with different acquisition parameters.

Results

The proposed dictionary‐based fat–water separation technique (DBFW MRF) successfully provided fat fraction, water, and fat T₁, B₀, and B₁₊ maps both in vitro and in vivo. The fat fraction and water T₁ values obtained with DBFW MRF show excellent agreement with DIXON MRF as well as with the reference values obtained using a Cartesian IDEAL with a long TR (concordance correlation coefficient: 0.97/0.99 for fat fraction–water T₁). Whereas fat fraction values with Cartesian IDEAL were degraded in the presence of T₁ saturation, MRF methods successfully estimated and accounted for T₁ in the fat fraction estimates.

Conclusion

The DBFW MRF technique can successfully provide T₁ and fat fraction quantification in under 20 s per slice, intrinsically correcting T₁ biases typical of fast Dixon techniques. These features could improve the diagnostic quality and use of images in presence of fat.

Response to short-term deprivation of the human adult visual cortex measured with 7T BOLD

Sensory deprivation during the post-natal ‘critical period’ leads to structural reorganization of the developing visual cortex. In adulthood, the visual cortex retains some flexibility and adapts to sensory deprivation. Here we show that short-term (2 hr) monocular deprivation in adult humans boosts the BOLD response to the deprived eye, changing ocular dominance of V1 vertices, consistent with homeostatic plasticity. The boost is strongest in V1, present in V2, V3 and V4 but absent in V3a and hMT+. Assessment of spatial frequency tuning in V1 by a population Receptive-Field technique shows that deprivation primarily boosts high spatial frequencies, consistent with a primary involvement of the parvocellular pathway. Crucially, the V1 deprivation effect correlates across participants with the perceptual increase of the deprived eye dominance assessed with binocular rivalry, suggesting a common origin. Our results demonstrate that visual cortex, particularly the ventral pathway, retains a high potential for homeostatic plasticity in the human adult.

The world around us changes all the time, and the brain must adapt to these changes. This process, known as neuroplasticity, peaks during development. Abnormal sensory input early in life can therefore cause lasting changes to the structure of the brain. One example of this is amblyopia or ‘lazy eye’. Infants who receive insufficient input to one eye – for example, because of cataracts – can lose their sight in that eye, even if the cataracts are later removed. This is because the brain reorganizes itself to ignore messages from the affected eye.

Does the adult visual system also show neuroplasticity? To explore this question, Binda, Kurzawski et al. asked healthy adult volunteers to lie inside a high-resolution brain scanner with a patch covering one eye. At the start of the experiment, roughly half of the brain’s primary visual cortex responded to sensory input from each eye. But when the volunteers removed the patch two hours later, this was no longer the case.

Some areas of the visual cortex that had previously responded to stimuli presented to the non-patched eye now responded to stimuli presented to the patched eye instead. The patched eye had also become more sensitive to visual stimuli. Indeed, these changes in visual sensitivity correlated with changes in brain activity in a pathway called the ventral visual stream. This pathway processes the fine details of images. Groups of neurons within this pathway that responded to stimuli presented to the patched eye were more sensitive to fine details after patching than before.

Visual regions of the adult brain thus retain a high degree of neuroplasticity. They adapt rapidly to changes in the environment, in this case by increasing their activity to compensate for a lack of input. Notably, these changes are in the opposite direction to those that occur as a result of visual deprivation during development. This has important implications because lazy eye syndrome is currently considered untreatable in adulthood.

Vascular Function Is Improved After an Environmental Enrichment Program

Environmental enrichment may slow cognitive decay possibly acting through an improvement in vascular function. Aim of the study was to assess the effects of a 7-month cognitive, social, and physical training program on cognitive and vascular function in patients with mild cognitive impairment. In a single-center, randomized, parallel-group study, 113 patients (age, 65–89 years) were randomized to multidomain training (n=55) or usual care (n=58). All participants underwent neuropsychological tests and vascular evaluation, including brachial artery flow-mediated dilation, carotid–femoral pulse wave velocity, carotid distensibility, and assessment of circulating hematopoietic CD34+ and endothelial progenitor cells. At study entry, an age-matched control group (n=45) was also studied. Compared with controls, patients had at study entry a reduced flow-mediated dilation (2.97±2.14% versus 3.73±2.06%; P =0.03) and hyperemic stimulus (shear rate area under the curve, 19.1±15.7 versus 25.7±15.1×10 −3 ; P =0.009); only the latter remained significant after adjustment for confounders ( P =0.03). Training improved Alzheimer disease assessment scale cognitive (training, 14.0±4.8 to 13.1±5.5; nontraining, 12.1±3.9 to 13.2±4.8; P for interaction visit×training=0.02), flow-mediated dilation (2.82±2.19% to 3.40±1.81%, 3.05±2.08% to 2.24±1.59%; P =0.006; P =0.023 after adjustment for diameter and shear rate area under the curve), and circulating hematopoietic CD34 + cells and prevented the decline in carotid distensibility (18.4±5.3 to 20.0±6.6, 23.9±11.0 to 19.5±7.1 Pa −1 ; P =0.005). The only clinical predictor of improvement of cognitive function after training was established hypertension. There was no correlation between changes in measures of cognitive and vascular function. In conclusion, a multidomain training program slows cognitive decline, especially in hypertensive individuals. This effect is accompanied by improved systemic endothelial function, mobilization of progenitor CD34 + cells, and preserved carotid distensibility.

Clinical Trial Registration—

URL: http://www.clinicaltrials.gov . Unique identifier: NCT01725178.

Cortical and Subcortical Brain Morphometry Differences Between Patients With Autism Spectrum Disorder and Healthy Individuals Across the Lifespan: Results From the ENIGMA ASD Working Group

OBJECTIVE

Neuroimaging studies show structural differences in both cortical and subcortical brain regions in children and adults with autism spectrum disorder (ASD) compared with healthy subjects. Findings are inconsistent, however, and it is unclear how differences develop across the lifespan. The authors investigated brain morphometry differences between individuals with ASD and healthy subjects, cross-sectionally across the lifespan, in a large multinational sample from the Enhancing Neuroimaging Genetics Through Meta-Analysis (ENIGMA) ASD working group.

METHOD

The sample comprised 1,571 patients with ASD and 1,651 healthy control subjects (age range, 2-64 years) from 49 participating sites. MRI scans were preprocessed at individual sites with a harmonized protocol based on a validated automated-segmentation software program. Mega-analyses were used to test for case-control differences in subcortical volumes, cortical thickness, and surface area. Development of brain morphometry over the lifespan was modeled using a fractional polynomial approach.

RESULTS

The case-control mega-analysis demonstrated that ASD was associated with smaller subcortical volumes of the pallidum, putamen, amygdala, and nucleus accumbens (effect sizes [Cohen's d], 0.13 to -0.13), as well as increased cortical thickness in the frontal cortex and decreased thickness in the temporal cortex (effect sizes, -0.21 to 0.20). Analyses of age effects indicate that the development of cortical thickness is altered in ASD, with the largest differences occurring around adolescence. No age-by-ASD interactions were observed in the subcortical partitions.

CONCLUSIONS

The ENIGMA ASD working group provides the largest study of brain morphometry differences in ASD to date, using a well-established, validated, publicly available analysis pipeline. ASD patients showed altered morphometry in the cognitive and affective parts of the striatum, frontal cortex, and temporal cortex. Complex developmental trajectories were observed for the different regions, with a developmental peak around adolescence. These findings suggest an interplay in the abnormal development of the striatal, frontal, and temporal regions in ASD across the lifespan.

TRIMAGE: A dedicated trimodality (PET/MR/EEG) imaging tool for schizophrenia

Simultaneous PET/MR/EEG (Positron Emission Tomography - Magnetic Resonance - Electroencephalography), a new tool for the investigation of neuronal networks in the human brain, is presented here within the framework of the European Union Project TRIMAGE. The trimodal, cost-effective PET/MR/EEG imaging tool makes use of cutting edge technology both in PET and in MR fields. A novel type of magnet (1.5T, non-cryogenic) has been built together with a PET scanner that makes use of the most advanced photodetectors (i.e., SiPM matrices), scintillators matrices (LYSO) and digital electronics. The combined PET/MR/EEG system is dedicated to brain imaging and has an inner diameter of 260 mm and an axial Field-of-View of 160 mm. It enables the acquisition and assessment of molecular metabolic information with high spatial and temporal resolution in a given brain simultaneously. The dopaminergic system and the glutamatergic system in schizophrenic patients are investigated via PET, the same physiological/pathophysiological conditions with regard to functional connectivity, via fMRI, and its electrophysiological signature via EEG. In addition to basic neuroscience questions addressing neurovascular-metabolic coupling, this new methodology lays the foundation for individual physiological and pathological fingerprints for a wide research field addressing healthy aging, gender effects, plasticity and different psychiatric and neurological diseases. The preliminary performances of two components of the imaging tool (PET and MR) are discussed. Initial results of the search of possible candidates for suitable schizophrenia biomarkers are also presented as obtained with PET/MR systems available to the collaboration.