Large angle spin-echo imaging

Mouse Brain Extractor: Brain segmentation of mouse MRI using global positional encoding and SwinUNETR

In spite of the great progress that has been made towards automating brain extraction in human magnetic resonance imaging (MRI), challenges remain in the automation of this task for mouse models of brain disorders. Researchers often resort to editing brain segmentation results manually when automated methods fail to produce accurate delineations. However, manual corrections can be labor-intensive and introduce interrater variability. This motivated our development of a new deep-learning-based method for brain segmentation of mouse MRI, which we call Mouse Brain Extractor. We adapted the existing SwinUNETR architecture (Hatamizadeh et al., 2021) with the goal of making it more robust to scale variance. Our approach is to supply the network model with supplementary spatial information in the form of absolute positional encoding. We use a new scheme for positional encoding, which we call Global Positional Encoding (GPE). GPE is based on a shared coordinate frame that is relative to the entire input image. This differs from the positional encoding used in SwinUNETR, which solely employs relative pairwise image patch positions. GPE also differs from the conventional absolute positional encoding approach, which encodes position relative to a subimage rather than the entire image. We trained and tested our method on a heterogeneous dataset of N=223 mouse MRI, for which we generated a corresponding set of manually-edited brain masks. These data were acquired previously in other studies using several different scanners and imaging protocols and included in vivo and ex vivo images of mice with heterogeneous brain structure due to different genotypes, strains, diseases, ages, and sexes. We evaluated our method's results against those of seven existing rodent brain extraction methods and two state-of-the art deep-learning approaches, nnU-Net (Isensee et al., 2018) and SwinUNETR. Overall, our proposed method achieved average Dice scores on the order of 0.98 and average HD95 measures on the order of 100 μm when compared to the manually-labeled brain masks. In statistical analyses, our method significantly outperformed the conventional approaches and performed as well as or significantly better than the nnU-Net and SwinUNETR methods. These results suggest that Global Positional Encoding provides additional contextual information that enables our Mouse Brain Extractor to perform competitively on datasets containing multiple resolutions.

Alpha-180 spin-echo based line-scanning method for high resolution laminar-specific fMRI

Laminar-specific functional magnetic resonance imaging (fMRI) has been widely used to study circuit-specific neuronal activity by mapping spatiotemporal fMRI response patterns across cortical layers. Hemodynamic responses reflect indirect neuronal activity given limit of spatial and temporal resolution. Previous gradient-echo based line-scanning fMRI (GELINE) method was proposed with high temporal (50 ms) and spatial (50 µm) resolution to better characterize the fMRI onset time across cortical layers by employing 2 saturation RF pulses. However, the imperfect RF saturation performance led to poor boundary definition of the reduced region of interest (ROI) and aliasing problems outside of the ROI. Here, we propose α (alpha)-180 spin-echo-based line-scanning fMRI (SELINE) method to resolve this issue by employing a refocusing 180° RF pulse perpendicular to the excitation slice. In contrast to GELINE signals peaked at the superficial layer, we detected varied peaks of laminar-specific BOLD signals across deeper cortical layers with the SELINE method, indicating the well-defined exclusion of the large drain-vein effect with the spin-echo sequence. Furthermore, we applied the SELINE method with 200 ms TR to sample the fast hemodynamic changes across cortical layers with a less draining vein effect. In summary, this SELINE method provides a novel acquisition scheme to identify microvascular-sensitive laminar-specific BOLD responses across cortical depth.

In vivo irreversible and reversible transverse relaxation rates in human cerebral cortex via line scans at 7 T with 250 micron resolution perpendicular to the cortical surface

Understanding how and why MR signals and their associated relaxation rates vary with cortical depth could ultimately enable the noninvasive investigation of the laminar architecture of cerebral cortex in the living human brain. However, cortical gray matter is typically only a few millimeters thick, making it challenging to sample many cortical depths with the voxel sizes commonly used in MRI studies. Line-scan techniques provide a way to overcome this challenge and here we implemented a novel line-scan GESSE pulse sequence that allowed us to measure irreversible and reversible transverse relaxation rates-R₂ and R₂´, respectively-with extremely high resolution (250 μm) in the radial direction, perpendicular to the cortical surface. Eight healthy human subjects were scanned at 7 T using this sequence, with primary visual cortex (V1) targeted in three subjects and primary motor (M1) and somatosensory cortex (S1) targeted in the other five. In all three cortical areas, a peak in R₂ values near the central depths was seen consistently across subjects-an observation that has not been made before, to our knowledge. On the other hand, no consistent pattern was apparent for R₂´ values as a function of cortical depth. The intracortical R₂ peak reported here is unlikely to be explained by myelin content or by deoxyhemoglobin in the microvasculature; however, this peak is in accord with the laminar distribution of non-heme iron in these cortical areas, known from prior histology studies. Obtaining information about tissue microstructure via measurements of transverse relaxation (and other quantitative MR contrast mechanisms) at the extremely high radial resolutions achievable through the use of line-scan techniques could therefore bring us closer to being able to perform "in vivo histology" of the cerebral cortex.

Brain orchestration of pregnancy and maternal behavior in mice: A longitudinal morphometric study

Reproduction induces changes within the brain to prepare for gestation and motherhood. However, the dynamic of these central changes and their relationships with the development of maternal behavior remain poorly understood. Here, we describe a longitudinal morphometric neuroimaging study in female mice between pre-gestation and weaning, using new magnetic resonance imaging (MRI) resources comprising a high-resolution brain template, its associated tissue priors (60-µm isotropic resolution) and a corresponding mouse brain atlas (1320 regions of interest). Using these tools, we observed transient hypertrophies not only within key regions controlling gestation and maternal behavior (medial preoptic area, bed nucleus of the stria terminalis), but also in the amygdala, caudate nucleus and hippocampus. Additionally, unlike females exhibiting lower levels of maternal care, highly maternal females developed transient hypertrophies in somatosensory, entorhinal and retrosplenial cortices among other regions. Therefore, coordinated and transient brain modifications associated with maternal performance occurred during gestation and lactation.

Effect of radiofrequency inhomogeneity on water-content based electrical properties tomography and its correction by flip angle maps

Water-content based electrical properties tomography (wEPT) can retrieve electrical properties (EPs) from water-content maps. B₁⁺ field information is not involved in the traditional magnetic resonance electrical properties tomography approach. wEPT can be performed through conventional MR scanning, such as T₁-weighted spin-echo imaging, which provides convenient access to multiple clinical applications. However, the inhomogeneous radiofrequency (RF) field induced by RF coils would cause inaccuracy in wEPT reconstructions during MR scanning. We conducted a detailed investigation to evaluate the effect of inhomogeneous RF field on wEPT reconstructions to guarantee that EP mapping is desired for clinical practice. Two important considerations are involved, namely, multiple typical coil configurations and various flip angles (FAs). We proposed a correction scheme with actual FA mapping to calibrate the RF inhomogeneity and finally validated it by using human imaging at 3 T. This study illustrates a detailed evaluation for wEPT under imperfect RF homogeneity and further provides a feasible correction procedure to mitigate it. The profound knowledge of wEPT provided in our work will benefit its performance in clinical applications.

In Vivo 3D Digital Atlas Database of the Adult C57BL/6J Mouse Brain by Magnetic Resonance Microscopy

In this study, a 3D digital atlas of the live mouse brain based on magnetic resonance microscopy (MRM) is presented. C57BL/6J adult mouse brains were imaged in vivo on a 9.4 Tesla MR instrument at an isotropic spatial resolution of 100 μm. With sufficient signal-to-noise (SNR) and contrast-to-noise ratio (CNR), 20 brain regions were identified. Several atlases were constructed including 12 individual brain atlases, an average atlas, a probabilistic atlas and average geometrical deformation maps. We also investigated the feasibility of using lower spatial resolution images to improve time efficiency for future morphological phenotyping. All of the new in vivo data were compared to previous published in vitro C57BL/6J mouse brain atlases and the morphological differences were characterized. Our analyses revealed significant volumetric as well as unexpected geometrical differences between the in vivo and in vitro brain groups which in some instances were predictable (e.g. collapsed and smaller ventricles in vitro) but not in other instances. Based on these findings we conclude that although in vitro datasets, compared to in vivo images, offer higher spatial resolutions, superior SNR and CNR, leading to improved image segmentation, in vivo atlases are likely to be an overall better geometric match for in vivo studies, which are necessary for longitudinal examinations of the same animals and for functional brain activation studies. Thus the new in vivo mouse brain atlas dataset presented here is a valuable complement to the current mouse brain atlas collection and will be accessible to the neuroscience community on our public domain mouse brain atlas website.

In Vivo 3D Digital Atlas Database of the Adult C57BL/6J Mouse Brain by Magnetic Resonance Microscopy

In this study, a 3D digital atlas of the live mouse brain based on magnetic resonance microscopy (MRM) is presented. C57BL/6J adult mouse brains were imaged in vivo on a 9.4 Tesla MR instrument at an isotropic spatial resolution of 100 mum. With sufficient signal-to-noise (SNR) and contrast-to-noise ratio (CNR), 20 brain regions were identified. Several atlases were constructed including 12 individual brain atlases, an average atlas, a probabilistic atlas and average geometrical deformation maps. We also investigated the feasibility of using lower spatial resolution images to improve time efficiency for future morphological phenotyping. All of the new in vivo data were compared to previous published in vitro C57BL/6J mouse brain atlases and the morphological differences were characterized. Our analyses revealed significant volumetric as well as unexpected geometrical differences between the in vivo and in vitro brain groups which in some instances were predictable (e.g. collapsed and smaller ventricles in vitro) but not in other instances. Based on these findings we conclude that although in vitro datasets, compared to in vivo images, offer higher spatial resolutions, superior SNR and CNR, leading to improved image segmentation, in vivo atlases are likely to be an overall better geometric match for in vivo studies, which are necessary for longitudinal examinations of the same animals and for functional brain activation studies. Thus the new in vivo mouse brain atlas dataset presented here is a valuable complement to the current mouse brain atlas collection and will be accessible to the neuroscience community on our public domain mouse brain atlas website.

Fully balanced steady-state 3D-spin-echo (bSSSE) imaging at 3 Tesla

In this work a 3D fully balanced steady-state spin-echo sequence (bSSSE) was designed and implemented on a 3-Tesla whole-body MR scanner. In contrast to previously introduced nonbalanced SSSE (nbSSSE) sequences, bSSSE features fully rewound imaging gradients. Additionally, all phase encodings and the readout prephasing are performed after the refocusing pulse. Its performance was compared with previously used gradient-echo (GE) and nbSSSE sequences. The signal response of each sequence was simulated by numerical solution of the Bloch equation. All sequences were then optimized based on these simulations for high-resolution (HR) imaging of tissues with short relaxation times and applied to in vivo HR-MRI of trabecular bone. bSSSE outperformed nbSSSE sequences in terms of SNR and SNR efficiency. Measurements of structural bone parameters revealed thinner trabeculae and a lower bone/bone-marrow fraction for the bSSSE sequence compared to the nbSSSE sequence. By applying a parallel imaging technique based on generalized autocalibrating partially parallel acquisition (GRAPPA) with a reduction factor of 2, we were able to maintain a clinically feasible scan time using bSSSE. Comparisons of structural bone parameters revealed a difference of less than 3% for all structural parameters between parallel and conventional imaging.

Automated segmentation of neuroanatomical structures in multispectral MR microscopy of the mouse brain

We present the automated segmentation of magnetic resonance microscopy (MRM) images of the C57BL/6J mouse brain into 21 neuroanatomical structures, including the ventricular system, corpus callosum, hippocampus, caudate putamen, inferior colliculus, internal capsule, globus pallidus, and substantia nigra. The segmentation algorithm operates on multispectral, three-dimensional (3D) MR data acquired at 90-microm isotropic resolution. Probabilistic information used in the segmentation is extracted from training datasets of T2-weighted, proton density-weighted, and diffusion-weighted acquisitions. Spatial information is employed in the form of prior probabilities of occurrence of a structure at a location (location priors) and the pairwise probabilities between structures (contextual priors). Validation using standard morphometry indices shows good consistency between automatically segmented and manually traced data. Results achieved in the mouse brain are comparable with those achieved in human brain studies using similar techniques. The segmentation algorithm shows excellent potential for routine morphological phenotyping of mouse models.

Measurement and correction of transmitter and receiver induced nonuniformities in vivo

Signal intensity nonuniformities in high field MR imaging limit the ability of MRI to provide quantitative information and can negatively impact diagnostic scan quality. In this paper, a simple method is described for correcting these effects based on in vivo measurement of the transmission field B1+ and reception sensitivity maps. These maps can be obtained in vivo with either gradient echo (GE) or spin echo (SE) imaging sequences, but the SE approach exhibits an advantage over the GE approach for correcting images over a range of flip angles. In a uniform phantom, this approach reduced the ratio of the signal SD to its mean from around 30% before correction to approximately 6% for the SE approach and 9% for the GE approach after correction. The application of the SE approach for correcting intensity nonuniformities is demonstrated in vivo with human brain images obtained using a conventional spin echo sequence at 3.0 T. Furthermore, it is also shown that this in vivo B1+ and reception sensitivity mapping can be performed using segmented echo planar imaging sequences providing acquisition times of less than 2 min. Although the correction presented here is demonstrated with a simultaneous transmit and receive volume coil, it can be extended to the case of separate transmission and reception coils, including surface and phase array coils.

Optimization of Signal Intensity and T1-Dependent Contrast with Nonstandard Flip Angles in Spin-Echo and Inversion-Recovery MR Imaging

Theoretical expressions of flip angle values maximizing signal intensity and T1-dependent contrast are derived for spin-echo and inversion-recovery sequences. Experimental data and theoretical predictions are closely correlated for experiments carried out on phantoms, despite the nonideal shape of the RF refocusing pulse used for slice selection. The use of nonstandard angles is justified when rapid MR acquisitions are needed and/or large T1 species must be imaged with refocused sequences. Copyright 1998 Academic Press. Copyright 1998 Academic Press

Fast 3D large-angle spin-echo imaging (3D FLASE)

A rapid steady-state 3D spin-echo imaging pulse sequence, based on the principle of nutating the spins by an angle greater than 90 degrees, has been designed and implemented on a clinical 1.5-T whole-body MR scanner. The pulse sequence, denoted fast large-angle spin-echo (FLASE), has been optimized for high-resolution imaging of tissues with short T2 and T2*. Features of FLASE include a minimum-phase Shinnar-Le Roux excitation pulse and distribution of phase- and slice-encoding gradients before and after the 180 degrees refocusing pulse to minimize the critical time delay between inversion and restoration of the residual longitudinal magnetization and for minimizing echo time. A Bloch equation analysis, corroborated by experimental data, shows FLASE signal-to-noise to be superior to its closest analog, 3D rapid spin-echo excitation (RASEE) (Jara et al., Magn Reson Medicine 29, 528 (1993)), and 3D gradient-recalled acquisition in steady state (GRASS). It is demonstrated that with judicious RF phase-cycling and steady state operation, FLASE can produce high-quality microimages free of intravoxel phase dispersion from susceptibility-induced background gradients. The performance of the method is exemplified with ultra high-resolution images of trabecular bone in vitro and in vivo in the human calcaneus and wrist at voxel sizes as low as 98 x 98 x 200 microns3. Finally, the contrast behavior of refocused FLASE can be altered by disrupting the steady state analogous to gradient echo imaging.