Three-dimensional texture analyses of multi-quantitative relaxation time maps for evaluating cartilage repair with the treatment of allogeneic human adipose-derived mesenchymal progenitor cells.

BACKGROUND: To explore the ability of three-dimensional texture analyses based on gray-level run-length matrix (GLRLM) for examining the spatial distribution of pixel values on magnetic resonance imaging (MRI) relaxation time maps and detecting the compositional variation of cartilage repair following treatment with allogeneic human adipose-derived mesenchymal progenitor cells (haMPCs). METHODS: Participants with knee osteoarthritis were randomly divided into three groups with intra-articular haMPCs injections: low-, medium-, and high-dose groups. We analyzed five GLRLM parameters in the T1rho, T2 and T2star maps, including run length non-uniformity (RLNonUni), gray-level non-uniformity (GLevNonU), long run emphasis (LngREmph), short run emphasis (ShrtREmp), and fraction of images in runs. We used the relative D values (the ratio of difference values to baseline) as the objective to avoid errors caused by individual differences. We calculated the two-tailed Pearson's linear correlation coefficient (r) to investigate the correlations of the texture parameters with the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) scores. RESULTS: Compared with the base time, significant reduction of WOMAC score was observed in both high and medium doses groups at terminal time, indicating relief of pain symptoms in high and medium groups with the treatment of allogeneic haMPCs. Significant differences were observed in the GLRLM parameters of cartilage MR relaxation time maps in different doses groups. In both T1rho and T2 relaxation time maps, the high-dose group showed significant increases in relative D values of RLNonUni, GLevNonU, LngREmph and ShrtREmp, which indicated significant changes in the uniformity of relaxation time maps. For T2star map, GLRLM parameters such as GLevNonU and ShrtREmp, especially LngREmph, showed significant increases in relative D values in high-dose group. Among all GLRLM features, LngREmph of three relaxation time maps had performed excellent linear correlations with WOMAC scores. CONCLUSIONS: Texture analysis of the cartilage may allow the detection of compositional variation in cartilage repair with the treatment of allogeneic haMPCs. This technique displays potential applications in understanding the mechanism of stem cell repair of the cartilage and assessing the treatment response.

Cell-type-specific optogenetic fMRI on basal forebrain reveals functional network basis of behavioral preference.

The basal forebrain (BF) is a complex structure that plays key roles in regulating various brain functions. However, it remains unclear how cholinergic and non-cholinergic BF neurons modulate large-scale functional networks and their relevance in intrinsic and extrinsic behaviors. With an optimized awake mouse optogenetic fMRI approach, we revealed that optogenetic stimulation of four BF neuron types evoked distinct cell-type-specific whole-brain BOLD activations, which could be attributed to BF-originated low-dimensional structural networks. Additionally, optogenetic activation of VGLUT2, ChAT, and PV neurons in the BF modulated the preference for locomotion, exploration, and grooming, respectively. Furthermore, we uncovered the functional network basis of the above BF-modulated behavioral preference through a decoding model linking the BF-modulated BOLD activation, low-dimensional structural networks, and behavioral preference. To summarize, we decoded the functional network basis of differential behavioral preferences with cell-type-specific optogenetic fMRI on the BF and provided an avenue for investigating mouse behaviors from a whole-brain view.

Differential Effect of Global Signal Regression Between Awake and Anesthetized Conditions in Mice.

Introduction: In resting-state functional magnetic resonance imaging (rs-fMRI) studies, global signal regression (GSR) is a controversial preprocessing strategy. It effectively eliminates global noise driven by motion and respiration but also can introduce artifacts and remove functionally relevant metabolic information. Most preclinical rs-fMRI studies are performed in anesthetized animals, and anesthesia will alter both metabolic and neuronal activity. Methods: In this study, we explored the effect of GSR on rs-fMRI data collected under anesthetized and awake state in mice (n = 12). We measured global signal amplitude, and also functional connectivity (FC), functional connectivity density (FCD) maps, and brain modularity, all commonly used data-driven analysis methods to quantify connectivity patterns. Results: We found that global signal amplitude was similar between the awake and anesthetized states. However, GSR had a different impact on connectivity networks and brain modularity changes between states. We demonstrated that GSR had a more prominent impact on the anesthetized state, with a greater decrease in functional connectivity and increased brain modularity. We classified mice using the change in amplitude of brain modularity coefficient (ΔQ) before and after GSR processing. The results revealed that, when compared with the largest ΔQ group, the smallest ΔQ group had increased FCD in the cortex region in both the awake and anesthetized states. This suggests differences in individual mice may affect how GSR differentially affects awake versus anesthetized functional connectivity. Discussion: This study suggests that, for rs-fMRI studies which compare different physiological states, researchers should use GSR processing with caution.

Sleep fMRI with simultaneous electrophysiology at 9.4 T in male mice.

Sleep is ubiquitous and essential, but its mechanisms remain unclear. Studies in animals and humans have provided insights of sleep at vastly different spatiotemporal scales. However, challenges remain to integrate local and global information of sleep. Therefore, we developed sleep fMRI based on simultaneous electrophysiology at 9.4 T in male mice. Optimized un-anesthetized mouse fMRI setup allowed manifestation of NREM and REM sleep, and a large sleep fMRI dataset was collected and openly accessible. State dependent global patterns were revealed, and state transitions were found to be global, asymmetrical and sequential, which can be predicted up to 17.8 s using LSTM models. Importantly, sleep fMRI with hippocampal recording revealed potentiated sharp-wave ripple triggered global patterns during NREM than awake state, potentially attributable to co-occurrence of spindle events. To conclude, we established mouse sleep fMRI with simultaneous electrophysiology, and demonstrated its capability by revealing global dynamics of state transitions and neural events.

Aberrant brain functional and structural developments in MECP2 duplication rats.

Transgenic animal models with homologous etiology provide a promising way to pursue the neurobiological substrates of the behavioral deficits in autism spectrum disorder (ASD). Gain-of-function mutations of MECP2 cause MECP2 duplication syndrome, a severe neurological disorder with core symptoms of ASD. However, abnormal brain developments underlying the autistic-like behavioral deficits of MECP2 duplication syndrome are rarely investigated. To this end, a human MECP2 duplication (MECP2-DP) rat model was created by the bacterial artificial chromosome transgenic method. Functional and structural magnetic resonance imaging (MRI) with high-field were performed on 16 male MECP2-DP rats and 15 male wildtype rats at postnatal 28 days, 42 days, and 56 days old. Multimodal fusion analyses guided by locomotor-relevant metrics and social novelty time separately were applied to identify abnormal brain networks associated with diverse behavioral deficits induced by MECP2 duplication. Aberrant functional developments of a core network primarily composed of the dorsal medial prefrontal cortex (dmPFC) and retrosplenial cortex (RSP) were detected to associate with diverse behavioral phenotypes in MECP2-DP rats. Altered developments of gray matter volume were detected in the hippocampus and thalamus. We conclude that gain-of-function mutations of MECP2 induce aberrant functional activities in the default-mode-like network and aberrant volumetric changes in the brain, resulting in autistic-like behavioral deficits. Our results gain critical insights into the biomarker of MECP2 duplication syndrome and the neurobiological underpinnings of the behavioral deficits in ASD.

Whole-brain mapping of mouse CSF flow via HEAP-METRIC phase-contrast MRI.

PURPOSE: CSF plays important roles in clearing brain waste and homeostasis. However, mapping whole-brain CSF flow in the rodents is difficult, primarily due to its assumed very low velocity. Therefore, we aimed to develop a novel phase-contrast MRI method to map whole-brain CSF flow in the mouse brain. METHODS: A novel generalized Hadamard encoding-based multi-band scheme (dubbed HEAP-METRIC, Hadamard Encoding APproach of Multi-band Excitation for short TR Imaging aCcelerating) using complex Hadamard matrix was developed and incorporated into conventional phase contrast (PC)-MRI to significantly increase SNR. RESULTS: Slow flow phantom imaging validated HEAP-METRIC PC-MRI's ability to achieve fast and accurate mapping of slow flow velocities (~102  µm/s). With the SNR gain afforded by HEAP-METRIC scheme, high-resolution (0.08 × 0.08 mm in-plane resolution and 36 0.4 mm slices) PC-MRI was completed in 21 min for whole-brain CSF flow mapping in the mouse. Using this novel method, we provide the first report of whole-brain CSF flow in the awake mouse brain with an average flow velocity of ~200 µm/s. Furthermore, HEAP-METRIC PC-MRI revealed CSF flow was reduced by isoflurane anesthesia, accompanied by reduction of glymphatic function as measured by dynamic contrast-enhanced MRI. CONCLUSION: We developed and validated a generalized HEAP-METRIC PC-MRI for mapping low velocity flow. With this method, we have achieved the first whole-brain mapping of awake mouse CSF flow and have further revealed that anesthesia reduces CSF flow velocity.

High-resolution relaxometry-based calibrated fMRI in murine brain: Metabolic differences between awake and anesthetized states.

Functional magnetic resonance imaging (fMRI) techniques using the blood-oxygen level-dependent (BOLD) signal have shown great potential as clinical biomarkers of disease. Thus, using these techniques in preclinical rodent models is an urgent need. Calibrated fMRI is a promising technique that can provide high-resolution mapping of cerebral oxygen metabolism (CMRO2). However, calibrated fMRI is difficult to use in rodent models for several reasons: rodents are anesthetized, stimulation-induced changes are small, and gas challenges induce noisy CMRO2 predictions. We used, in mice, a relaxometry-based calibrated fMRI method which uses cerebral blood flow (CBF) and the BOLD-sensitive magnetic relaxation component, R2', the same parameter derived in the deoxyhemoglobin-dilution model of calibrated fMRI. This method does not use any gas challenges, which we tested on mice in both awake and anesthetized states. As anesthesia induces a whole-brain change, our protocol allowed us to overcome the former limitations of rodent studies using calibrated fMRI. We revealed 1.5-2 times higher CMRO2, dependent upon brain region, in the awake state versus the anesthetized state. Our results agree with alternative measurements of whole-brain CMRO2 in the same mice and previous human anesthesia studies. The use of calibrated fMRI in rodents has much potential for preclinical fMRI.

Automating Perforator Flap MRA and CTA Reporting.

Surgical breast reconstruction after mastectomy requires precise perforator coordinates/dimensions, perforator course, and fat volume in a radiology report. Automatic perforator reporting software was implemented as an OsiriX Digital Imaging and Communications in Medicine (DICOM) viewer plugin. For perforator analysis, the user identifies a reference point (e.g., umbilicus) and marks each perforating artery/vein bundle with multiple region of interest (ROI) points along its course beginning at the muscle-fat interface. Computations using these points and analysis of image data produce content for the report. Post-processing times were compared against conventional/manual methods using de-identified images of 26 patients with surgically confirmed accuracy of perforator locations and caliber. The time from loading source images to completion of report was measured. Significance of differences in mean processing times for this automated approach versus the conventional/manual approach was assessed using a paired t test. The mean conventional reporting time for our radiologists was 76 ± 27 min (median 65 min) compared with 25 ± 6 min (median 25 min) using our OsiriX plugin (p < 0.01). The conventional approach had three reports with transcription errors compared to none with the OsiriX plugin. Otherwise, the reports were similar. In conclusion, automated reporting of perforator magnetic resonance angiography (MRA) studies is faster compared with the standard, manual approach, and transcription errors which are eliminated.

Left ventricle: fully automated segmentation based on spatiotemporal continuity and myocardium information in cine cardiac magnetic resonance imaging (LV-FAST).

CMR quantification of LV chamber volumes typically and manually defines the basal-most LV, which adds processing time and user-dependence. This study developed an LV segmentation method that is fully automated based on the spatiotemporal continuity of the LV (LV-FAST). An iteratively decreasing threshold region growing approach was used first from the midventricle to the apex, until the LV area and shape discontinued, and then from midventricle to the base, until less than 50% of the myocardium circumference was observable. Region growth was constrained by LV spatiotemporal continuity to improve robustness of apical and basal segmentations. The LV-FAST method was compared with manual tracing on cardiac cine MRI data of 45 consecutive patients. Of the 45 patients, LV-FAST and manual selection identified the same apical slices at both ED and ES and the same basal slices at both ED and ES in 38, 38, 38, and 41 cases, respectively, and their measurements agreed within -1.6 ± 8.7 mL, -1.4 ± 7.8 mL, and 1.0 ± 5.8% for EDV, ESV, and EF, respectively. LV-FAST allowed LV volume-time course quantitatively measured within 3 seconds on a standard desktop computer, which is fast and accurate for processing the cine volumetric cardiac MRI data, and enables LV filling course quantification over the cardiac cycle.

Phase-corrected bipolar gradients in multi-echo gradient-echo sequences for quantitative susceptibility mapping.

OBJECTIVE: Large echo spacing of unipolar readout gradients in current multi-echo gradient-echo (GRE) sequences for mapping fields in quantitative susceptibility mapping (QSM) can be reduced using bipolar readout gradients thereby improving acquisition efficiency. MATERIALS AND METHODS: Phase discrepancies between odd and even echoes in the bipolar readout gradients caused by non-ideal gradient behaviors were measured, modeled as polynomials in space and corrected for accordingly in field mapping. The bipolar approach for multi-echo GRE field mapping was compared with the unipolar approach for QSM. RESULTS: The odd-even-echo phase discrepancies were approximately constant along the phase encoding direction and linear along the readout and slice-selection directions. A simple linear phase correction in all three spatial directions was shown to enable accurate QSM of the human brain using a bipolar multi-echo GRE sequence. Bipolar multi-echo acquisition provides QSM in good quantitative agreement with unipolar acquisition while also reducing noise. CONCLUSION: With a linear phase correction between odd-even echoes, bipolar readout gradients can be used in multi-echo GRE sequences for QSM.

Algorithm for fast monoexponential fitting based on Auto-Regression on Linear Operations (ARLO) of data.

PURPOSE: To develop a fast and accurate monoexponential fitting algorithm based on Auto-Regression on Linear Operations (ARLO) of data, and to validate its accuracy and computational speed by comparing it with the conventional Levenberg-Marquardt (LM) and Log-Linear (LL) algorithms. METHODS: ARLO, LM, and LL performances for T2* mapping were evaluated in simulation and in vivo imaging of liver (n=15) and myocardial (n=1) iron overload patients and the brain (two healthy volunteers). RESULTS: In simulations, ARLO consistently delivered accuracy similar to LM and significantly superior to LL. In in vivo mapping of T2 * values, ARLO showed excellent agreement with LM, while LL showed only limited agreements with ARLO and LM. Compared with LM and LL in the liver, ARLO was 125 and 8 times faster using our Matlab implementations, and 156 and 13 times faster using our C++ implementations. In C++ implementations, ARLO reduced the online whole-brain processing time from 9 min 15 s of LM and 35 s of LL to 2.7 s, providing T2 * maps approximately in real time. CONCLUSION: Due to comparable accuracy and significantly higher speed, ARLO can be considered as a valid alternative to the conventional LM algorithm for online T2 * mapping.

Quantitative susceptibility mapping of multiple sclerosis lesions at various ages.

PURPOSE: To assess multiple sclerosis (MS) lesions at various ages by using quantitative susceptibility mapping (QSM) and conventional magnetic resonance (MR) imaging. MATERIALS AND METHODS: Retrospectively selected were 32 clinically confirmed MS patients (nine men and 23 women; 39.3 years ± 10.9) who underwent two MR examinations (interval, 0.43 years ± 0.16) with three-dimensional gradient-echo sequence from August 2011 to August 2012. To estimate the ages of MS lesions, MR examinations performed 0.3-10.6 years before study examinations were studied. Hyperintensity on T2-weighted images was used to define MS lesions. QSM images were reconstructed from gradient-echo data. Susceptibility of MS lesions and temporal rates of change were obtained from QSM images. Lesion susceptibilities were analyzed by t test with intracluster correlation adjustment and Bonferroni correction in multiple comparisons. RESULTS: MR imaging of 32 patients depicted 598 MS lesions, of which 162 lesions (27.1%) in 23 patients were age measurable and six (1.0%) were only visible at QSM. The susceptibilities relative to normal-appearing white matter (NAWM) were 0.53 ppb ± 3.34 for acute enhanced lesions, 38.43 ppb ± 13.0 (positive; P < .01) for early to intermediately aged nonenhanced lesions, and 4.67 ppb ± 3.18 for chronic nonenhanced lesions. Temporal rates of susceptibility changes relative to cerebrospinal fluid were 12.49 ppb/month ± 3.15 for acute enhanced lesions, 1.27 ppb/month ± 2.31 for early to intermediately aged nonenhanced lesions, and -0.004 ppb/month ± 0 for chronic nonenhanced lesions. CONCLUSION: Magnetic susceptibility of MS lesions increased rapidly as it changed from enhanced to nonenhanced, it attained a high susceptibility value relative to NAWM during its initial few years (approximately 4 years), and it gradually dissipated back to susceptibility similar to that of NAWM as it aged, which may provide new insight into pathophysiologic features of MS lesions. Online supplemental material is available for this article.

Reducing the object orientation dependence of susceptibility effects in gradient echo MRI through quantitative susceptibility mapping.

This study demonstrates the dependence of non-local susceptibility effects on object orientation in gradient echo MRI and the reduction of non-local effects by deconvolution using quantitative susceptibility mapping. Imaging experiments were performed on a 3T MRI system using a spoiled 3D multi-echo GRE sequence on phantoms of known susceptibilities, and on human brains of healthy subjects and patients with intracerebral hemorrhages. Magnetic field measurements were determined from multiple echo phase data. To determine the quantitative susceptibility mapping, these field measurements were deconvolved through a dipole inversion kernel under a constraint of consistency with the magnitude images. Phantom and human data demonstrated that the hypointense region in GRE magnitude image corresponding to a susceptibility source increased in volume with TE and varied with the source orientation. The induced magnetic field extended beyond the susceptibility source and varied with its orientation. In quantitative susceptibility mapping, these blooming artifacts, including their dependence on object orientation, were reduced, and the material susceptibilities were quantified.

Graded contribution of hippocampus to multifeature binding across temporal delay.

It is well accepted that all the elements of an episode are combined together by the hippocampus to form an episodic representation. In this study, we further investigated whether the activation pattern of the hippocampus during multifeatural episodic encoding with temporal discontinuities varies with recollective detail in a graded way. Specifically, in the present functional magnetic resonance imaging study, we manipulated two associative features, color and size, and presented the features and the item sequentially. Right hippocampal activation increased with the number of features successfully bound to the item, supporting a graded role of hippocampal activation in bridging temporal discontinuities for integrating multiple episodic features.