Quantitative Neuroimaging
Neuro Monday, 17 May 2021

Oral Session - Quantitative Neuroimaging
Neuro
Monday, 17 May 2021 14:00 - 16:00
  • Mapping of Thalamic Matrix and Core Nuclei using QSM at 9.4 Tesla
    Vinod Jangid Kumar1, Klaus Scheffler1,2, Gisela E Hagberg1,2, and Wolfgang Grodd1
    1Max Planck Institute for Biological Cybernetics, Tuebingen, Germany, 2Biomedical Magnetic Resonance, University Hospital and Eberhard-Karl’s University, Tuebingen, Germany
    The preliminary results show lower diamagnetic and paramagnetic sources for the core nuclei. In contrast, the matrix nuclei showed observable higher values.
    Fig. 4: Histogram of the mean values for the diamagnetic and paramagnetic sources of all nuclei. Each Histogram shows values for the matrix nuclei in blue color and the core nuclei in red color. Note the lower values for the core nuclei in contrast to matrix nuclei indicating the more substantial contribution of diamagnetic and paramagnetic sources, i.e., iron, myelin, and calcium in the matrix nuclei.
    Fig. 1: Thalamus QSM at 9.4 T. A) QSM map (scaled between: -215.0546 to 208.2119) B1) Negative values: diamagnetic sources (i.e., Myelin, Calcium, etc.) B2) Positive values (Paramagnetic sources, i.e., Iron, etc.) Both the values are in parts per billion, i.e., ppb (magnetic strength irrelevant). The depicted color code for all views is hot. The negative Myelin map shows a homogeneous distribution of negative QSM values across the thalamus; in contrast, the positive map shows increased QSM values mainly at the posterior, lateral, and intralaminar, midline nuclei.
  • On Comparability and Reproducibility of Myelin Sensitive Imaging Techniques
    Tom Hilbert1,2,3, Lucas Soustelle4, Gian Franco Piredda1,2,3, Thomas Troalen5, Stefan Sommer6,7, Arun Joseph8,9,10, Reto Meuli2, Jean-Philippe Thiran2,3, Guillaume Duhamel4, Olivier M. Girard4, and Tobias Kober1,2,3
    1Advanced Clinical Imaging Technology (ACIT), Siemens Healthcare, Lausanne, Switzerland, 2Department of Radiology, Lausanne University Hospital (CHUV), Lausanne, Switzerland, 3LTS5, École Polytechnique Fédérale de Lausanne, Lausanne, Switzerland, 4Aix Marseille Univ, CNRS, CRMBM, Marseille, France, 5Siemens Healthcare SAS, Saint-Denis, France, 6Siemens Healthcare, Zurich, Switzerland, 7Swiss Center for Musculoskeletal Imaging (SCMI), Balgrist Campus, Zurich, Switzerland, 8Advanced Clinical Imaging Technology (ACIT), Siemens Healthcare, Bern, Switzerland, 9Translational Imaging Center, Sitem-Insel, Bern, Switzerland, 10Departments of Radiology and Biomedical Research, University of Bern, Bern, Switzerland
    Longitudinal relaxation rate and magnetization transfer-based methods showed good correlation between each other. Myelin water fraction and ultra-short-echo-time imaging appears to show a different contrast. All methods showed good reproducibility.
    Figure 1: Example images from one subject in the three orthogonal views. The red arrows indicate the genu of the corpus callosum where in some methods a common hyperintensity was observed. R1 – longitudinal relaxation rate, MTR – Magnetization Transfer Ratio, MPF – Macromolecular Proton Fraction, ihMTSat – Inhomogeneous magnetization transfer saturation, MWF – Myelin Water Fraction, IR-UTE – Inversion Recovery Ultra-Short Time-to-Echo.
    Figure 3: Agreement and Bland-Altman plots for each method that show the reproducibility of each method. Colours indicate different tissue types: blue - brain structures, red - ventricles, yellow - gray matter, purple - white matter.
  • Delineating perfusion and the effects of vascularisation patterns across the hippocampal subfields at 7T
    Roy AM Haast1, Sriranga Kashyap2, Mohamed D Yousif1, Dimo Ivanov2, Benedikt A Poser2, and Ali R Khan1
    1Centre for Functional and Metabolic Mapping, Western University, London, ON, Canada, 2Maastricht University, Maastricht, Netherlands
    In this study we used arterial spin labeling combined with time-of-flight angiography acquired at 7T to delineate perfusion and the effects of vascularisation patterns across the hippocampal subfields.
    Across subjects average (A) perfusion-weighted (left) and perfusion (right) flatmaps. Subfield-specific averages, per subject and hemisphere, are shown in B for both perfusion-weighted (left) and perfusion metrics (right). Semi-transparent points show vertex-wise averages. Across subjects average perfusion-weighted signal (left) and standard variation (right), as a function of number of included runs and separated per subfield (colored lines) are displayed in C.
    Correlation analyses between vertex-wise perfusion-weighted average across-subjects (A, y-axis), standard variation (B, y-axis) and distance (x-axes). Marginal plots show distribution for the metric displayed along that axis. Points are color-coded based on the diameter of most nearby vessel. Overlaid solid and dashed lines indicate corresponding best fits for each category of diameter.
  • In vivo human T2* imaging at 0.35 mm reveals up to 15 ms of local variations within gray matter across depths at 7T
    Omer Faruk Gulban1,2, Saskia Bollman3, Renzo Huber1, Kendrick Kay4, Benedikt Poser1, Federico De Martino1, and Dimo Ivanov1
    1Department of Cognitive Neuroscience , Faculty of Psychology and Neuroscience, Maastricht University, Maastricht, Netherlands, 2Brain Innovation, Maastricht, Netherlands, 3Centre for Advanced Imaging, The University of Queensland, Brisbane, Australia, 4Center for Magnetic Resonance Research, University of Minnesota, Minneapolis, MN, United States
    • We measured in vivo human brain T2* values using 7T MRI at 0.35 × 0.35 × 0.35 mm3 using MRI at 7 T.
    • Gray matter T2* locally varies up to 15 ms within Heschl’s gyrus & Calcarine sulcus across cortical depths.
    • T2* seems more affected by biological tissue compositions than local cortex alignment to B0.
    Figure 2. Top row shows a zoomed-in image around Heschl's gyrus. Bottom row shows a zoomed-in image around the calcarine sulcus. Left column shows MP2RAGE UNI contrasts and right column shows T2* measurements obtained from ME-GRE.
    Figure 3. A) Showing steps for generating our regions of interest centered around the primary auditory and visual cortices. B) showing steps for the depth metrics and B0 alignment.
  • Construction of Spatiotemporal Cortical Surface Atlases for Fetal Brains
    Zhengwang Wu1, Yuchen Pei1, Ya Wang1, Tao Zhong1, Fenqiang Zhao1, Li Wang1, He Zhang2, and Gang Li1
    1Department of Radiology and BRIC, UNC-Chapel Hill, Chapel Hill, NC, United States, 2Department of Radiology, Obstetrics and Gynecology Hospital, Fu Dan University, Shanghai, China
    We constructed a set of temporally-densely sampled cortical surface atlases for characterizing the dynamic fetal brain cortex expanding from the 22 to 36 gestational weeks.
    Fig. 1. The constructed fetal cortical surface atlases (left hemisphere) at different gestational weeks (GW). The number in the bracket indicates the subject number for construction the atlas at the time point. (a) The sulcal depth on the spherical surface. (b) The average convexity on the spherical surface. (c) The mean curvature on the spherical surface. (d) The sulcal depth on the inner cortical surface. (e) The average convexity on the inner cortical surface. (f) The mean curvature on the inner cortical surface. (g) The Desikan parcellation on the inner cortical surface.
    Fig. 2. The cortical surface area of the atlases at different gestational weeks, in both the left and right hemisphere.
  • Combined blood flow and CO2-mediated effects underlie the tissue-specific response to hypercapnia: Insight from 7T MR-based imaging
    Allen A Champagne1,2 and Alex A Bhogal3
    1School of Medicine, Queen's University, Kingston, ON, Canada, 2Center for Neuroscience Studies, Queen's University, Kingston, ON, Canada, 3Radiology, University Medical Center Utrecht, Utrecht, Netherlands
    The cerebrovascular response to hypercapnia is driven by factors including CO2 sensitivity, blood flow (re)distribution effects, which together, emphasize the importance temporal analysis for improving the management of vascular brain diseases.

    Figure 2. Tissue-based distribution of lag parameters for each respiratory design

    (A) Refence anatomical axial slices in MNI space. (B-D) Averaged lag maps (seconds). The cumulative percent frequency (normalized to 100%) for the distribution of lag (seconds) is shown for each tissue which was extracted using the grey- (black) and white- (grey) matter mask displayed in (B), bottom right corner. A dotted red line was added to each histogram (B-D) at 40 seconds, for reference and comparison.

    Figure 1. Grey- and white-matter average CVR lag times for each stimulus
  • 7T QSM guided Histologically Consistent Thalamic Sub-nucleus Parcellation in 3T QSM Atlas Space
    Weimin Zhang1, Chenyu He1, Xiaojun Guan2, Xiaojun Xu2, Hongjiang Wei3, and Yuyao Zhang1,4,5
    1School of Information Science and Technology, ShanghaiTech University, Shanghai, China, 2Department of Radiology, The Second Affiliated Hospital, Zhejiang, China, 3Institute for Medical Imaging Technology, School of Biomedical Engineering, Shanghai, China, 4Shanghai Engineering Research Center of Intelligent Vision and Imaging, ShanghaiTech University, Shanghai, China, 5iHuman Institute, ShanghaiTech University, Shanghai, China
    we present a multi-atlas label fusion framework to delineate a total number of 64 thalamic nuclei using both 7T and 3T QSM images.   By combining group-wise registration, label propagation, and label fusion processes, we achieve a 64 thalamic nuclei parcellation map in 3T QSM atlas space.
    Fig. 1. Flowchart of multi-atlas label fusion framework. Firstly, a non-linear registration is performed to align 3T individual to 7T QSM image space; then inverted deformation fields are applied onto 7T thalamic sub-nucleus parcellation map to generate individual estimation maps; finally, joint label fusion method fuses these parcellation map into 3T atlas space.
    Fig. 3. Thalamic parcellation map on a representative slice of 3T atlas image shown in axial (top row), sagittal (middle row), and coronal (bottom row) views. The Schaltenbrand and Wahren histology atlas references are shown on the right.
  • Whole-Brain 3D Quantitative BOLD Mapping With Preliminary Estimation of R2' and Venous Blood Volume
    Hyunyeol Lee1 and Felix W Wehrli1
    1Radiology, University of Pennsylvania, Philadelphia, PA, United States
    Here, we develop a new whole-brain 3D quantitative BOLD parameter mapping method. The confounders  (R2, non-heme iron) were addressed by employing prior estimates of R2’ and cerebral venous blood volume. Results suggest feasibility of the proposed, prior-based whole-brain 3D qBOLD mapping.
    Figure 4. Whole-brain 3D maps of Yv, DBV, R′2,nb, and χnb in the three orthogonal planes, obtained using the proposed qBOLD method. Note the expected contrast and physiologically plausible range in all four parameters, i.e., a clear depiction of gray/white matter differentiation, relatively homogeneous Yv, and elevated R′2,nb and χnb in the deep brain structures relative to cortical regions.
    Figure 1. Schematic of the proposed, preliminary estimates-based qBOLD parameter mapping procedure.
  • 10min Whole-Brain MRI Solution – Comprehensive Quantification of MR Relaxometry and Susceptibility Plus Synthetic Contrast-Weighted Images
    Sen Ma1, Tianle Cao1,2, Nan Wang1, Anthony G. Christodoulou1, Zhaoyang Fan1, Yibin Xie1, and Debiao Li1
    1Biomedical Imaging Research Institute, Cedars-Sinai Medical Center, Los Angeles, CA, United States, 2Department of Bioengineering, University of California, Los Angeles, Los Angeles, CA, United States
    We propose an integrated and efficient solution to clinical brain MRI in a single 10min sequence, producing co-registered, quantitative PD, T1, T2, T1ρ, T2*, QSM, and ΔB0 plus synthetic weighted images including PDw, T1w, T2w, T2*w, FLAIR, SWI, true-SWI, mIP, and true-SWI mIP simultaneously.
    Figure 2. Axial PD, T1, T2, T1ρ, T2*, QSM, and ∆B0 maps of three slices using the proposed method.
    Figure 4. Synthetic PDw, T1w, T2w, T2*w, FLAIR, SWI, tSWI of three axial slices corresponding to Figure 2, as well as mIP and tmIP calculated from SWI and tSWI respectively using 16 axial slices. Contrast-weighted images are synthesized with good image quality, image contrasts, and venous structures.
  • Magnetic resonance recording of local neuronal firings (mrLNF) in the human brain: A proof of concept
    Yongxian Qian1, Karthik Lakshmanan1, Anli Liu2, Yvonne W. Lui1, and Fernando E. Boada1
    1Radiology, New York University, New York, NY, United States, 2Neurology, New York University, New York, NY, United States
    This human study at 3T tests a novel idea of magnetic resonance recording of local neuronal firings (mrLNF). This technique has a very high temporal resolution (0.25ms) and can non-invasively detect fast action potentials and slow postsynaptic neuronal currents at the firing sources.
    Fig. 1. Conceptual schematic of magnetic resonance recording of local neuronal firings (mrLNF). A) Localization of local neuronal firing volumes defined by multiple individual coil sensitivities (shading areas) or by selective single-voxel excitation (dashed line). B) Intrinsic micro quantum sensors of nuclear spins (short arrows) immersed in ionic flow (neuronal current) inside neurons. C) Modulation of neuronal firing magnetic field Bn onto MR main field B0. D) Recording of neuronal firing via MRI scanner.
    Fig. 3. The mrLNF recordings on a 35-year-old female subject’s head at resting state at 3T. A) Localization of the mrLNF recordings on sodium MRI images. B) Raw FID signal at Channel 0. C) The mrLNF recordings from single FID acquisition, showing a fast action potential firing at Channel 0 (arrow, 2ms in duration, 648.6 μT in intensity). D) The mrLNF recordings from three consecutive FID acquisitions.
Back to Top
Digital Poster Session - Quantitative Neuro: Technical Advances
Neuro
Monday, 17 May 2021 15:00 - 16:00
  • Precise localization of Deep brain nuclei in MNI-space guided by Hybrid Multi-modal MRI Brain Atlas
    Chenyu He1, Xiaojun Guan2, Boliang Yu1, Hongjiang Wei3, and Yuyao Zhang1,4,5
    1School of Information Science and Technology, ShanghaiTech University, Shanghai, China, 2Department of Radiology, The Second Affiliated Hospital, Zhejiang, University School of Medicine, Hangzhou, China, 3Institute for Medical Imaging Technology, School of Biomedical Engineering, Shanghai, Jiao Tong University, Shanghai, China, 4Shanghai Engineering Research Center of Intelligent Vision and Imaging, ShanghaiTech University, Shanghai, China, 5iHuman Institute, ShanghaiTech University, Shanghai, China
    Caudate nucleus (CN), putamen (Pu), ventral pallidum (VP), globus pallidus (GP, GPi & GPe), nucleus accumbens (NAC), subthalamic nucleus (STN), substantia nigra (SN, SNc & SNr) and red nucleus (RN) are preciously depicted in the MNI space in our work.
    Visualization of hand-crafted nuclei depictions. (a) Overall visualization of the parcellation map in axial, sagittal and coronal views along with a 3D rendered view (label information only given on right hemisphere); (b) Labeled CN, NAC and VP on the sections of three views; (c) Labeled Pu, GPe and GPi on the sections of three views; (d) Labeled STN, RN, SNr and SNc on the sections of three views.
    Pipeline of MNI space QSM atlas generation. (a) Special normalization for individual T1w and QSM image; and hybrid image generation. (b) Atlas space construction. (c) Guiding by T1w template, QSM atlas is projected into MNI space.
  • Multi-contrast Brain Image Registration for Low-contrast Paediatric Magnetic Resonance Imaging
    Cassandra M.J. Wannan1, Shane Tonissen2, Bruce Tonge3, Efstratios Skafidas 1,2, Christos Pantelis1,2, and Warda T. Syeda1
    1Melbourne Neuropsychiatry Centre, The University of Melbourne, Parkville, Australia, 2Department of Electrical and Electronic Engineering, The University of Melbourne, Parkville, Australia, 3Psychiatry Monash Health, Faculty of Medicine, Nursing and Health Sciences School of Clinical Sciences, Monash University, Melbourne, Australia

    - A multi-contrast image registration technique leverages information from multi-sequence MRI to non-linearly register an image to a reference space.

    - The proposed technique outperforms single-contrast registration methods in a sample of 6-months old healthy infants.  

    Figure 2: Exemplar skull-stripped T1w template images (top row) and corresponding single subject coronal slices from images registered using multi-contrast registration (second row), single-contrast ANTs (third row), FLIRT (fourth row) and FNIRT (last row). Single- and Multi-contrast images match template neuro-anatomy more closely in skull-stripped data compared to other methods.
    Figure 1: Exemplar T1w template coronal images in a single subject (top row) to visually compare multi-contrast registration (second row), single-contrast ANTs (third row) and FLIRT (bottom row) in whole-head data (reference space). Overlay: T1w template brain mask. Multi-contrast technique shows improved overlap between template and registered images. White arrows show underestimation of brain surface in single-contrast ANTs.
  • Comprehensive Human Brain Imaging Enhancement with the Single-frequency Excitation Wideband MRI (SE-WMRI) Technique
    Jordan Wang1, Po-Wei Cheng1, Ming-Jang Chiu2, Tzi-Dar Chiueh3, and Jyh-Horng Chen3
    1Graduate institute of biomedical electronics and bioinformatics, National Taiwan University, Taipei, Taiwan, 2National Taiwan University Hospital, Taipei, Taiwan, 3Department of electrical engineering, National Taiwan University, Taipei, Taiwan
    A comprehensive human brain imaging enhancement is achieved with the Single-frequency Excitation Wideband MRI technique on a Siemens PRISMA 3T human MRI system in terms of acceleration, SNR improvement, and resolution enhancement.
    Fig. 1 Conventional GRE and 2-time accelerated SE-GRE image comparison under the same resolution with the zoomed-in comparison of the center region. The two images have similar SNR, and their SSIM value is above 0.99, indicating an extremely high resemblance for image quality in terms of signal strength, contrast, and structural similarity.
    Fig.3 Standard GRE (Res. = 1 mm) and SE-GRE (Res. = 0.5 mm) images acquired within the same scan time and the zoomed-in comparisons of the bottom-left region and the top-right region. The comparison shows significant improvement for clarity of structural details, visibility of vessels, and edge sharpness between different tissues.
  • Deep Learning Enables 60% Accelerated Volumetric Brain MRI While Preserving Quantitative Performance – A Prospective, Multicenter Trial
    Suzie Bash1, Long Wang2, Chris Airriess3, Sara Dupont2, Greg Zaharchuk4, Enhao Gong2, Tao Zhang2, Ajit Shankaranarayanan2, and Lawrence Tanenbaum5
    1Neuroradiology, RadNet, Woodland Hills, CA, United States, 2Subtle Medical, Menlo Park, CA, United States, 3Cortechs.ai, San Diego, CA, United States, 4Stanford University, Stanford, CA, United States, 5RadNet, New York, NY, United States
    Deep learning can enable 60% faster brain MR examinations with matched clinical disease status predictability and statistically superior perceived image quality while maintaining high quantitative accuracy when compared with the longer standard of care exams.  
    FIG 3. Representative 3D T1W multiplanar images with volumetric segmentation on a 3T scanner. [Left to right]: Sagittal, coronal, axial T1W images with SOC (scan time 5:01 min) on the top row and FAST-DL (scan time 2:37 min) on bottom row.
    FIG 4. Representative axial 3D T1W images on a 3T scanner. [Left to right]: SOC (scan time 9:13 min), FAST (scan time 4:36 min), FAST-DL (scan time 4:36 min).
  • T2 Mapping of the Cranial Nerves with Multi-Interleaved X-prepared Turbo-spine Echo with Intuitive Relaxometry (MIXTURE) FLAIR
    Hajime Yokota1, Takayuki Sakai2, Masami Yoneyama3, Yansong Zhao4, and Takashi Uno1
    1Department of Diagnostic Radiology and Radiation Oncology, Chiba University, Chiba, Japan, 2Department of Radiology, Eastern Chiba Medical Center, Togane, Japan, 3Philips Japan, Tokyo, Japan, 4Philips Healthcare, Cleveland, OH, United States
    MIXTURE provided high-resolution morphology images and T2 map simultaneously and can measure T2 values with high reproducibility for the cranial nerve, which was impossible in conventional multi-echo TSE and MIXTURE T2.

    Figure 1. Scheme of the MIXTURE (Multi-Interleaved X-prepared tse with inTUitive RElaxometry)

    T2-mapping was performed using T2-prepared 3D segmented turbo spin-echo (TSE) with variable refocusing pulse trains. Two images with different TE (TE = 0 and 50ms) were acquired with interleaved acquisition. To obtain FLAIR contrast, inversion time (TI) is adjusted to suppress cerebrospinal fluid.

    Figure 3. Trigeminal nerve presentation of multi-echo TSE, MIXURE T2, and MIXTURE FLAIR

    Arrows indicates bilateral trigeminal nerves. The prepontine cistern around the trigeminal nerve looks noisy on T2 maps of multi-echo TSE and MIXURE T2, whereas the noise is reduced on the T2 map of MIXURE FLAIR.

  • Interleaved Black- and Bright-Blood Acquisition for Automatic Brain Metastasis Detection using Deep Learning Convolutional Neural Network
    Makoto Obara1, Yoshitomo Kikuchi2, Akio Hiwatashi2, Alexander Fischer 3, Yuta Akamine1, Tetsuo Ogino1, Masami Yoneyama1, Ronee Asad1, Yu Ueda1, Jihun Kwon1, and Marc Van Cauteren4
    1Philips Japan, Tokyo, Japan, 2Departments of Clinical Radiology, Graduate School of Medical Sciences, Kyushu University, Fukuoka, Japan, 3Philips GmbH Innovative Technologies, Aachen, Germany, 4Philips Healthcare, Tokyo, Japan
    We developed deep learning model for automatic brain metastasis detection using black- and bright-blood images for input data. The results suggest the usefulness of using two contrasts images as input, compared to single contrast input.

    Figure 3: Four representative segmentation cases by VISIBLE-DM:

    Segmented regions are shown in pink. Ground truth created by radiologists are also shown for comparison.

    Figure 2: Model development: a) Annotation examples (white arrows), b) Models are developed with three kinds of input data, black and bright, black only and bright only images.
  • Performance of data driven learned sampling patterns for accelerating brain 3D-T1ρ MRI
    Rajiv G Menon1, Marcelo V.W. Zibetti1, and Ravinder R. Regatte1
    1New York University Langone Health, New York, NY, United States
    Data-driven optimization of sampling pattern provided significant improvement in the performance of 3D-T mapping for brain applications. A significant reduction in the time required for the acquisition of 3D-T1ρ MRI data can be achieved.
    Figure 3: Comparison of SPs at different AFs. (a) shows FS k-space and SENSE reconstruction. (b) shows Poisson-disc SP at different AFs (4, 10, 20, 30) and the corresponding Low-rank reconstruction. (c) shows the optimized SP at the same AFs as (b) and resulting low-rank reconstructions. The improvement in performance is highlighted by the arrows.
    Figure 5: Comparison shows improved performance for optimized SP. (a) shows NRMSE errors across iterations for k-space training and validation sets (b) shows NRMSE errors across iterations for image training and validation sets (c) shows lower k-space NRMSE errors at different AFs for the optimized SP compared to PD SP, both SP using low-rank reconstruction (d) shows lower image NRMSE errors at different AFs for optimized SP vs PD (e) the improvement of the optimized SP vs PD in T mapping compared to FS reference at different AFs.
  • A Multiblock Partial Least Squares Correlation Framework for Covariate Adjustment and Interpretation of Latent Associations in Multimodal Data
    Warda T. Syeda1, Bjørn H. Ebdrup1,2,3, Cassandra M.J. Wannan1, Micah Cearns1, Rigas Soldatos1, Antonia Merritt1, Mahesh Jayaram 1, Andrew Zalesky 1, Jayachandra M. Raghava2,4, Birgitte Fagerlund 2, Egill Rostrup 2, Birte Glenthøj 2,3, Leigh A. Johnston5,6, Chad Bousman 1,7, Ian Everall8, Efstratios Skafidas 1,9, and Christos Pantelis1,9
    1Melbourne Neuropsychiatry Centre, The University of Melbourne, Parkville, Australia, 2Center for Neuropsychiatric Schizophrenia Research and Center for Clinical Intervention and Neuropsychiatric Schizophrenia Research, Mental Health Centre Glostrup, Copenhagen University Hospital, Glostrup, Denmark, 3Faculty of Health and Medical Sciences, Department of Clinical Medicine, University of Copenhagen, Denmark, Denmark, 4Functional Imaging Unit, Department of Clinical Physiology, Nuclear Medicine and PET, Rigshospitalet, Glostrup, Denmark, 5Department of Biomedical Engineering, The University of Melbourne, Melbourne, Australia, 6Melbourne Brain Centre Imaging Unit, The University of Melbourne, Parkville, Australia, 7Departments of Medical Genetics, Psychiatry, and Physiology & Pharmacology, University of Calgary, Calgary, AB, Canada, 8Institute of Psychiatry, Psychology & Neuroscience, King's College London, London, United Kingdom, 9Electrical and Electronic Engineering Department, The University of Melbourne, Parkville, Australia

    A multiblock PLS-C framework with covariate representation to perform multivariate statistical modelling is presented. 

    Latent structure-cognition associations in treatment resistant schizophrenia provide insights into widespread patterns of structural and cognitive deficits. 

    Figure 4: LV1 cognition-cortical volume saliences. A) Cortical volume saliences in TRS (red) and controls (green). Black lines: confidence intervals (CIs). Unreliable regions (CIs crossing zero) are grayed out. B) Cognitive saliences. y-axis: salience strength, x-axis: cognitive variables: PAL total errors (PAL-TE), stages completed (PAL-SC), IED interdimensional (IED-IS) and extradimensional shift (IED-ES), spatial-span length (SSP-SL), spatial working-memory strategy (SWM-S), total errors (SWM-TE). C) Reliable volume saliences projected onto a glass-brain.
    Figure 2: Decomposition of multi-block cross-correlation matrix into additive components using MB-PLS-C framework (x-axis: cognitive measures, y-axis: regional cortical volumes). The first two components corresponding to the significant latent variables, LV1 and LV2, describe correlations in the latent space between structure and cognitive measures across patients and healthy controls. Four disjoint data blocks with multivariate measures of brain structure, cognition and covariates.
  • Concordance of Regional Hypoperfusion by ASL MRI and 15O-water PET in Frontotemporal Dementia: Is ASL an Efficacious Alternative?
    Tracy Ssali1, Lucas Narciso1,2, Justin Hicks1,2, Matthais Günther3, Frank Prato1,2, Udunna Anazodo1,2, Elizabeth Finger4, and Keith St Lawrence1,2
    1Lawson Health Research Institute, London, ON, Canada, 2Department of Medical Biophysics, Western University, London, ON, Canada, 3Fraunhofer Institute for Medical Image Computing MEVIS, Bremen, Germany, 4Department of Clinical Neurological Sciences, Western University, London, ON, Canada
    This work highlights the potential of ASL for identifying regional hypoperfusion in FTD patients. While 15O-water PET data showed greater sensitivity, similar areas of hypoperfusion were identified by ASL, particularly for relative CBF maps which reduced inter-subject variability. 
    Figure 1: Mean control perfusion maps measured by ASL and 15O-water (whole brain CBF = 61.6 ± 12.4 ml/100g/min).
    Figure 2: T-maps generated using relative and absolute CBF measured by ASL and 15O-water in 4 exemplary FTD patients. Areas in red (ASL) and blue (15O-water) indicate regions of significant hypoperfusion in the patient participant compared to the control group. sv = semantic variant, bv=behavioural variant, nfPPA = non-fluent primary progressive aphasia, psp = progressive supranuclear palsy.
  • The relationship between Cerebrovascular Reactivity and baseline Cerebral Blood Flow: the effect of acquisition and analysis choices
    Rachael C Stickland1, Kristina M Zvolanek1,2, Stefano Moia3,4, Apoorva Ayyagari1,2, César Caballero-Gaudes3, and Molly G Bright1,2
    1Physical Therapy and Human Movement Sciences, Feinberg School of Medicine, Northwestern University, Chicago, IL, United States, 2Biomedical Engineering, McCormick School of Engineering, Northwestern University, Evanston, IL, United States, 3Basque Center on Cognition, Brain and Language, Donostia, Spain, 4University of the Basque Country EHU/UPV, Donostia, Spain
    Our results suggest that the relationship between baseline Cerebral Blood Flow and BOLD Cerebral Vascular Reactivity (CVR), across space and across people, is stronger when CVR is optimized for vascular delay and modelled with data including breathing tasks.
    Example data files for one subject, in T1-weighted space. Top panel shows the Harvard-Oxford Cortical atlas: 48 regions split into left and right hemisphere parcels. The middle panel shows CVR maps, not optimized for lag (No-Opt) and optimized for lag (Lag-Opt), for the BH+REST data segment. The bottom panel shows baseline CBF maps. Average CBF and CVR were obtained from each parcel.
    The association between bCBF and CVR values in 96 cortical parcels, compared across each of the five data segments (columns), for CVR values with no lag optimization (No-Opt) and CVR values with lag optimization (Opt). The thick lines that connect filled dots represent the group mean, and the thin lines represent each individual subject. The shaded grey box represents which correlations were not significant at the single subject level, based on a critical value of 1.96 (p<0.05, two-tailed).
  • Image-based self-gating for motion artifact free imaging of the eye
    Kilian Stumpf1, Hanna Frantz1, Patrick Metze1, Thomas Hüfken1, Tobias Speidel1, and Volker Rasche1
    1Department of Internal Medicine II, Ulm University Medical Center, Ulm, Germany
    An image-based self-gating approach is applied to MR acquisitions of the eye using tiny golden angle profile ordering, allowing the reconstruction of images free of motion artifacts without the use of external tracking devices.
    Figure 1: A) SG signal derived from the tracking of the ON position and the ON diameter for both eyes. The bins for two ON positions are marked in yellow with the reconstructed images being presented in Fig. 3. B) Scout image in axial orientation with visible ON (solid arrow), EOMS (dashed arrows) and orbital fat (dotted arrows) and two exemplary lines (green) on which the location of the ON is calculated and tracked.
    Figure 3: Images resulting from reconstruction of profiles (A) without SG, (B) with SG from the bin for centre location of the ON and (C) with SG from the bin for ‘looking left’ as determined by the SG signal in Fig 1. While motion blur of the ON (solid arrow) is clearly visible in A, the image-based SG approach enables the reconstruction of images in various motion phases with reduced motion blur (solid arrows in B,C) and allows e.g. the visualization of an EOM contraction (striped arrow in C).
  • Combined cluster analysis of time evolution and tissue type with total variation denoising (CCTV) for QQ-based oxygen extract fraction mapping
    Junghun Cho1, Pascal Spincemaille1, Thanh D Nguyen1, Ajay Gupta1, and Yi Wang1,2
    1Radiology, Weill Cornell Medicine, New York, NY, United States, 2Biomedical Engineering, Cornell University, Ithaca, NY, United States
    A Combined Cluster analysis and Total Variation denoising algorithm (CCTV) for oxygen extraction fraction (OEF) maps based on QSM+qBOLD (QQ) improved detection of lesion OEF abnormality in stroke patients compared to cluster analysis only without total variation denoising.
    Comparison between the OEF obtained by QQ-CAT and QQ-CCTV in 6 stroke patients imaged between 6hrs and 10 days post stroke onset. QQ-CCTV generally shows more uniform OEF maps compared to QQ-CAT. In 7 and 9 days post-onset patients, low OEF areas in QQ-CCTV agree better with DWI-defined lesions.
    Figure 1. Comparison between the OEF obtained by QQ-CAT and QQ-CCTV in a simulated stroke dataset. The numbers indicate RMSE (yellow) and MSD (white). On average, QQ-CCTV provide greater accuracy (MAE: 4.5 and 3.5%) and precision (MSD: 5.1 vs 1.7%). The OEFavg and OEFstd indicates the average and standard deviation OEF map among 5 trials, respectively.
  • Impact of b-value on the estimation of white matter fiber orientation dependent R2*
    Melanie Bauer1,2, Celine Berger1,2, Claudia Lenz1,2, Eva Scheurer1,2, and Christoph Birkl3
    1Institute of Forensic Medicine, Biomedical Engineering, University of Basel, Basel, Switzerland, 2Institute of Forensic Medicine, Health Department Basel-Stadt, Basel, Switzerland, 3Department of Neuroradiology, Medical University of Innsbruck, Innsbruck, Austria
    The orientation dependency of R2* is not influenced by the chosen b-values for post mortem in situ and in vivo MRI scans.
    Figure 1: R2* behavior in relation to increasing angle of white matter fibers. Post mortem data recorded with b = 2000 s/mm2 (red), b = 6000 s/mm2 (green) and b = 9000 s/mm2 (yellow) are marked with dots and in vivo data recorded with b = 800 s/mm2 (blue) and b = 2000 s/mm2 (red) are marked with triangles. The shaded areas represent the 95% confidence interval of the respective data.
    Figure 2: Boxplots for MD (left) and FA (right) of different b-values used for in vivo (red) and post mortem (blue) scans. Statistical significant differences < 0.05 are marked with *, < 0.01 with + and < 0.001 with x.
  • NORDIC denoising before image reconstruction.
    Steen Moeller1, Cheryl Olman2, Luca vizioli1, Logan Dowdle1, Essa Yacoub1, Mehmet Akcakaya1,3, and Kamil Ugurbil1
    1University of Minnesota, MINNEAPOLIS, MN, United States, 2Psychology, University of Minnesota, MINNEAPOLIS, MN, United States, 3ELECTRICAL AND COMPUTER ENGINEERING, University of Minnesota, Minneapolis, MN, United States
     NORDIC (LLR PCA) is directly compatible with undersampled k-space data
    Figure 4. GLM for a 4 block-design scan. Voxels significantly modulated by stimulus presentation in the occipital cortex with (p < 0.001) after cluster-based correction for multiple comparisons. Statistical analysis was performed with a GLM from 4 block-design scans (11 12-sec blocks of visual stim interleaved with 10 12-sec blocks of rest). The GLM regressors were formed by convolution of boxcar representing stimulus presentation with a simple gamma function $$$g(t) = t^q \cdot e^{-t} / (q^q \cdot e^{-q}), q=4 $$$
    Figure 3. Channel dependent impact of NORDIC applied before GRAPPA reconstruction with R=3 for a single channel covering a sagittal slice through the occipital lobe.
  • Influence of equipment changes on a longitudinal trial
    Ken Sakaie1, Janel Fedler2, Jon Yankey2, Kunio Nakamura1, Josef Debbins3, Mark J. Lowe1, Paola Raska1, and Robert J. Fox1
    1The Cleveland Clinic, Cleveland, OH, United States, 2University of Iowa, Iowa City, IA, United States, 3Barrow Neurological Institute, Phoenix, AZ, United States
    We find that accounting for the specific type of hardware change is optimal, but trends differ between different imaging measures (BPF and TD). We expect these results to be useful in planning future clinical trials.
    Table 2. Effect of adjusting for hardware changes. BPF is unitless, TD units are 10-3 mm2/sec. Values are change per year, with 95% CI’s in parentheses. Models are: no adjustment for hardware changes (Original), data acquired after a hardware change excluded (Exclude), hardware change treated as a binary yes/no time-dependent covariate (Binary) and type of hardware change is a time-dependent covariate (Type). Lower AIC indicates a better fit, with a change of 2 being substantial.
    Table 4. Impact of scanner change on overall outcomes. Values are differences in treatment effect for BPF or TD by a particular type of hardware change versus no change. The GE HDxt to GE MR750 change could not be analyzed because only one subject was affected.
  • Rapid Whole-Brain Myelin Mapping via Selective Inversion Recovery and Compressed SENSE
    Ping Wang1, Nicholas Sisco1, and Richard Dortch1
    1Barrow Neurological Institute, Phoenix, AZ, United States
    For quantitative magnetization transfer using selective inversion recovery, the combination of compressed sensing and parallel imaging yields significantly reduced scan times with little effect on the precision and accuracy of myelin quantification in white matter. 
    Fig. 3: (A) - (C): PSR maps of a healthy volunteer (male, 36 yrs) under the CS-SENSE factors of 3, 8, and 12, respectively; (D) - (F) show the corresponding R1f maps.
    Fig. 1: (A) - (C) PSR maps of the BSA phantoms under the CS-SENSE factor of 3, 8, and 12 respectively, the phantoms' concentrations are labeled on (A). The corresponding R1f maps under each CE-SENSE factor are shown in (D) - (F).
  • Reproducibility and Multi-vendor Accuracy Comparison of T1- and T2- Mapping Using 2D and 3D Synthetic MRI with GRAPPA, SENSE and Compressed Sense
    Maarten Naeyaert1, Tim Vanderhasselt1, Marcel Warntjes2, and Hubert Raeymaekers1
    1Radiology, Universitair Ziekenhuis Brussel, Brussels, Belgium, 2SyntheticMR AB, Linköping, Sweden
    Comparison of 2D and 3D simultaneously acquired T1, T2 and PD maps on 3T scanners of different vendors shows almost-general agreement between scanners, sequence and with and without GRAPPA, SENSE, or compressed sensing. T1 and T2 values approach the reference values and are reproducible.
    Figure 3: Measured versus reference T1 (left) and T2 (right). Colour of the symbols indicates vendor, the shape indicates acceleration used, and filled vs hollow indicates 2D vs 3D respectively.
    Figure 4: Difference between measured and reference T1 (left) and T2 (right) values. Colour of the symbols indicates vendor, the shape indicates acceleration used, and filled vs hollow indicates 2D vs 3D respectively.
  • Improving T1-weighted MRI brain images by optimizing differential identifiability
    Bradley Fitzgerald1, Kausar Abbas2, Thomas M. Talavage1,3, and Joaquin Goni2
    1Electrical & Computer Engineering, Purdue University, West Lafayette, IN, United States, 2Industrial Engineering, Purdue University, West Lafayette, IN, United States, 3Biomedical Engineering, University of Cincinnati, Cincinnati, OH, United States
    Differential identifiability of scan-rescan subject images can be maximized via principal component reconstruction of a T1 anatomical brain MRI dataset. This results in increased similarity between same subject repeated scans, as well as reducing apparent intersession noise in images.
    Figure 1. Differential identifiability and distance metrics for original T1 data and for T1 data with added noise. (A) displays computed Idiff values for original T1 images (black) and T1 images with added noise (red), reconstructed with the full range of possible principle components (PCs). (B) displays Dself (blue) and Dothers (green) values for original T1 images and T1 images with added noise, reconstructed with varying number of PCs.
    Figure 3. T1 anatomical images associated with original and added noise datasets at full (36 PC) reconstruction and optimized Idiff (18 PC) reconstruction.
  • GAN-based analysis for investigation of disease specific image pattern in SWI data of patients suffering from multiple sclerosis
    Alina Lopatina1,2, Stefan Ropele3, Renat Sibgatulin1, Jürgen R Reichenbach1,2,4, and Daniel Güllmar1
    1Medical Physics Group / IDIR, Jena University Hospital, Jena, Germany, 2Michael-Stifel-Center for Data-Driven and Simulation Science, Jena, Germany, 3Department of Neurology, Medical University of Graz, Graz, Austria, 4Center of Medical Optics and Photonics Jena, Jena, Germany
    In this study, we proposed a method to translate susceptibility-weighted images of multiple sclerosis patients to healthy using generative adversarial networks. The method identified voxel information around the central veins and ventricles as MS-specific.
    Figure 1. Image-to-image translation using FPG for two MS subjects for three slice positions (in rows). The input images (MS) are translated to the healthy domain (HC’). The absolute difference between the input and the translated images reveals the disease-related voxels, which are overlaid on the original input images. Green indicates an increase and pink a decrease in intensity.
    Figure 2. Image-to-image translation using FPG on two MS subjects for one slice position. The input images (MS) are translated to the healthy domain (HC’) (upper row) and to the diseased domain (MS’, lower row). The absolute difference between the input and the translated images is threshold-limited to highlight disease-related voxels, which are overlaid on the input images. Green indicates an increase and pink a decrease in intensity.
  • 3D Quantitative MRI of the Brain: Effects of B1 Inhomogeneity in 3D-QALAS
    Anders Tisell1,2, Peter Lundberg1,2, Marcel Jan Bertus Warntjes1, and Frederik Testud3
    1CMIV, Linköping University, Linköping, Sweden, 2Medical Radiation Physics, Linköping University, Linköping, Sweden, 3Siemens Healthcare AB, Malmö, Sweden
    In this work we validate an implementation of 3D QALAS for 3T.  In conclusion we see a high accuracy for estimation of R1 and R2 compared to gold standard measurements. However we also see a small effect of B1+ inhomogeneities in the estimated R1.  
    Figure 1. (top panel) the signal model for the QALAS sequence is shown with M0=1, describing the effect of different flip angles (different colours), and the corresponding measured data. (middle panel) QALAS data from an MS patient showing the signal maps for the individual dynamics. (bottom row) Corresponding QALAS data of a relaxation phantom.
    Figure 3. Bland-Altman plots of the differences between R1 and R2-measurements using the implemented QALAS sequence, compared to gold standard methods. QALAS measurements with FLASH flip angle of 4° (yellow) and green for FLASH flip angle of 10° (green).
Back to Top
Digital Poster Session - Quantitative Neuro: Translational Studies
Neuro
Monday, 17 May 2021 15:00 - 16:00
  • Hippocampal segmentation from 7T images showed reduced subfields volume in Sickle Cell Disease subjects
    Tales Santini1, Minseok Koo1, Nadim Farhat1, Vinicius P. Campos2, Salem Alkhateeb1, Marcelo A. C. Vieira2, Meryl A Butters1, Caterina Rosano1, Howard J Aizenstein1, Joseph Mettenburg1, Enrico M. Novelli1, and Tamer S Ibrahim1
    1University of Pittsburgh, Pittsburgh, PA, United States, 2University of Sao Paulo, Sao Carlos, Brazil
    We collected and processed 7T MRI images from individuals with sickle cell disease and matched controls. Individuals with SCD have significantly smaller volumes of the DG+CA2+CA3 hippocampal region. Other subregions also showed a trend towards smaller volumes
    Neuroimaging characteristics. Mean volume (in mm3), standard deviation and percent difference.
    Example of hippocampal subfields segmentation in a subject with SCD. a: coronal slice of the T2-weighted image, acquired at 7T with resolution 0.375x0.375x1.5mm2; b,c: zoom-in images showing details of the hippocampus structure, subject right and left, respectively; d,e: hippocampus subfield segmentations overlaying the T2-weighted image, subject right and left, respectively; f: 3D reconstruction of the hippocampal subfield segmentations.
  • Altered intrinsic brain functional network dynamics in patients with end-stage renal disease undergoing maintenance hemodialysis
    Baolin Wu1, Feifei Zhang1, Zhiyun Jia1,2, and Qiyong Gong1,3
    1Huaxi MR Research Center (HMRRC), Department of Radiology, West China Hospital of Sichuan University, Chengdu, China, 2Department of Nuclear Medicine, West China Hospital of Sichuan University, Chengdu, China, 3Psychoradiology Research Unit of Chinese Academy of Medical Sciences (2018RU011), Chengdu, China
    Our study demonstrated altered dynamic functional connectivity properties in patients with end-stage renal disease, and the number of transitions were correlated with cognitive performance in those patients.  
    Figure 3. Results of the k-means clustering analysis per state. A. Cluster centroids for each state. B. The strongest 5% of the functional connectivity matrix in each state. Red lines represent positive functional connectivity, and blue lines represent negative functional connectivity. VN = visual network; SMN = sensorimotor network; DMN = default mode network; AN = auditory network; LFPN = left frontoparietal network; RFPN = right frontoparietal network; CER = cerebellum; DAN = dorsal attention network; ECN = executive control network.
    Figure 5. Between-group comparison in temporal properties of functional connectivity states and the partial correlation analysis results. Dot plots show individual data points (squares), averages (transverse lines), and standard deviations (vertical lines) of the mean fractional windows (A), mean dwell time (B), and total number of transitions (C). The transverse lines represent averages, and the vertical lines represent standard deviations. D-E. Partial correlation analysis results. *p < 0.05. TMT-A = Trail Making Test A; SDMT = Symbol Digit Modalities Test
  • Diagnostic accuracy of ASL in comparison with DSC perfusion in the surveillance of different types of brain tumors
    Anna Lavrova1, Wouter Teunissen2, Esther Warnert2, Martin van den Bent3, Vladimir Cheremisin1, and Marion Smits2
    1Radiology, Saint Petersburg University, Saint Petersburg, Russian Federation, 2Radiology, Erasmus MC, Rotterdam, Netherlands, 3Neurology, Erasmus MC, Rotterdam, Netherlands

    Our findings so far suggest using ASL instead of DSC in glioma at 3T, whereas lymphoma and metastases require further investigation in the follow-up analysis.

    Figure 1. Pearson’s correlation coefficients of ASL CBF ratio and uncorrected DSC rCBV in enhancing glioma.
    Figure 2. Pearson’s correlation coefficients of ASL CBF ratio and uncorrected DSC rCBV in non-enhancing glioma.
  • Regional Brain Perfusion Changes in Cognition and Mood Regulatory Sites in Patients with Type 2 Diabetes Mellitus
    Bhaswati Roy1, Sarah Choi2, Matthew J. Freeby 3, and Rajesh Kumar1,4,5,6
    1Anesthesiology, University of California Los Angeles, Los Angeles, CA, United States, 2School of Nursing, University of California Los Angeles, Los Angeles, CA, United States, 3Medicine, Endocrinology - Diabetes and Metabolism, University of California Los Angeles, Los Angeles, CA, United States, 4Bioengineering, University of California Los Angeles, Los Angeles, CA, United States, 5Radiological Sciences, University of California Los Angeles, Los Angeles, CA, United States, 6Brain Research Institute, University of California Los Angeles, Los Angeles, CA, United States
    T2DM patients show cognitive and mood changes, and brain tissue injury in those regions; however, the underlying cause of tissue injury is unknown. We found significantly reduced CBF in multiple brain areas and its associations with functional deficits in those control sites in the condition. 
    Figure 1: Brain regions with reduced cerebral blood flow in T2DM patients over control subjects. The sites with reduced cerebral blood flow in T2DM patients were included the bilateral prefrontal cortices (a, b), right insula (c), bilateral cingulate (d, h), bilateral cerebellum (e, g), and right thalamus (f). All images are in neurological convention (L = left; R = right). Color bar indicates t-statistic values.
    Figure 3: Cognition showed positive associations with cerebral blood flow in T2DM patients. Positive correlations appeared between MoCA scores and cerebral blood flow of the bilateral prefrontal cortices (a, c), right thalamus (b), right putamen (d), bilateral hippocampus (e, i), bilateral cerebellum (f, g), right amygdala (h), and bilateral basal forebrain (j, k). Figure conventions are same as in Figure 1.
  • Application of quantitative susceptibility mapping in assessment of iron content in brain regions of normal children
    Shilong Tang1 and Lisha Nie2
    1Children's Hospital of Chongqing Medical University, Chongqing, China, 2GE Healthcare, MR Research China, Beijing, Beijing, China
    Quantitative susceptibility mapping can evaluate the iron content in each brain region of normal children and facilitate the diagnosis and treatment of clinical diseases.
    Fig. 1 Schematic diagram of the ROI measurement in various regions of the brain
    Fig. 2 Bar chart comparison of the iron content measurement results of various brain regions in each age group
  • Regional Brain Growth in Fetuses with Congenital Diaphragmatic Hernia
    Fedel Machado-Rivas1,2, Lina Acosta Buitrago3, Jungwhan J Choi1,2, Onur Afacan1,2, Clemente Velasco-Annis1, Simon K Warfield1,2, Ali Gholipour1,2, and Camilo Jaimes1,2
    1Radiology, Boston Children's Hospital, Boston, MA, United States, 2Radiology, Harvard Medical School, Boston, MA, United States, 3Universidad del Rosario, Bogota, Colombia
    Brain growth trajectories of fetal subjects with congenital diaphragmatic hernia (CDH) show global and localized volume loss when compared to typically developing fetuses. Volume differences between controls and CDH subjects are associated with hernia morphometric characteristics.
    Axial (A) and corononal (B) heat map of segments differences between controls and CDH subjects. Percent change for CDH subjects when compared to controls is displayed. Segments in blue represent volume gain, while segments in red correspond to volume loss. White segments were not significantly different from controls.
    (A)T2-weighted super-resolution volume reconstruction with inter-slice motion correction, and intensity normalization of a T2-weighted fetal acquisition (29 week-old fetus). (B) Propagation of fetal atlas labels. (C) 3D volumetric rendering of propagated labels.
  • Inter-rater reliability of sciatic nerve evaluation with MR neurography: a comparison study of multiple sequences
    Ryuna Kurosawa1, Hajime Yokota2, Takafumi Yoda1, Takayuki Sada1, Koji Matsumoto1, Takashi Namiki3, Masami Yoneyama3, Yoshitada Masuda1, and Takashi Uno2
    1Department of Radiology, Chiba University Hospital, Chiba, Japan, 2Diagnostic Radiology and Radiation Oncology, Graduate School of Medicine, Chiba University, Chiba, Japan, 3Philips Japan, Tokyo, Japan
    To evaluate inter-rater reliability of sciatic nerve evaluation with various MRN sequences. 3D-iNerveVIEW had high inter-rater reliability for signal measurements. For MRN measurements, it is important to understand the characteristics of each sequence.
    Figure 4. Intraclass coefficient correlations between two examiners for the sciatic nerve size and signal measurements. Bold represents ICC> 0.601. T2-FFE, DN, and DTI show substantial agreements in the long diameter, while T2WI, 3D-iNerveVIEW, and DTI show substantial agreements in the short diameter. Regarding the signal measurements, 3D-iNerveVIEW shows substantial to almost perfect agreements in signal-to-noise ratio (SNR), contrast-to-noise ratio (CNR), and contrast ratio (CR).
    Figure 5. Problems with 3D-iNerveVIEW. Case1: By suppressing the signal of the vessel, the boundary between nerves and blood vessels becomes unclear. In addition, the vessel walls running in parallel in the long diameter direction of the sciatic nerve become high signals, which reduces visibility. Case2: The vessels are misidentified and measured as nerves when the signal of the vessels running parallel to the long diameter of the sciatic nerve cannot be suppressed.
  • White Matter fiber orientation dependent R2*: comparison between post mortem in situ and in vivo
    Celine Berger1,2, Claudia Lenz1,2, Melanie Bauer1,2, Eva Scheurer1,2, and Christoph Birkl3
    1Institute of Forensic Medicine, Department of Biomedical Engineering, Basel, Switzerland, 2Institute of Forensic Medicine, Health Department Basel-Stadt, Basel, Switzerland, 3Department of Neuroradiology, Medical University of Innsbruck, Austria
    R2* increased with increasing fiber angle in both post mortem in situ and in vivo, whereby a decreased R2* orientation dependency was observed post mortem compared to in vivo.
    Figure 2: R2* as a function of the WM fiber angle averaged over the in vivo subjects (red) and the post mortem subjects (blue) fitted with the absolute values (A) and normalized to the global mean WM R2* (B). The shaded areas represent the 95% CI of the measured data.
    Figure 1: Representative maps of R2* (A) and fiber angle θ (B) of two post mortem subjects with brain temperatures of 5.6°C (top row) and 14.9 °C (middle row) and one in vivo subject shown in the bottom row.
  • Reproducibility of brain ultrashort-T2* component measurements in healthy volunteers
    Nikhil Deveshwar1,2 and Peder E. Z. Larson1,2
    1Department of Radiology and Biomedical Imaging, University of California, San Francisco, San Francisco, CA, United States, 2UC Berkeley - UCSF Graduate Program in Bioengineering, Berkeley and San Francisco, CA, United States
    Ultrashort-T2* fraction parameter maps generated at different scans show similar looking structures in the same volunteer. The distributions of the ultrashort-T2* fractional component in various brain ROIs show similar distributions suggesting this technique is reproducible.
    Figure 1: The proposed image processing pipeline. UTE scans at three different flip angles (18, 12, and 6 degrees) are fit with the following signal model. The resulting parameter map corresponding to the ultrashort-T2* fractional component is then used to isolate values in various brain ROIs.
    Figure 3: Split violin plots comparing the distributions of the ultrashort-T2* component fraction in various brain ROIs across 4 healthy volunteers. The dashed lines represent the interquartile range of each distribution.
  • Correlating Concussion-Related Symptoms to the Personalized MRI Assessment of Brain Abnormalities in Children
    Ethan Danielli1,2, David Stillo1,2, Rachelle Ho3,4, Carol DeMatteo3,5, Geoffrey B Hall4, Nicholas A Bock4, John F Connolly1,4,5,6, and Michael D Noseworthy1,2,5,7,8
    1School of Biomedical Engineering, McMaster University, Hamilton, ON, Canada, 2Imaging Research Centre, St. Joseph's Healthcare Hamilton, Hamilton, ON, Canada, 3School of Rehabilitation Sciences, McMaster University, Hamilton, ON, Canada, 4Department of Psychology, Neuroscience & Behaviour, McMaster University, Hamilton, ON, Canada, 5ARiEAL Research Centre, McMaster University, Hamilton, ON, Canada, 6Department of Linguistics, McMaster University, Hamilton, ON, Canada, 7Department of Electrical and Computer Engineering, McMaster University, Hamilton, ON, Canada, 8Department of Radiology, McMaster University, Hamilton, ON, Canada
    Decreased FA was significantly correlated with younger age and higher PCSS scores. Increased RD was significantly correlated with younger age and the interaction between time to scan and Post-Concussion Symptom Scale (PCSS) score.
    Figure 1. Plots of significant findings between DTI metrics and demographic information. Top left: Total FA injury burden versus Age; Top Right: Total FA injury burden versus PCSS Score; Bottom Left: Total RD injury burden versus Age; Bottom Right: A paired matrix plot indicating the distribution of each metric and the correlation between metrics.
    Table 1. Results for the multilinear regression of demographic and DTI metrics.
  • Cerebral perfusion network analysis to understand cognition in old age: a principal component analysis of ASL-MRI
    Jodi Karlyn Watt1,2,3, Stefan Pszczolkowski1,2,3, Yue Xing1,2,3, Christopher Tench1,2,3, Dorothee Auer1,2,3, and Alzheimer's Disease Neuroimaging Initiative4
    1Division of Clinical Neuroscience, University of Nottingham, Nottingham, United Kingdom, 2Sir Peter Mansfield Imaging Centre, University of Nottingham, Nottingham, United Kingdom, 3NIHR Nottingham Biomedical Research Centre, University of Nottingham, Nottingham, United Kingdom, 4Alzheimer's Disease Neuroimaging Initiative, Los Angeles, CA, United States
    The relationship between CBF and tests of global cognition remains unknown. Principal Component Analysis was used to investigate the MoCA/cerebral perfusion relationship, providing preliminary support for a perfusion pattern which partially explains old-age cognitive decline.
    Figure 2: Combined surviving components (components 1, 4 and 42, 98% of r2 retained) overlaid on the MNI152, 2mm brain template. Negative associations with MoCA are depicted in blue, and positive in red, with the ECN map11 in yellow.
    Figure 4: Colour-coded graph depicting the interrelatedness between any pairs of nodes. Positive associations are depicted with red edges, and negative with blue. The circle represents the mid-axial slice of the brain.
  • Altered structural connectivity and impairment of brain network-cognition relationship in obstructive sleep apnea (OSA)
    Tasfiya Islam1, Mengting Liu1, Dae Lim Koo2, Ryan Cabeen1, Eunyeon Joo3, and Hosung Kim1
    1USC Stevens Neuroimaging and Informatics Institute, Keck School of Medicine of University of Southern California, Los Angeles, CA, United States, 2Department of Neurology, Boramae Medical Center, Seoul National University College of Medicine, Seoul, Korea, Republic of, 33Department of Neurology, Samsung Medical Center, Sungkyunkwan University School of Medicine, Samsung Biomedical Research Institute, Seoul, Korea, Republic of
    OSA patients were found to have globally less efficient brain network. In OSA, the higher modularity linked with network segregation and less efficient network may further lead to the failure in coordinating cognitive resources and the impaired visual memory function.
    Figure 2: Network based statistics (NBS) under the contrast [-1, 1] for OSA and healthy control group connectivity comparison.
    Figure 3: Scatter plot of significant effects found in data collected from behavioral versus brain measure comparison for both healthy and OSA subjects under threshold value of 0.35.
  • QSM detects early alterations of brain venous blood oxygenation in fetuses with complex congenital heart diseases
    Cong Sun1, Aocai Yang1, Jiaguang Song2, Minhui Ouyang3, Jinxia Zhu4, Lei Xue5, Hao Huang3,6, and Gunagbin Wang1
    1Radiology, Shandong Medical Imaging Research Institute, Cheeloo College of Medicine, Shandong University, Jinan, China, 2Ultrasound, Shandong Provincial Hospital Affiliated to Shandong University, Jinan, China, 3Radiology Research, Children’s Hospital of Philadelphia, Philadelphia, PA, United States, 4MR Collaboration, Healthcare Siemens Ltd., Beijing, China, 5MR Application, Siemens Healthineers Ltd., Jinan, China, 6Radiology, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
    Fetuses with various kinds of complex congenital heart disease had higher values of venous blood oxygen saturation than did normal fetuses of the same gestational age, as measured with quantitative susceptibility mapping in utero. 
    Figure 4a: The venous blood oxygen saturation (SvO2) alterations across gestational age (GA) of the normal group: SvO2 = -0.004223*GA + 0.8783. Figure 4b: SvO2 alterations across gestational age between the congenital heart disease (CHD) group (SvO2 = -0.00514*GA + 0.942) and GA-matched normal group (SvO2 = -0.01205*GA + 1.087) and their corresponding fitting curve. Red circles and red curves represent the CHD group; blue circles and blue curves represent the normal group.
    Figure 5: Venous blood oxygen saturation (SvO2) values in the congenital heart disease (CHD) group and gestational age (GA)-matched normal group, along with their corresponding mean and 95% CI SvO2 values. In the analysis of covariance, the SvO2 value showed a significant difference between the CHD (80.8%±4.6%) group and the GA-matched normal group (75.7%±8.0%) after the effects of GA were excluded (p=0.038).
  • Quantitative cerebral oxygenation mapping by MRI with whole brain coverage compared to PET
    Jan Kufer1, Christine Preibisch1,2, Samira Epp1, Jens Goettler1,3, Kilian Weiss4, Mikkel Bo Hansen5, Claus Zimmer1, Kim Mouridsen5, Fahmeed Hyder3, and Stephan Kaczmarz1,3
    1Department of Neuroradiology, School of Medicine, Technical University of Munich (TUM), Munich, Germany, 2Department of Neurology, School of Medicine, Technical University of Munich (TUM), Munich, Germany, 3Department of Radiology & Biomedical Imaging (MRRC), Yale University, New Haven, CT, United States, 4Philips Healthcare, Hamburg, Germany, 5Department of Clinical Medicine, Aarhus University, Aarhus, Denmark
    Comparison of two MRI-based cerebral oxygenation biomarkers OEF from mqBOLD and OEC modelled from DSC showed good agreement with PET reference data in cortical gray matter in young healthy volunteers. Some variation between MRI-based OEF and OEC was attributed to different underlying models.
    Figure 3: Spatial correlation across 28 gray matter VOIs. Each cross indicates the mean parameter value in a specific Brodmann area averaged across subjects. Strongest spatial correlation was found between PET- and mqBOLD-OEF (A) with r=0.69 (p<0.05). Correlation between PET-OEF and DSC-OEC (B) was also strong, but slightly weaker with r=0.50 (p<0.05). Similarly, good correlation was also found between both MRI techniques (C) with r=0.57 (p<0.05). Note that MRI and PET data were acquired in similar cohorts, but different subjects.
    Figure 2: Exemplary PET-OEF, MRI-OEF and MRI-OEC data in four different slices. PET-OEF (left column) is quite homogenous across the entire brain, without contrast between gray (GM) vs. white matter (WM). MRI-based OEF (central column) and OEC (right column) show some similarity with PET. However, both parameter maps show an artifactual GM/WM contrast with elevated WM values. While OEF values in GM seem similar between PET and mqBOLD, OEC values based on the CFIN-model are generally lower.
  • Carbogen-based Cerebrovascular Reserve Using BOLD-based fMRI of Human Brain : tissue, territorial and cortical specificity
    Tzu-chen Yeh1,2, Chou-ming Cheng3, and Chi-che Chou3
    1Department of Radiology, Taipei Veterans General Hospital, Taipei, Taiwan, 2Institue of Brain Science, National Yang-Ming University, Taipei, Taiwan, 3Department of Medical Research, Taipei Veterans General Hospital, Taipei, Taiwan
    A custom-designed system was applied for inhalation of various gas mixtures and air. Independent component analysis was applied to obtain the carbogen response function for mapping of CO2-CVR. And we demonstrated the tissue, territorial and cortical specificities of CO2-CVRof human brain.
    The custom designed automatic delivery system of CO2-CVR included mass flow controllers with medical graded gas mixtures and air, LabView for device control, CO2 analyzer with PoweLab and Chart 5 with synchronization with radiofrequency pulse of MRI.
    Cortical specificity of CO2-CVR showed highest t values of one-sample t test using carbon response function (see text for details) as 14.53 at ventral portion of Brodmann area 23 (blue). And cortices with CO2-CVR larger than 13 were labelled as green and red for 13 < t values < 14 and t values >14, respectively.
  • Clinical Whole-Brain R2* and Quantitative Susceptibility Maps at 3T – Reproducibility and Parameter Optimization Towards Millimetric 5min Scan
    Thomas Troalen1, Arnaud Le Troteur2, Sylviane Confort-Gouny2, Patrick Vioux2, Claire Costes2, Lauriane Pini2, Jean-Philippe Ranjeva2, Maxime Guye2, and Ludovic de Rochefort2
    1Siemens Healthcare SAS, Saint-Denis, France, 2CRMBM UMR7339 CNRS Aix-Marseille Université, Marseille, France
    This work demonstrates the ability to accelerate multi gradient echo sequences for a joint R2* and QSM in the brain. Combining this sequence with a state-of-the art automatic post-processing pipeline, we propose here a standardized whole brain clinical protocol of 5 minutes.
    Figure 1: Automatic processing pipeline: The ICBM 2009a nonlinear symmetric MNI template was non-linearly registered to the subject’s T1w space (using Vida_64Ch as reference). Inter/Intra-run co-registrations were performed (T1wRef-to-T1w, as well as T1w-to-R2*). Tissue segmentation was achieved using SPM12 software using the default brain probability maps. Labels were propagated to R2*/QSM space and restricted to WM/GM tissue types. WM/GM lobes and cerebellum were extracted, as well as four deep grey nuclei and the brainstem prior ROI analysis on the quantitative maps.
    Reproducibility across MR system types (Verio and Vida), receive head coils (20, 32 and 64 channels) and sequence parameters. Note that the measured variability on the Vida with different head coils was not different than the intra-run metric variations on the Verio system. The signal drop caused by the increased acceleration factor from 3 to 4 was not reflected in the R2* and QSM estimations, with values in the same range as the other protocols. The central table reports mean and standard deviation per segmented region across all measurements.
  • Interhemispheric Functional Connectivity and 18F-fluatmetamol PET differentiate AD, amyloid and non-amyloid Mild cognitive decline
    Eva YW Cheung1, Patrick KC Chiu2, YF Shea2, Joseph SK Kwan3, and Henry KF Mak1
    1The University of Hong Kong, Hong Kong, Hong Kong, 2Queen Mary Hospital, Hong Kong, Hong Kong, 3Imperial College London, London, United Kingdom
    Interhemispheric functional connectivity and 18F-flutametamol PET 
    AD
    amyMCI
  • Associations of Musical Aptitude with High Angular Resolution Diffusion Imaging (HARDI) derived structural connectivity
    Archith Rajan1, Apurva Shah1, Madhura Ingalhalikar1, and Nandini C Singh2
    1Symbiosis Centre for Medical Image Analysis, Symbiosis International(Deemed) University, Pune, India, 2Language,Literacy and Music Lab, National Brain Research Centre, Gurgaon, India
    • Prevalence of inter-hemispheric than intra-hemispheric structural connectivity that showed positive linear associations with better music perception abilities.
    • Sequential and not sensory sub-scores showed associations with whole brain structural connectivity

    Fig.1 Significant edge list for positive associations with d’-total scores. The nodes with its total number of connections are also indicated.
    Fig.2 Significant edge list for positive associations with d’s computed from sequential processing component (comprising Melody, Standard Rhythm, Embedded Rhythm and Accent). The nodes with its total number of connections are also indicated.
  • Imaging intracortical structure using navigator-based, motion and B0-corrected T2*-weighted MRI at 7 T
    Jiaen Liu1, Peter van Gelderen1, Jacco A. de Zwart1, and Jeff H. Duyn1
    1AMRI, LFMI, NINDS, National Institutes of Health, Bethesda, MD, United States
    High-resolution T2*-weighted 7 T MRI can be used to delineate intracortical structure owning to its high signal-to-noise and contrast-to-noise ratio. Reliable ultrahigh resolution T2*-weighted MRI can be reliably obtained using a navigator-based sequence with motion and B0 correction.
    Fig. 1 Top row: magnitude of corrected, uncorrected and lower-resolution corrected (isotropic 0.6 mm) T2*-weighted MRI from one subject. The yellow arrow points to the hypointensive Line of Gennari in the primary visual cortex. Bottom row: the corresponding susceptibility-induced off-resonance frequency maps at 7 T. The Line of Gennari appears as a positive-frequency band (note darker equals more positive).
    Fig. 2 Normalized root mean square error (NRMSE) between two repeated T2*-weighted scans for each subject. Red line represents the expected NRMSE caused by thermal noise.
  • Detecting normal Fetal brain development with T1Mapping Imaging Technique
    Yan-Chao Liu1, Bo-Hao Zhang2, De-Sheng Xuan2, Xue-Yuan Wang2, Kai-Yu Wang3, Xin Zhao2, and Xiao-An Zhang2
    1Department of Radiology, The Third Affiliated Hospital of Zhengzhou University, Zhengzhou, China, 2Department of Radiology, the Third Affiliated Hospital of Zhengzhou University, Zhengzhou, 450052, China, Zhengzhou, China, 3MR Research China, GE Healthcare, Beijing 100000, PR China, Beijing, China
    Multiple MRI quantitative techniques have confirmed that postnatal brain development is an ongoing maturation process. Fetal brain development is an ongoing process.  In this work, T1Mapping allowed quantitative assessment of fetal brain development.
    T1Mapping maps in the fetal brain (A–D). Examples of ROIs for DTI are shown in white color (1–2). The area of the ROI was adjusted appropriately according to the gestational week and anatomical structures. Regions of interest: 1, Thalamus; 2, Corticospinal fibers. A and B are T1Mapping images of the same fetus; And gestation ages was 28 week. C and D are T1Mapping images of the same fetus; And gestation ages was 36 week.
    T1 relaxation time trajectories in Thalamus for 17 subjects. T1 had good correlation with gestational age. A regression line is shown along with 95% confidence intervals (dotted lines).