BRAHMA: Population specific T1, T2, and FLAIR weighted brain templates and their impact in structural and functional imaging studies

Differences in brain morphology across population groups necessitate creation of population-specific Magnetic Resonance Imaging (MRI) brain templates for interpretation of neuroimaging data. Variations in the neuroanatomy in a genetically heterogeneous population make the development of a population-specific brain template for the Indian subcontinent imperative. A dataset of high-resolution 3D T1, T2-weighted, and FLAIR images acquired from a group of 113 volunteers (M/F - 56/57, mean age-28.96 ± 7.80 years) are used to construct T1, T2-weighted, and FLAIR templates, collectively referred to as Indian Brain Template, “BRAHMA”. A processing pipeline is developed and implemented in a MATLAB based toolbox for template construction and generation of tissue probability maps and segmentation atlases, with additional labels for deep brain regions such as the Substantia Nigra generated from the T2-weighted and FLAIR templates. The use of BRAHMA template for analysis of structural and functional neuroimaging data obtained from Indian participants, provides improved accuracy with statistically significant results over that obtained using the ICBM-152 (International Consortium for Brain Mapping) template. Our results indicate that segmentations generated on structural images are closer in volume to those obtained from registration to the BRAHMA template than to the ICBM-152. Furthermore, functional MRI data obtained for Working Memory and Finger Tapping paradigms processed using the BRAHMA template show a significantly higher percentage of the activation area than ICBM-152 in relevant brain regions, i.e. the left middle frontal gyrus, and the left and right precentral gyri, respectively. The availability of different image contrasts, tissue maps, and segmentation atlases makes the BRAHMA template a comprehensive tool for multi-modal image analysis in laboratory and clinical settings.     

A deep-learning semantic segmentation approach to fully automated MRI-based left-ventricular deformation analysis in cardiotoxicity

Left-ventricular (LV) strain measurements with the Displacement Encoding with Stimulated Echoes (DENSE) MRI sequence provide accurate estimates of cardiotoxicity damage related to breast cancer chemotherapy. This study investigated an automated LV chamber quantification tool via segmentation with a supervised deep convolutional neural network (DCNN) before strain analysis with DENSE images. Segmentation for chamber quantification analysis was conducted with a custom DeepLabV3+ DCNN with ResNet-50 backbone on 42 female breast cancer datasets (22 training-sets, eight validation-sets and 12 independent test-sets). Parameters such as LV end-diastolic diameter (LVEDD) and ejection fraction (LVEF) were quantified, and myocardial strains analyzed with the Radial Point Interpolation Method (RPIM). Myocardial classification was validated against ground-truth with sensitivity-specificity analysis, the metrics of Dice, average perpendicular distance (APD) and Hausdorff-distance. Following segmentation, validation was conducted with the Cronbach's Alpha (C-Alpha) intraclass correlation coefficient between LV chamber quantification results with DENSE and Steady State Free Precession (SSFP) acquisitions and a vendor tool-based method to segment the DENSE data, and similarly for myocardial strain analysis in the chambers. The results of myocardial classification from segmentation of the DENSE data were accuracy = 97%, Dice = 0.89 and APD = 2.4 mm in the test-set. The C-Alpha correlations from comparing chamber quantification results between the segmented DENSE and SSFP data and vendor tool-based method were 0.97 for LVEF (56 ± 7% vs 55 ± 7% vs 55 ± 6%, p = 0.6) and 0.77 for LVEDD (4.6 ± 0.4 cm vs 4.5 ± 0.3 cm vs 4.5 ± 0.3 cm, p = 0.8). The validation metrics against ground-truth and equivalent parameters obtained from the SSFP segmentation and vendor tool-based comparisons show that the DCNN approach is applicable for automated LV chamber quantification and subsequent strain analysis in cardiotoxicity.     

An energy minimization method for MS lesion segmentation from T1-w and FLAIR images

In this paper, we extend the multiplicative intrinsic component optimization (MICO) algorithm to multichannel MR image segmentation, with focus on segmentation of multiple sclerosis (MS) lesions. The MICO algorithm was originally proposed by Li et al. in Ref. [1] for normal brain tissue segmentation and intensity inhomogeneity correction of a single channel MR image, which exhibits desirable advantages over other methods for MR image segmentation and intensity inhomogeneity correction in terms of segmentation accuracy and robustness. In this paper, we extend the MICO algorithm to multi-channel MR image segmentation and enable the segmentation of MS lesions. We assign different weights for different channels to control the impact of each channel. The weighted channels allow the enhancement of the impact of the FLAIR image on the segmentation of MS lesions by assigning a larger weight to the FLAIR image channel than the other channels. With the inherent mechanism of estimation of the bias field, our method is able to deal with the intensity inhomogeneity in the input multi-channel MR images. In the application of our method, we only use T1-w and FLAIR images as the input two channel MR images. Experimental results show promising result of our method.     

Meningiomas: Preoperative predictive histopathological grading based on radiomics of MRI

PurposeWe aimed to develop a radiomics model to predict the histopathological grading of meningiomas by magnetic resonance imaging (MRI) before surgery.MethodsWe recruited 131 patients with pathological diagnosis of meningiomas. All the patients had undergone MRI before surgery on a 3.0 T MRI scanner to obtain T1 fluid- attenuated inversion recovery (T1 FLAIR) images, T2-weighted images (T2WI) and T1 FLAIR with contrast enhancement (CE-T1 FLAIR) images covering the whole brain. The removing features with low variance, univariate feature selection, and least absolute shrinkage and selection operator (LASSO) were used to select radiomics features. Six classifiers were used to train the models (logistic regression (LR), k-nearest neighbor (KNN), decision tree (DT), support vector machine (SVM), random forests (RF), and XGBoost), and then 24 models were established using a random verification method to differentiate low-grade from high-grade meningiomas. The performance was assessed by receiver-operating characteristic (ROC) analysis, the f1-score, sensitivity, and specificity.ResultsThe radiomics features were significantly associated with the histopathological grading. Quantitative imaging features (n = 1409) were extracted, and nine features were selected to predict the grades of meningiomas. The best performance of the radiomics model for the degree of differentiation was obtained by SVM (area under the curve (AUC), 0.956; 95% confidence interval (CI), 0.83–1.00; sensitivity, 0.87; specificity, 0.92; f1-score, 0.90).ConclusionThe radiomics models are of great value in predicting the histopathological grades of meningiomas, and have broad prospects in radiology and clinics.     

Deep unregistered multi-contrast MRI reconstruction

Multiple magnetic resonance images of different contrasts are normally acquired for clinical diagnosis. Recently, research has shown that the previously acquired multi-contrast (MC) images of the same patient can be used as anatomical prior to accelerating magnetic resonance imaging (MRI). However, current MC-MRI networks are based on the assumption that the images are perfectly registered, which is rarely the case in real-world applications. In this paper, we propose an end-to-end deep neural network to reconstruct highly accelerated images by exploiting the shareable information from potentially misaligned reference images of an arbitrary contrast. Specifically, a spatial transformation (ST) module is designed and integrated into the reconstruction network to align the pre-acquired reference images with the images to be reconstructed. The misalignment is further alleviated by maximizing the normalized cross-correlation (NCC) between the MC images. The visualization of feature maps demonstrates that the proposed method effectively reduces the misalignment between the images for shareable information extraction when applied to the publicly available brain datasets. Additionally, the experimental results on these datasets show the proposed network allows the robust exploitation of shareable information across the misaligned MC images, leading to improved reconstruction results.     

Brain imaging: Comparison of T1W FLAIR BLADE with conventional T1W SE

IntroductionAlthough T1 weighted spin echo (T1W SE) images are widely used to study anatomical details and pathologic abnormalities of the brain, its role in delineation of lesions and reduction of artifacts has not been thoroughly investigated. BLADE is a fairly new technique that has been reported to reduce motion artifacts and improve image quality.ObjectiveThe primary objective of this study is to compare the quality of T1-weighted fluid attenuated inversion recovery (FLAIR) images with BLADE technique (T1W FLAIR BLADE) and the quality of T1W SE images in the MR imaging of the brain. The goal is to highlight the advantages of the two sequences as well as which one can better reduce flow and motion artifacts so that the imaging of the lesions will not be impaired.Materials and methodsBrain examinations with T1W FLAIR BLADE and T1W SE sequences were performed on 48 patients using a 1.5 T scanner. These techniques were evaluated by two radiologists based on: a) a qualitative analysis i.e. overall image quality, presence of artifacts, CSF nulling; and b) a quantitative analysis of signal-to-noise ratios (SNR), contrast-to-noise ratios (CNR) and Relative Contrast. The statistical analysis was performed using the Kruskal-Wallis non-parametric system.ResultsIn the qualitative analysis, BLADE sequences had a higher scoring than the conventional sequences in all the cases. The overall image quality was better on T1W FLAIR BLADE. Motion and flow-related artifacts were lower in T1W FLAIR BLADE. Regarding the SNR measurements, T1W SE appeared to have higher values in the majority of cases, whilst T1W-FLAIR BLADE had higher values in the CNR and Relative Contrast measurements.ConclusionT1W FLAIR BLADE sequence appears to be superior to T1W SE in overall image quality and reduction of motion and flow-pulsation artifacts as well as in nulling CSF and has been preferred by the clinicians. T1W FLAIR BLADE may be an alternative approach in brain MRI imaging.     

Automatic knee cartilage segmentation from multi-contrast MR images using support vector machine classification with spatial dependencies

Accurate segmentation of knee cartilage is required to obtain quantitative cartilage measurements, which is crucial for the assessment of knee pathology caused by musculoskeletal diseases or sudden injuries. This paper presents an automatic knee cartilage segmentation technique which exploits a rich set of image features from multi-contrast magnetic resonance (MR) images and the spatial dependencies between neighbouring voxels. The image features and the spatial dependencies are modelled into a support vector machine (SVM)-based association potential and a discriminative random field (DRF)-based interaction potential. Subsequently, both potentials are incorporated into an inference graphical model such that the knee cartilage segmentation is cast into an optimal labelling problem which can be efficiently solved by loopy belief propagation. The effectiveness of the proposed technique is validated on a database of multi-contrast MR images. The experimental results show that using diverse forms of image and anatomical structure information as the features are helpful in improving the segmentation, and the joint SVM-DRF model is superior to the classification models based solely on DRF or SVM in terms of accuracy when the same features are used. The developed segmentation technique achieves good performance compared with gold standard segmentations and obtained higher average DSC values than the state-of-the-art automatic cartilage segmentation studies.     

Are multi-contrast magnetic resonance images necessary for segmenting multiple sclerosis brains? A large cohort study based on deep learning

BackgroundMagnetic resonance images with multiple contrasts or sequences are commonly used for segmenting brain tissues, including lesions, in multiple sclerosis (MS). However, acquisition of images with multiple contrasts increases the scan time and complexity of the analysis, possibly introducing factors that could compromise segmentation quality.ObjectiveTo investigate the effect of various combinations of multi-contrast images as input on the segmented volumes of gray (GM) and white matter (WM), cerebrospinal fluid (CSF), and lesions using a deep neural network.MethodsU-net, a fully convolutional neural network was used to automatically segment GM, WM, CSF, and lesions in 1000 MS patients. The input to the network consisted of 15 combinations of FLAIR, T1-, T2-, and proton density-weighted images. The Dice similarity coefficient (DSC) was evaluated to assess the segmentation performance. For lesions, true positive rate (TPR) and false positive rate (FPR) were also evaluated. In addition, the effect of lesion size on lesion segmentation was investigated.ResultsHighest DSC was observed for all the tissue volumes, including lesions, when the input was combination of all four image contrasts. All other input combinations that included FLAIR also provided high DSC for all tissue classes. However, the quality of lesion segmentation showed strong dependence on the input images. The DSC and TPR values for inputs with the four contrast combination and FLAIR alone were very similar, but FLAIR showed a moderately higher FPR for lesion size <100 μl. For lesions smaller than 20 μl all image combinations resulted in poor performance. The segmentation quality improved with lesion size.ConclusionsBest performance for segmented tissue volumes was obtained with all four image contrasts as the input, and comparable performance was attainable with FLAIR only as the input, albeit with a moderate increase in FPR for small lesions. This implies that acquisition of only FLAIR images provides satisfactory tissue segmentation. Lesion segmentation was poor for very small lesions and improved rapidly with lesion size.     

Prognosticating brain tumor patient survival after laser thermotherapy: Comparison between neuroradiological reading and semi-quantitative analysis of MRI data

Background and purposeGiven increasing interest in laser interstitial thermotherapy (LITT) to treat brain tumor patients, we explored if examining multiple MRI contrasts per brain tumor patient undergoing surgery can impact predictive accuracy of survival post-LITT.Materials and methodsMRI contrasts included fluid-attenuated inversion recovery (FLAIR), T1 pre-gadolinium (T1pre), T1 post-gadolinium (T1Gd), T2, diffusion-weighted imaging (DWI), apparent diffusion coefficient (ADC), susceptibility weighted images (SWI), and magnetization-prepared rapid gradient-echo (MPRAGE). The latter was used for MRI data registration across preoperative to postoperative scans. Two ROIs were identified by thresholding preoperative FLAIR (large ROI) and T1Gd (small ROI) images. For each MRI contrast, a numerical score was assigned based on changing image intensity of both ROIs (vs. a normal ROI) from preoperative to postoperative stages. The fully-quantitative method was based on changing image intensity across scans at different stages without any human intervention, whereas the semi-quantitative method was based on subjective criteria of cumulative trends across scans at different stages. A fully-quantitative/semi-quantitative score per patient was obtained by averaging scores for each MRI contrast. A standard neuroradiological reading score per patient was obtained from radiological interpretation of MRI data. Scores from all 3 methods per patient were compared against patient survival, and re-examined for comorbidity and pathology effects.ResultsPatient survival correlated best with semi-quantitative scores obtained from T1Gd, ADC, and T2 data, and these correlations improved when biopsy and comorbidity were included.ConclusionThese results suggest interfacing neuroradiological readings with semi-quantitative image analysis can improve predictive accuracy of patient survival.     

Signal characteristics of FLAIR related to water content: comparison with conventional spin echo imaging in infarcted rat brain

To examine the correlation between tissue water content and signal intensity on fluid-attenuated inversion recovery (FLAIR) images, we analyzed infarcted rat brain, verified the results by theoretical simulation, and compared them with conventional spin-echo images. We produced brain infarction with cavitation in five rats by middle cerebral artery occlusion. After in vivo MRI, histologic sections of the MRI plane were obtained. We measured the signal intensity of regions on FLAIR and spin-echo images, and measured the area of cavitation on histologic sections. We plotted curves of cavity percentage to signal intensity. Theoretical values were calculated using a two-compartment model. On the curve of cavity area to signal intensity, the signal on FLAIR images peaked in tissues with 20% to 30% area of cavitation. On the theoretical curve, the signal on FLAIR images peaked at 90% tissue water content. These results seem to be characteristic of FLAIR.     

PDE-based spatial smoothing: a practical demonstration of impacts on MRI brain extraction, tissue segmentation and registration

Spatial smoothing is typically used to denoise magnetic resonance imaging (MRI) data. Gaussian smoothing kernels, associated with heat equations or isotropic diffusion (ISD), are widely adopted for this purpose because of their easy implementation and efficient computation, but despite these advantages, Gaussian smoothing kernels blur the edges, curvature and texture of images. To overcome these issues, researchers have proposed anisotropic diffusion (ASD) and non-local means [i.e., diffusion (NLD)] kernels. However, these new filtering paradigms are rarely applied to MRI analyses. In the current study, using real degraded MRI data, we demonstrated the effect of denoising using ISD, ASD and NLD kernels. Furthermore, we evaluated their impact on three common preprocessing steps of MRI data analysis: brain extraction, segmentation and registration. Results suggest that NLD-based spatial smoothing is most effective at improving the quality of MRI data preprocessing and thus should become the new standard method of smoothing in MRI data processing.     

First-principles and Monte Carlo combinational study on Zn1−xCoxO diluted magnetic semiconductor

The electronic structures and magnetic properties of Zn_1−xCo_xO (x=5.55%,8.33%,12.5%) are studied using first-principles calculations in combination with Monte Carlo (MC) simulation. The combinational method makes possible a complete simulation from the microscopic magnetic interaction to macroscopic magnetic behavior. The calculated results from first principles indicate that the ferromagnetic ground state is stabilized by a half-metallic electronic structure which originates from the strong hybridization between Co 3d electrons and O 2p electrons. With the magnetic coupling strengths obtained from first-principles calculations, the MC simulation predicts the ferromagnetism of Zn_1−xCo_xO (x=5.55%,8.33%,12.5%) with , which is consistent with the experimental facts.     

Automatic white matter lesion segmentation using an adaptive outlier detection method

White matter (WM) lesions are diffuse WM abnormalities that appear as hyperintense (bright) regions in cranial magnetic resonance imaging (MRI). WM lesions are often observed in older populations and are important indicators of stroke, multiple sclerosis, dementia and other brain-related disorders. In this paper, a new automated method for WM lesions segmentation is presented. In the proposed method, the presence of WM lesions is detected as outliers in the intensity distribution of the fluid-attenuated inversion recovery (FLAIR) MR images using an adaptive outlier detection approach. Outliers are detected using a novel adaptive trimmed mean algorithm and box-whisker plot. In addition, pre- and postprocessing steps are implemented to reduce false positives attributed to MRI artifacts commonly observed in FLAIR sequences. The approach is validated using the cranial MRI sequences of 38 subjects. A significant correlation (R=0.9641, P value=3.12×10(-3)) is observed between the automated approach and manual segmentation by radiologist. The accuracy of the proposed approach was further validated by comparing the lesion volumes computed using the automated approach and lesions manually segmented by an expert radiologist. Finally, the proposed approach is compared against leading lesion segmentation algorithms using a benchmark dataset.     

Tsallis-Entropy Segmentation through MRF and Alzheimer anatomic reference for Brain Magnetic Resonance Parcellation

Quantifying the intracranial tissue volume changes in magnetic resonance imaging (MRI) assists specialists to analyze the effects of natural or pathological changes. Since these changes can be subtle, the accuracy of the automatic compartmentalization method is always criticized by specialists. We propose and then evaluate an automatic segmentation method based on modified q-entropy (Mqe) through a modified Markov Random Field (MMRF) enhanced by Alzheimer anatomic reference (AAR) to provide a high accuracy brain tissues parcellation approach (Mqe-MMRF). We underwent two strategies to evaluate Mqe-MMRF; a simulation of different levels of noise and non-uniformity effect on MRI data (7 subjects) and a set of twenty MRI data available from MRBrainS13 as patient brain tissue segmentation challenge. We accessed eleven quality metrics compared to reference tissues delineations to evaluate Mqe-MMRF. MRI segmentation scores decreased by only 4.6% on quality metrics after noise and non-uniformity simulations of 40% and 9%, respectively. We found significant mean improvements in the metrics of the five training subjects, for whole-brain 0.86%, White Matter 3.20%, Gray Matter 3.99%, and Cerebrospinal Fluid 4.16% (p-values < 0.02) when Mqe-MMRF compared to the other reference methods. We also processed the Mqe-MMRF on 15 evaluation subjects group from MRBrainS13 online challenge, and the results held a higher rank than the reference tools; FreeSurfer, SPM, and FSL. Since the proposed method improved the precision of brain segmentation, specifically, for GM, and thus one can use it in quantitative and morphological brain studies.     

Use of random matrix theory in the discovery of resting state brain networks

Connectomics identifies brain networks in vivo in resting state functional MRI. However, the presence of noise produces spurious identification of brain networks, which have low test-retest reliability. A Network Based Statistics approach to network identification has been previously proposed that affords much better statistical power relative to Bonferroni method but nevertheless provides a sufficiently conservative, family-wise control for false positives. We propose the use of Random Matrix Theory (RMT) to discover brain networks and to associate those networks with demographic and clinical variables. We parcellated the brain into cortical and subcortical regions using either an anatomical or a functional brain atlas. We applied RMT to study functional connectivity across brain regions by first computing the correlation matrix for time courses in those brain regions and then identifying eigenvalues that deviate from the theoretical random distribution that RMT predicts, on the assumption that real brain networks would produce eigenvalues that differ significantly from the random distribution. We assessed the specificity and test-retest reliability of identified networks through application of this RMT-based approach to (1) synthetic data generated under the null-hypothesis, (2) resting state functional MRI data from 4 real-world cohorts of patients and healthy controls, and (3) synthetic data generated by the addition of increasing amounts of noise to real-world datasets. Our findings showed that RMT method was robust to the atlas used for parcellating the brain and did not discover a brain network in synthetic data when in fact a network was not present (i.e., specificity was high); RMT-identified networks in the real-world dataset had high test-retest reliability; and RMT-based method consistently discovered the same network in the presence of increasing noise in the real-world dataset.     

A workflow for automated segmentation of the liver surface,hepatic vasculature and biliary tree anatomy from multiphase MR images

Accurate assessment of 3D models of patient-specific anatomy of the liver, including underlying hepatic and biliary tree, is critical for preparation and safe execution of complex liver resections, especially due to high variability of biliary and hepatic artery anatomies. Dynamic MRI with hepatospecific contrast agents is currently the only type of diagnostic imaging that provides all anatomical information required for generation of such a model, yet there is no information in the literature on how the complete 3D model can be generated automatically. In this work, a new automated segmentation workflow for extraction of patient-specific 3D model of the liver, hepatovascular and biliary anatomy from a single multiphase MRI acquisition is developed and quantitatively evaluated. The workflow incorporates course 4D k-means clustering estimation and geodesic active contour refinement of the liver boundary, based on organ’s characteristic uptake of gadolinium contrast agents overtime. Subsequently, hepatic vasculature and biliary ducts segmentations are performed using multiscale vesselness filters. The algorithm was evaluated using 15 test datasets of patients with liver malignancies of various histopathological types. It showed good correlation with expert manual segmentation, resulting in an average of 1.76 ± 2.44 mm Hausdorff distance for the liver boundary, and 0.58 ± 0.72 and 1.16 ± 1.98 mm between centrelines of biliary ducts and liver veins, respectively. A workflow for automatic segmentation of the liver, hepatic vasculature and biliary anatomy from a single diagnostic MRI acquisition was developed. This enables automated extraction of 3D models of patient-specific liver anatomy, and may facilitating better perception of organ’s anatomy during preparation and execution of liver surgeries. Additionally, it may help to reduce the incidence of intraoperative biliary duct damage due to an unanticipated variation in the anatomy.     

Optimizing Image Reconstruction in SENSE Using GPU

Parallel magnetic resonance imaging (MRI) (pMRI) uses multiple receiver coils to reduce the MRI scan time. To accelerate the data acquisition process in MRI, less amount of data is acquired from the scanner which leads to artifacts in the reconstructed images. SENSitivity Encoding (SENSE) is a reconstruction algorithm in pMRI to remove aliasing artifacts from the undersampled multi coil data and recovers fully sampled images. The main limitation of SENSE is computing inverse of the encoding matrix. This work proposes the inversion of encoding matrix using Jacobi singular value decomposition (SVD) algorithm for image reconstruction on GPUs to accelerate the reconstruction process. The performance of Jacobi SVD is compared with Gauss–Jordan algorithm. The simulations are performed on two datasets (brain and cardiac) with acceleration factors 2, 4, 6 and 8. The results show that the graphics processing unit (GPU) provides a speed up to 21.6 times as compared to CPU reconstruction. Jacobi SVD algorithm performs better in terms of acceleration in reconstructions on GPUs as compared to Gauss–Jordan method. The proposed algorithm is suitable for any number of coils and acceleration factors for SENSE reconstruction on real time processing systems.     

In-vivo diffusion MRI protocol optimization for the chimpanzee brain and examination of aging effects on the primate optic nerve at 3T

BackgroundDiffusion MRI (dMRI) data acquisition protocols are well-established on modern high-field clinical scanners for human studies. However, these protocols are not suitable for the chimpanzee (or other large-brained mammals) because of its substantial difference in head geometry and brain volume compared with humans. Therefore, an optimal dMRI data acquisition protocol dedicated to chimpanzee neuroimaging is needed.MethodsA multi-shot (4 segments) double spin-echo echo-planar imaging (MS-EPI) sequence and a single-shot double spin-echo EPI (SS-EPI) sequence were optimized separately for in vivo dMRI data acquisition of chimpanzees using a clinical 3T scanner. Correction for severe susceptibility-induced image distortion and signal drop-off of the chimpanzee brain was performed and evaluated using FSL software. DTI indices in different brain regions and probabilistic tractography were compared. A separate DTI data set from n=34 chimpanzees (13 to 56 years old) was collected using the optimal protocol. Age-related changes in diffusivity indices of optic nerve fibers were evaluated.ResultsThe SS-EPI sequence acquired dMRI data of the chimpanzee brain with approximately doubled the SNR as the MS-EPI sequence given the same scan time. The quality of white matter fiber tracking from the SS-EPI data was much higher than that from MS-EPI data. However, quantitative analysis of DTI indices showed no difference in most ROIs between the SS-EPI and MS-EPI sequences. The progressive evolution of diffusivity indices of optic nerves indicated mild changes in fiber bundles of chimpanzees aged 40 years and above.ConclusionThe single-shot EPI-based acquisition protocol provided better image quality of dMRI for chimpanzee brains and is recommended for in vivo dMRI study or clinical diagnosis of chimpanzees (or other large animals) using a clinical scanner. Also, the tendency of FA decrease or diffusivity increase in the optic nerve of aged chimpanzees was seen but did not show significant age-related changes, suggesting aging may have less impact on optic nerve fiber integrity of chimpanzees, in contrast to previous results for both macaque monkeys and humans.     

Stable Atlas-based Mapped Prior (STAMP) machine-learning segmentation for multicenter large-scale MRI data

Machine learning (ML)-based segmentation methods are a common technique in the medical image processing field. In spite of numerous research groups that have investigated ML-based segmentation frameworks, there remains unanswered aspects of performance variability for the choice of two key components: ML algorithm and intensity normalization. This investigation reveals that the choice of those elements plays a major part in determining segmentation accuracy and generalizability. The approach we have used in this study aims to evaluate relative benefits of the two elements within a subcortical MRI segmentation framework. Experiments were conducted to contrast eight machine-learning algorithm configurations and 11 normalization strategies for our brain MR segmentation framework. For the intensity normalization, a Stable Atlas-based Mapped Prior (STAMP) was utilized to take better account of contrast along boundaries of structures. Comparing eight machine learning algorithms on down-sampled segmentation MR data, it was obvious that a significant improvement was obtained using ensemble-based ML algorithms (i.e., random forest) or ANN algorithms. Further investigation between these two algorithms also revealed that the random forest results provided exceptionally good agreement with manual delineations by experts. Additional experiments showed that the effect of STAMP-based intensity normalization also improved the robustness of segmentation for multicenter data sets. The constructed framework obtained good multicenter reliability and was successfully applied on a large multicenter MR data set (n > 3000). Less than 10% of automated segmentations were recommended for minimal expert intervention. These results demonstrate the feasibility of using the ML-based segmentation tools for processing large amount of multicenter MR images. We demonstrated dramatically different result profiles in segmentation accuracy according to the choice of ML algorithm and intensity normalization chosen.     

A unified framework for mapping individual interregional high-order morphological connectivity based on regional cortical features from anatomical MRI

Building individual brain networks form the single volume of anatomical MRI is a challenging task. Furthermore, the high-order connectivity of morphological networks remains unexplored. This paper aimed to investigate the individual high-order morphological connectivity from anatomical MRI. Towards this goal, a unified framework based on six feature distances (euclidean, seuclidean, mahalanobis, cityblock, minkowski, and chebychev) was proposed to derive high-order interregional morphological features. The test-retest datasets and the healthy aging datasets were applied to analyze the reliability and the inter-subject variability of the novel features. In addition, the predictive models based on these novel features were established for age estimation. The proposed six neuroanatomical features exhibited significant high-to-excellent reliability. Certain connections were significantly correlated to biological age based on the six novel metrics (p < .05, FDR corrected). Moreover, the predicted age were significantly correlated to the original age in each regression task (r > 0.5, p < 10⁻⁶). The results suggested that the novel high-order metrics were reliable and could reflect individual differences, which could be beneficial for current methods of individual brain connectomes.