Towards automated in vivo parcellation of the human cerebral cortex using supervised classification of magnetic resonance fingerprinting residuals.

BACKGROUND: Accurate parcellation of the cerebral cortex in an individual is a guide to its underlying organization. The most promising in vivo quantitative magnetic resonance (MR)-based microstructural cortical mapping methods are yet to achieve a level of parcellation accuracy comparable to quantitative histology. METHODS: We scanned 6 participants using a 3D echo-planar imaging MR fingerprinting (EPI-MRF) sequence on a 7T Siemens scanner. After projecting MRF signals to the individual-specific inflated model of the cortical surface, normalized autocorrelations of MRF residuals of vertices of 8 microstructurally distinct areas (BA1, BA2, BA4a, BA6, BA44, BA45, BA17, and BA18) from 3 cortical regions were used as feature vector inputs into linear support vector machine (SVM), radial basis function SVM (RBF-SVM), random forest, and k-nearest neighbors supervised classification algorithms. The algorithms' prediction performance was compared using: (i) features from each vertex or (ii) features from neighboring vertices. RESULTS: The neighborhood-based RBF-SVM classifier achieved the highest prediction score of 0.85 for classification of MRF residuals in the central region from a held-out participant. CONCLUSIONS: We developed an automated method of cortical parcellation using a combination of MR fingerprinting residual analysis and machine learning classification. Our findings provide the basis for employing unsupervised learning algorithms for whole-cortex structural parcellation in individuals.

Generalisation of continuous time random walk to anomalous diffusion MRI models with an age-related evaluation of human corpus callosum.

Diffusion MRI measures of the human brain provide key insight into microstructural variations across individuals and into the impact of central nervous system diseases and disorders. One approach to extract information from diffusion signals has been to use biologically relevant analytical models to link millimetre scale diffusion MRI measures with microscale influences. The other approach has been to represent diffusion as an anomalous transport process and infer microstructural information from the different anomalous diffusion equation parameters. In this study, we investigated how parameters of various anomalous diffusion models vary with age in the human brain white matter, particularly focusing on the corpus callosum. We first unified several established anomalous diffusion models (the super-diffusion, sub-diffusion, quasi-diffusion and fractional Bloch-Torrey models) under the continuous time random walk modelling framework. This unification allows a consistent parameter fitting strategy to be applied from which meaningful model parameter comparisons can be made. We then provided a novel way to derive the diffusional kurtosis imaging (DKI) model, which is shown to be a degree two approximation of the sub-diffusion model. This link between the DKI and sub-diffusion models led to a new robust technique for generating maps of kurtosis and diffusivity using the sub-diffusion parameters ßSUB and DSUB. Superior tissue contrast is achieved in kurtosis maps based on the sub-diffusion model. 7T diffusion weighted MRI data for 65 healthy participants in the age range 19-78 years was used in this study. Results revealed that anomalous diffusion model parameters α and ß have shown consistent positive correlation with age in the corpus callosum, indicating α and ß are sensitive to tissue microstructural changes in ageing.

Instant tissue field and magnetic susceptibility mapping from MRI raw phase using Laplacian enhanced deep neural networks.

Quantitative susceptibility mapping (QSM) is an MRI post-processing technique that produces spatially resolved magnetic susceptibility maps from phase data. However, the traditional QSM reconstruction pipeline involves multiple non-trivial steps, including phase unwrapping, background field removal, and dipole inversion. These intermediate steps not only increase the reconstruction time but accumulates errors. This study aims to overcome existing limitations by developing a Laplacian-of-Trigonometric-functions (LoT) enhanced deep neural network for near-instant quantitative field and susceptibility mapping (i.e., iQFM and iQSM) from raw MRI phase data. The proposed iQFM and iQSM methods were compared with established reconstruction pipelines on simulated and in vivo datasets. In addition, experiments on patients with intracranial hemorrhage and multiple sclerosis were also performed to test the generalization of the proposed neural networks. The proposed iQFM and iQSM methods in healthy subjects yielded comparable results to those involving the intermediate steps while dramatically improving reconstruction accuracies on intracranial hemorrhages with large susceptibilities. High susceptibility contrast between multiple sclerosis lesions and healthy tissue was also achieved using the proposed methods. Comparative studies indicated that the most significant contributor to iQFM and iQSM over conventional multi-step methods was the elimination of traditional Laplacian unwrapping. The reconstruction time on the order of minutes for traditional approaches was shortened to around 0.1 s using the trained iQFM and iQSM neural networks.

Improving FLAIR SAR efficiency at 7T by adaptive tailoring of adiabatic pulse power through deep learning B1+ estimation.

PURPOSE: The purpose of this study is to demonstrate a method for specific absorption rate (SAR) reduction for 2D T2 -FLAIR MRI sequences at 7 T by predicting the required adiabatic radiofrequency (RF) pulse power and scaling the RF amplitude in a slice-wise fashion. METHODS: We used a time-resampled frequency-offset corrected inversion (TR-FOCI) adiabatic pulse for spin inversion in a T2 -FLAIR sequence to improve B1+ homogeneity and calculated the pulse power required for adiabaticity slice-by-slice to minimize the SAR. Drawing on the implicit B1+ inhomogeneity in a standard localizer scan, we acquired 3D AutoAlign localizers and SA2RAGE B1+ maps in 28 volunteers. Then, we trained a convolutional neural network (CNN) to estimate the B1+ profile from the localizers and calculated pulse scale factors for each slice. We assessed the predicted B1+ profiles and the effect of scaled pulse amplitudes on the FLAIR inversion efficiency in oblique transverse, sagittal, and coronal orientations. RESULTS: The predicted B1+ amplitude maps matched the measured ones with a mean difference of 9.5% across all slices and participants. The slice-by-slice scaling of the TR-FOCI inversion pulse was most effective in oblique transverse orientation and resulted in a 1 min and 30 s reduction in SAR induced delay time while delivering identical image quality. CONCLUSION: We propose a SAR reduction technique based on the estimation of B1+ profiles from standard localizer scans using a CNN and show that scaling the inversion pulse power slice-by-slice for FLAIR sequences at 7T reduces SAR and scan time without compromising image quality.

Can anomalous diffusion models in magnetic resonance imaging be used to characterise white matter tissue microstructure?

PURPOSE: During the time window of diffusion weighted magnetic resonance imaging experiments (DW-MRI), water diffusion in tissue appears to be anomalous as a transient effect, with a mean squared displacement that is not a linear function of time. A number of statistical models have been proposed to describe water diffusion in tissue, and parameters describing anomalous as well as Gaussian diffusion have previously been related to measures of tissue microstructure such as mean axon radius. We analysed the relationship between white matter tissue characteristics and parameters of existing statistical diffusion models. METHODS: A white matter tissue model (ActiveAx) was used to generate multiple b-value diffusion-weighted magnetic resonance imaging signals. The following models were evaluated to fit the diffusion signal: 1) Gaussian models - 1a) mono-exponential decay and 1b) bi-exponential decay; 2) Anomalous diffusion models - 2a) stretched exponential, 2b) continuous time random walk and 2c) space fractional Bloch-Torrey equation. We identified the best candidate model based on the relationship between the diffusion-derived parameters and mean axon radius and axial diffusivity, and applied it to the in vivo DW-MRI data acquired at 7.0 T from five healthy participants to estimate the same selected tissue characteristics. Differences between simulation parameters and fitted parameters were used to assess accuracy and in vivo findings were compared to previously reported observations. RESULTS: The space fractional Bloch-Torrey model was found to be the best candidate in characterising white matter on the base of the ActiveAx simulated DW-MRI data. Moreover, parameters of the space fractional Bloch-Torrey model were sensitive to mean axon radius and axial diffusivity and exhibited low noise sensitivity based on simulations. We also found spatial variations in the model parameter ß to reflect changes in mean axon radius across the mid-sagittal plane of the corpus callosum. CONCLUSION: Simulations have been used to define how the parameters of the most common statistical magnetic resonance imaging diffusion models relate to axon radius and diffusivity. The space fractional Bloch-Torrey equation was identified as the best model for the characterisation of axon radius and diffusivity. This model allows changes in mean axon radius and diffusivity to be inferred from spatially resolved maps of model parameters.

Assessment of microstructural signal compartments across the corpus callosum using multi-echo gradient recalled echo at 7 T.

Quantitative assessment of tissue microstructure is important in studying human brain diseases and disorders in which white matter is implicated, as it has been linked to demyelination, re-myelination, and axonal damage in clinical conditions. Ultra-high field magnetic resonance imaging data obtained using a multi-echo gradient echo sequence has been shown to contain information on myelin, axonal and extracellular compartments in white matter. In this study, we aimed to assess the sensitivity of a three-compartment model to estimate the variation of corresponding compartment parameters (water fraction, relaxation time and frequency shift) of the corpus callosum sub-regions, which are known to have different tissue structure. Additionally, we computed the g-ratio using myelin and axonal water fractions and performed a voxel-by-voxel analysis in the corpus callosum. Based on data acquired for ten participants, we show that the myelin compartment water fraction and T2∗ is consistent across the corpus callosum sub-regions, whilst myelin frequency shift varies. The results show that the variation in water fraction, T2∗ and frequency shift for the myelin signal compartment across the corpus callosum is smaller than for the axonal and extracellular signal compartments. The computed g-ratio was comparable to previously published studies in the corpus callosum. Our study suggests that a multi-echo GRE approach in vivo combined with a complex three-compartment model is sensitive to microstructural parameter variations across the human corpus callosum.

Signal compartments in ultra-high field multi-echo gradient echo MRI reflect underlying tissue microstructure in the brain.

Gradient recalled echo magnetic resonance imaging (GRE-MRI) at ultra-high field holds great promise for new contrast mechanisms and delineation of putative tissue compartments that contribute to the multi-echo GRE-MRI signal may aid structural characterization. Several studies have adopted the three water-pool compartment model to study white matter brain regions, associating individual compartments with myelin, axonal and extracellular water. However, the number and identifiability of GRE-MRI signal compartments has not been fully explored. We undertook this task for human brain imaging data. Multiple echo time GRE-MRI data were acquired in five healthy participants, specific anatomical structures were segmented in each dataset (substantia nigra, caudate, insula, putamen, thalamus, fornix, internal capsule, corpus callosum and cerebrospinal fluid), and the signal fitted with models comprising one to six signal compartments using a complex-valued plane wave formulation. Information criteria and cluster analysis methods were used to ascertain the number of distinct compartments within the signal from each structure and to determine their respective frequency shifts. We identified five principal signal compartments with different relative contributions to each structure's signal. Voxel-based maps of the volume fraction of each of these compartments were generated and demonstrated spatial correlation with brain anatomy.

The challenge of bias-free coil combination for quantitative susceptibility mapping at ultra-high field.

PURPOSE: Quantitative susceptibility mapping is a technique to estimate the magnetic property of tissue with particularly high sensitivity at ultra-high field. However, a key challenge at ultra-high field is the combination of phase data acquired using phased array receive coils. Several methods for combining phase data have been proposed, but the influence of coil combination choices on susceptibility quantitation has not been studied systematically. METHODS: We combined phase data using COMPOSER (COMbining Phase data using a Short Echo-time Reference scan) and a reference-free channel-by-channel method. We investigated the effect of the chosen combination method on susceptibility results in a group of 28 participants at 7 Tesla. RESULTS: Our results show that reference scans can bias susceptibility values. Although the proposed reference-free channel-by-channel method cannot remove transmit field phase, it shows comparable results to the COMPOSER method in which a high-resolution ultrashort echo-time reference scan was used. CONCLUSIONS: We conclude that ultrashort echo-time reference scans reduce quantitation bias and remove the transmit field phase when using COMPOSER to combine phase data, and not combining the phase data before susceptibility processing avoids this bias, resulting in comparable results. Magn Reson Med 79:97-107, 2018. © 2017 InternationalSociety for Magnetic Resonance in Medicine.

Utility of gradient recalled echo magnetic resonance imaging for the study of myelination in cuprizone mice treated with fingolimod.

The availability of high-field-strength magnetic resonance imaging (MRI) systems has brought about the development of techniques that aim to map myelination via the exploitation of various contrast mechanisms. Myelin mapping techniques have the potential to provide tools for the diagnosis and treatment of diseases, such as multiple sclerosis. In this study, we evaluated the sensitivity of T2 *, frequency shift and susceptibility measures to myelin levels in a cuprizone mouse model of demyelination. The model was supplemented with two different dosages of fingolimod, a drug known to positively affect demyelination. A decrease in grey-white matter contrast with the cuprizone diet was observed for T2 *, frequency shift and susceptibility measures, together with myelin basic protein antibody findings. These results indicate that T2 *, frequency shift and susceptibility measures have the potential to act as biomarkers for myelination. Susceptibility was found to be the most sensitive measure to changes in grey-white matter contrast. In addition, fingolimod treatment was found to reduce the level of demyelination, with a larger dosage exhibiting a greater reduction in demyelination for the in vivo MRI results. Overall, susceptibility mapping appears to be a more promising tool than T2 * or frequency shift mapping for the early diagnosis and treatment of diseases in which myelination is implicated.

Tissue microstructure features derived from anomalous diffusion measurements in magnetic resonance imaging.

OBJECTIVES: Tissue microstructure features, namely axon radius and volume fraction, provide important information on the function of white matter pathways. These parameters vary on the scale much smaller than imaging voxels (microscale) yet influence the magnetic resonance imaging diffusion signal at the image voxel scale (macroscale) in an anomalous manner. Researchers have already mapped anomalous diffusion parameters from magnetic resonance imaging data, but macroscopic variations have not been related to microscale influences. With the aid of a tissue model, we aimed to connect anomalous diffusion parameters to axon radius and volume fraction using diffusion-weighted magnetic resonance imaging measurements. EXPERIMENTAL DESIGN: An ex vivo human brain experiment was performed to directly validate axon radius and volume fraction measurements in the human brain. These findings were validated using electron microscopy. Additionally, we performed an in vivo study on nine healthy participants to map axon radius and volume fraction along different regions of the corpus callosum projecting into various cortical areas identified using tractography. PRINCIPAL OBSERVATIONS: We found a clear relationship between anomalous diffusion parameters and axon radius and volume fraction. We were also able to map accurately the trend in axon radius along the corpus callosum, and in vivo findings resembled the low-high-low-high behaviour in axon radius demonstrated previously. CONCLUSIONS: Axon radius and volume fraction measurements can potentially be used in brain connectivity studies and to understand the implications of white matter structure in brain diseases and disorders. Hum Brain Mapp 38:1068-1081, 2017. © 2016 Wiley Periodicals, Inc.

Characterization of anomalous relaxation using the time-fractional Bloch equation and multiple echo T2 *-weighted magnetic resonance imaging at 7 T.

PURPOSE: To study the utility of fractional calculus in modeling gradient-recalled echo MRI signal decay in the normal human brain. METHODS: We solved analytically the extended time-fractional Bloch equations resulting in five model parameters, namely, the amplitude, relaxation rate, order of the time-fractional derivative, frequency shift, and constant offset. Voxel-level temporal fitting of the MRI signal was performed using the classical monoexponential model, a previously developed anomalous relaxation model, and using our extended time-fractional relaxation model. Nine brain regions segmented from multiple echo gradient-recalled echo 7 Tesla MRI data acquired from five participants were then used to investigate the characteristics of the extended time-fractional model parameters. RESULTS: We found that the extended time-fractional model is able to fit the experimental data with smaller mean squared error than the classical monoexponential relaxation model and the anomalous relaxation model, which do not account for frequency shift. CONCLUSIONS: We were able to fit multiple echo time MRI data with high accuracy using the developed model. Parameters of the model likely capture information on microstructural and susceptibility-induced changes in the human brain. Magn Reson Med 77:1485-1494, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Echo time-dependent quantitative susceptibility mapping contains information on tissue properties.

PURPOSE: Magnetic susceptibility is a physical property of matter that varies depending on chemical composition and abundance of different molecular species. Interest is growing in mapping of magnetic susceptibility in the human brain using magnetic resonance imaging techniques, but the influences affecting the mapped values are not fully understood. METHODS: We performed quantitative susceptibility mapping on 7 Tesla (T) multiple echo time gradient recalled echo data and evaluated the trend in 10 regions of the human brain. Temporal plots of susceptibility were performed in the caudate, pallidum, putamen, thalamus, insula, red nucleus, substantia nigra, internal capsule, corpus callosum, and fornix. We implemented an existing three compartment signal model and used optimization to fit the experimental result to assess the influences that could be responsible for our findings. RESULTS: The temporal trend in susceptibility is different for different brain regions, and subsegmentation of specific regions suggests that differences are likely to be attributable to variations in tissue structure and composition. Using a signal model, we verified that a nonlinear temporal behavior in experimentally computed susceptibility within imaging voxels may be the result of the heterogeneous composition of tissue properties. CONCLUSIONS: Decomposition of voxel constituents into meaningful parameters may lead to informative measures that reflect changes in tissue microstructure. Magn Reson Med 77:1946-1958, 2017. © 2016 International Society for Magnetic Resonance in Medicine.

Selective channel combination of MRI signal phase.

PURPOSE: Signal magnitude can robustly be combined using the sum-of-squares approach. Methods have been developed to combine complex images. However, techniques based only on signal phase have not been developed and evaluated. THEORY AND METHODS: We performed simulations to demonstrate the effect of noise on coil combination. 32-channel 7 Tesla human gradient echo MRI brain data were collected. We combined phase images based on phase noise leading to spatially selective and coil selective combination of phase images. We compared our selective combination approach to optimal noise distribution and adaptive combination methods. RESULTS: We found that selective combination of signal phases leads to improved phase signal-to-noise ratio. Furthermore, a phase shift can be present in combined phase images introduced by the method used to combine multiple channel phases. CONCLUSION: Mapping of signal phase from ultra-high field MRI data undoubtedly provides a wealth of information about the ageing brain and the effects of neurodegenerative disorders. Measurement of signal phase is essential in frequency shift mapping and in quantitative susceptibility mapping. The method used to combine signal phase should be informed by an understanding of the noise distribution in signal phase at the individual channel level. Magn Reson Med 76:1469-1477, 2016. © 2015 International Society for Magnetic Resonance in Medicine.

MRI signal phase oscillates with neuronal activity in cerebral cortex: implications for neuronal current imaging.

Neuronal activity produces transient ionic currents that may be detectable using magnetic resonance imaging (MRI). We examined the feasibility of MRI-based detection of neuronal currents using computer simulations based on the laminar cortex model (LCM). Instead of simulating the activity of single neurons, we decomposed neuronal activity to action potentials (AP) and postsynaptic potentials (PSP). The geometries of dendrites and axons were generated dynamically to account for diverse neuronal morphologies. Magnetic fields associated with APs and PSPs were calculated during spontaneous and stimulated cortical activity, from which the neuronal current induced MRI signal was determined. We found that the MRI signal magnitude change (<0.1 ppm) is below currently detectable levels but that the signal phase change is likely to be detectable. Furthermore, neuronal MRI signals are sensitive to temporal and spatial variations in neuronal activity but independent of the intensity of neuronal activation. Synchronised neuronal activity produces large phase changes (in the order of 0.1 mrad). However, signal phase oscillates with neuronal activity. Consequently, MRI scans need to be synchronised with neuronal oscillations to maximise the likelihood of detecting signal phase changes due to neuronal currents. These findings inform the design of MRI experiments to detect neuronal currents.

Image synthesis of interictal SPECT from MRI and PET using machine learning.

Background: Cross-modality image estimation can be performed using generative adversarial networks (GANs). To date, SPECT image estimation from another medical imaging modality using this technique has not been considered. We evaluate the estimation of SPECT from MRI and PET, and additionally assess the necessity for cross-modality image registration for GAN training. Methods: We estimated interictal SPECT from PET and MRI as a single-channel input, and as a multi-channel input to the GAN. We collected data from 48 individuals with epilepsy and converted them to 3D isotropic images for consistence across the modalities. Training and testing data were prepared in native and template spaces. The Pix2pix framework within the GAN network was adopted. We evaluated the addition of the structural similarity index metric to the loss function in the GAN implementation. Root-mean-square error, structural similarity index, and peak signal-to-noise ratio were used to assess how well SPECT images were able to be synthesised. Results: High quality SPECT images could be synthesised in each case. On average, the use of native space images resulted in a 5.4% percentage improvement in SSIM than the use of images registered to template space. The addition of structural similarity index metric to the GAN loss function did not result in improved synthetic SPECT images. Using PET in either the single channel or dual channel implementation led to the best results, however MRI could produce SPECT images close in quality. Conclusion: Synthesis of SPECT from MRI or PET can potentially reduce the number of scans needed for epilepsy patient evaluation and reduce patient exposure to radiation.

A short 18F-FDG imaging window triple injection neuroimaging protocol for parametric mapping in PET.

BACKGROUND: In parametric PET, kinetic parameters are extracted from dynamic PET images. It is not commonly used in clinical practice because of long scan times and the requirement for an arterial input function (AIF). To address these limitations, we designed an 18F-fluorodeoxyglucose (18F-FDG) triple injection dynamic PET protocol for brain imaging with a standard field of view PET scanner using a 24-min imaging window and an input function modeled using measurements from a region of interest placed over the left ventricle. METHODS: To test the protocol in 6 healthy participants, we examined the quality of voxel-based maps of kinetic parameters in the brain generated using the two-tissue compartment model and compared estimated parameter values with previously published values. We also utilized data from a 36-min validation imaging window to compare (1) the modeled AIF against the input function measured in the validation window; and (2) the net influx rate ([Formula: see text]) computed using parameter estimates from the short imaging window against the net influx rate obtained using Patlak analysis in the validation window. RESULTS: Compared to the AIF measured in the validation window, the input function estimated from the short imaging window achieved a mean area under the curve error of 9%. The voxel-wise Pearson's correlation between [Formula: see text] estimates from the short imaging window and the validation imaging window exceeded 0.95. CONCLUSION: The proposed 24-min triple injection protocol enables parametric 18F-FDG neuroimaging with noninvasive estimation of the AIF from cardiac images using a standard field of view PET scanner.

Automated extraction of the arterial input function from brain images for parametric PET studies.

BACKGROUND: Accurate measurement of the arterial input function (AIF) is crucial for parametric PET studies, but the AIF is commonly derived from invasive arterial blood sampling. It is possible to use an image-derived input function (IDIF) obtained by imaging a large blood pool, but IDIF measurement in PET brain studies performed on standard field of view scanners is challenging due to lack of a large blood pool in the field-of-view. Here we describe a novel automated approach to estimate the AIF from brain images. RESULTS: Total body 18F-FDG PET data from 12 subjects were split into a model adjustment group (n = 6) and a validation group (n = 6). We developed an AIF estimation framework using wavelet-based methods and unsupervised machine learning to distinguish arterial and venous activity curves, compared to the IDIF from the descending aorta. All of the automatically extracted AIFs in the validation group had similar shape to the IDIF derived from the descending aorta IDIF. The average area under the curve error and normalised root mean square error across validation data were - 1.59 ± 2.93% and 0.17 ± 0.07. CONCLUSIONS: Our automated AIF framework accurately estimates the AIF from brain images. It reduces operator-dependence, and could facilitate the clinical adoption of parametric PET.

ParaPET: non-invasive deep learning method for direct parametric brain PET reconstruction using histoimages.

BACKGROUND: The indirect method for generating parametric images in positron emission tomography (PET) involves the acquisition and reconstruction of dynamic images and temporal modelling of tissue activity given a measured arterial input function. This approach is not robust, as noise in each dynamic image leads to a degradation in parameter estimation. Direct methods incorporate into the image reconstruction step both the kinetic and noise models, leading to improved parametric images. These methods require extensive computational time and large computing resources. Machine learning methods have demonstrated significant potential in overcoming these challenges. But they are limited by the requirement of a paired training dataset. A further challenge within the existing framework is the use of state-of-the-art arterial input function estimation via temporal arterial blood sampling, which is an invasive procedure, or an additional magnetic resonance imaging (MRI) scan for selecting a region where arterial blood signal can be measured from the PET image. We propose a novel machine learning approach for reconstructing high-quality parametric brain images from histoimages produced from time-of-flight PET data without requiring invasive arterial sampling, an MRI scan, or paired training data from standard field-of-view scanners. RESULT: The proposed is tested on a simulated phantom and five oncological subjects undergoing an 18F-FDG-PET scan of the brain using Siemens Biograph Vision Quadra. Kinetic parameters set in the brain phantom correlated strongly with the estimated parameters (K1, k2 and k3, Pearson correlation coefficient of 0.91, 0.92 and 0.93) and a mean squared error of less than 0.0004. In addition, our method significantly outperforms (p < 0.05, paired t-test) the conventional nonlinear least squares method in terms of contrast-to-noise ratio. At last, the proposed method was found to be 37% faster than the conventional method. CONCLUSION: We proposed a direct non-invasive DL-based reconstruction method and produced high-quality parametric maps of the brain. The use of histoimages holds promising potential for enhancing the estimation of parametric images, an area that has not been extensively explored thus far. The proposed method can be applied to subject-specific dynamic PET data alone.

The laminar cortex model: a new continuum cortex model incorporating laminar architecture.

Local field potentials (LFPs) are widely used to study the function of local networks in the brain. They are also closely correlated with the blood-oxygen-level-dependent signal, the predominant contrast mechanism in functional magnetic resonance imaging. We developed a new laminar cortex model (LCM) to simulate the amplitude and frequency of LFPs. Our model combines the laminar architecture of the cerebral cortex and multiple continuum models to simulate the collective activity of cortical neurons. The five cortical layers (layer I, II/III, IV, V, and VI) are simulated as separate continuum models between which there are synaptic connections. The LCM was used to simulate the dynamics of the visual cortex under different conditions of visual stimulation. LFPs are reported for two kinds of visual stimulation: general visual stimulation and intermittent light stimulation. The power spectra of LFPs were calculated and compared with existing empirical data. The LCM was able to produce spontaneous LFPs exhibiting frequency-inverse (1/ƒ) power spectrum behaviour. Laminar profiles of current source density showed similarities to experimental data. General stimulation enhanced the oscillation of LFPs corresponding to gamma frequencies. During simulated intermittent light stimulation, the LCM captured the fundamental as well as high order harmonics as previously reported. The power spectrum expected with a reduction in layer IV neurons, often observed with focal cortical dysplasias associated with epilepsy was also simulated.

Hippocampal demyelination is associated with increased magnetic susceptibility in a mouse model of concussion.

Structural and functional deficits in the hippocampus are a prominent feature of moderate-severe traumatic brain injury (TBI). In this work, we investigated the potential of Quantitative Susceptibility Imaging (QSM) to reveal the temporal changes in myelin integrity in a mouse model of concussion (mild TBI). We employed a cross-sectional design wherein we assigned 43 mice to cohorts undergoing either a concussive impact or a sham procedure, with QSM imaging at day 2, 7, or 14 post-injury, followed by Luxol Fast Blue (LFB) myelin staining to assess the structural integrity of hippocampal white matter (WM). We assessed spatial learning in the mice using the Active Place Avoidance Test (APA), recording their ability to use visual cues to locate and avoid zone-dependent mild electrical shocks. QSM and LFB staining indicated changes in the stratum lacunosum-molecular layer of the hippocampus in the concussion groups, suggesting impairment of this key relay between the entorhinal cortex and the CA1 regions. These imaging and histology findings were consistent with demyelination, namely increased magnetic susceptibility to MR imaging and decreased LFB staining. In the APA test, sham animals showed fewer entries into the shock zone compared to the concussed cohort. Thus, we present radiological, histological, and behavioral findings that concussion can induce significant and alterations in hippocampal integrity and function that evolve over time after the injury.