Objective Evaluation of Multiple Sclerosis Lesion Segmentation using a Data Management and Processing Infrastructure

We present a study of multiple sclerosis segmentation algorithms conducted at the international MICCAI 2016 challenge. This challenge was operated using a new open-science computing infrastructure. This allowed for the automatic and independent evaluation of a large range of algorithms in a fair and completely automatic manner. This computing infrastructure was used to evaluate thirteen methods of MS lesions segmentation, exploring a broad range of state-of-theart algorithms, against a high-quality database of 53 MS cases coming from four centers following a common definition of the acquisition protocol. Each case was annotated manually by an unprecedented number of seven different experts. Results of the challenge highlighted that automatic algorithms, including the recent machine learning methods (random forests, deep learning, …), are still trailing human expertise on both detection and delineation criteria. In addition, we demonstrate that computing a statistically robust consensus of the algorithms performs closer to human expertise on one score (segmentation) although still trailing on detection scores.

Semi-Automated Extraction of Crohns Disease MR Imaging Markers using a 3D Residual CNN with Distance Prior

We propose a 3D residual convolutional neural network (CNN) algorithm with an integrated distance prior for segmenting the small bowel lumen and wall to enable extraction of pediatric Crohns disease (pCD) imaging markers from T1-weighted contrast-enhanced MR images. Our proposed segmentation framework enables, for the first time, to quantitatively assess luminal narrowing and dilation in CD aimed at optimizing surgical decisions as well as analyzing bowel wall thickness and tissue enhancement for assessment of response to therapy. Given seed points along the bowel lumen, the proposed algorithm automatically extracts 3D image patches centered on these points and a distance map from the interpolated centerline. These 3D patches and corresponding distance map are jointly used by the proposed residual CNN architecture to segment the lumen and the wall, and to extract imaging markers. Due to lack of available training data, we also propose a novel and efficient semi-automated segmentation algorithm based on graph-cuts technique as well as a software tool for quickly editing labeled data that was used to train our proposed CNN model. The method which is based on curved planar reformation of the small bowel is also useful for visualizing, manually refining, and measuring pCD imaging markers. In preliminary experiments, our CNN network obtained Dice coefficients of 75 ± 18%, 81 ± 8% and 97 ± 2% for the lumen, wall and background, respectively.

Localization of stereo-electroencephalography signals using a finite difference complete electrode model

Surgical intervention in epilepsy aims to eliminate seizures in refractory patients by resecting the tissue responsible for seizure onset. Stereo-electroencephalography (sEEG) provides highly accurate but invasive electrophysiological measurements using narrow multi-contact electrodes implanted stereotactically through small holes in the skull. However, the three dimensional nature of sEEG measurements make observed seizure onsets difficult to associate with physical cortical regions. Three dimensional source localization from sEEG measurements can improve the interpretation of this data, but requires more accurate modeling as compared to localization from scalp EEG. Here, we present a finite difference approach that models the contact impedance and physical extent of each electrode (the so-called complete electrode model), to localize brain electrical activity from sEEG measurements. We applied this model to MRI and CT in a patient with intractable epilepsy, and reconstructed activity associated with multiple types of recurrent ictal spikes observed in sEEG. Independently, the neurosurgeon resected the clinically determined seizure focus, creating a resection cavity, and rendering the patient free of seizures. Our localization placed the seizure focus at a focal region in the occipital lobe, entirely contained within the resection region.

Head motion measurement and correction using FID navigators

PURPOSE

To develop a novel framework for rapid, intrinsic head motion measurement in MRI using FID navigators (FIDnavs) from a multichannel head coil array.

METHODS

FIDnavs encode substantial rigid-body motion information; however, current implementations require patient-specific training with external tracking data to extract quantitative positional changes. In this work, a forward model of FIDnav signals was calibrated using simulated movement of a reference image within a model of the spatial coil sensitivities. A FIDnav module was inserted into a nonselective 3D FLASH sequence, and rigid-body motion parameters were retrospectively estimated every readout time using nonlinear optimization to solve the inverse problem posed by the measured FIDnavs. This approach was tested in simulated data and in 7 volunteers, scanned at 3T with a 32-channel head coil array, performing a series of directed motion paradigms.

RESULTS

FIDnav motion estimates achieved mean absolute errors of 0.34 ± 0.49 mm and 0.52 ± 0.61° across all subjects and scans, relative to ground-truth motion measurements provided by an electromagnetic tracking system. Retrospective correction with FIDnav motion estimates resulted in substantial improvements in quantitative image quality metrics across all scans with intentional head motion.

CONCLUSIONS

Quantitative rigid-body motion information can be effectively estimated using the proposed FIDnav-based approach, which represents a practical method for retrospective motion compensation in less cooperative patient populations.

Presurgical language fMRI: Clinical practices and patient outcomes in epilepsy surgical planning

The goal of this study was to document current clinical practice and report patient outcomes in presurgical language functional MRI (fMRI) for epilepsy surgery. Epilepsy surgical programs worldwide were surveyed as to the utility, implementation, and efficacy of language fMRI in the clinic; 82 programs responded. Respondents were predominantly US (61%) academic programs (85%), and evaluated adults (44%), adults and children (40%), or children only (16%). Nearly all (96%) reported using language fMRI. Surprisingly, fMRI is used to guide surgical margins (44% of programs) as well as lateralize language (100%). Sites using fMRI for localization most often use a distance margin around activation of 10mm. While considered useful, 56% of programs reported at least one instance of disagreement with other measures. Direct brain stimulation typically confirmed fMRI findings (74%) when guiding margins, but instances of unpredicted decline were reported by 17% of programs and 54% reported unexpected preservation of function. Programs reporting unexpected decline did not clearly differ from those which did not. Clinicians using fMRI to guide surgical margins do not typically map known language-critical areas beyond Broca's and Wernicke's. This initial data shows many clinical teams are confident using fMRI not only for language lateralization but also to guide surgical margins. Reported cases of unexpected language preservation when fMRI activation is resected, and cases of language decline when it is not, emphasize a critical need for further validation. Comprehensive studies comparing commonly-used fMRI paradigms to predict stimulation mapping and post-surgical language decline remain of high importance.

AUTOMATIC RENAL SEGMENTATION IN DCE-MRI USING CONVOLUTIONAL NEURAL NETWORKS

Kidney function evaluation using dynamic contrast-enhanced MRI (DCE-MRI) images could help in diagnosis and treatment of kidney diseases of children. Automatic segmentation of renal parenchyma is an important step in this process. In this paper, we propose a time and memory efficient fully automated segmentation method which achieves high segmentation accuracy with running time in the order of seconds in both normal kidneys and kidneys with hydronephrosis. The proposed method is based on a cascaded application of two 3D convolutional neural networks that employs spatial and temporal information at the same time in order to learn the tasks of localization and segmentation of kidneys, respectively. Segmentation performance is evaluated on both normal and abnormal kidneys with varying levels of hydronephrosis. We achieved a mean dice coefficient of 91.4 and 83.6 for normal and abnormal kidneys of pediatric patients, respectively.

Cerebellar volume as an imaging marker of development in infants with tuberous sclerosis complex

OBJECTIVE

In this cohort analysis, we studied 1-year-old infants with tuberous sclerosis complex (TSC), correlating volumes of cerebellar structures with neurodevelopmental measures.

METHODS

We analyzed data from a prospective biomarker study in infants with TSC (ClinicalTrials.gov NCT01780441). We included participants aged 12 months with an identified mutation of TSC1 or TSC2. Using MRI segmentation performed with the PSTAPLE algorithm, we measured relative volumes (structure volume divided by intracranial contents volume) of the following structures: right/left cerebellar white matter, right/left cerebellar exterior, vermal lobules I-V, vermal lobules VI-VII, and vermal lobules VIII-X. We correlated relative volumes to Mullen Scales of Early Learning (MSEL) scores.

RESULTS

There were 70 participants (mean age 1.03 [0.11] years): n = 11 had a TSC1 mutation; n = 59 had a TSC2 mutation. For patients with TSC2 mutation, for every percentage increase in total cerebellar volume, there was an approximate 10-point increase in MSEL composite score (β = 10.47 [95% confidence interval 5.67, 15.27], p < 0.001). For patients with TSC1 mutation, the relationship between cerebellar volume and MSEL composite score was not statistically significant (β = -10.88 [95% confidence interval -22.16, 0.41], p = 0.06). For patients with TSC2 mutation, there were positive slopes when regressing expressive language and visual reception skills with volumes of nearly all cerebellar structures (p ≤ 0.29); there were also positive slopes when regressing receptive language skills, gross motor skills, and fine motor skills with volumes of cerebellar right/left exterior (p ≤ 0.014).

CONCLUSIONS

Cerebellar volume loss-perhaps reflecting Purkinje cell degeneration-may predict neurodevelopmental severity in patients with TSC2 mutations.

Fetal Neuropathology in Zika Virus-Infected Pregnant Female Rhesus Monkeys

The development of interventions to prevent congenital Zika syndrome (CZS) has been limited by the lack of an established nonhuman primate model. Here we show that infection of female rhesus monkeys early in pregnancy with Zika virus (ZIKV) recapitulates many features of CZS in humans. We infected 9 pregnant monkeys with ZIKV, 6 early in pregnancy (weeks 6-7 of gestation) and 3 later in pregnancy (weeks 12-14 of gestation), and compared findings with uninfected controls. 100% (6 of 6) of monkeys infected early in pregnancy exhibited prolonged maternal viremia and fetal neuropathology, including fetal loss, smaller brain size, and histopathologic brain lesions, including microcalcifications, hemorrhage, necrosis, vasculitis, gliosis, and apoptosis of neuroprogenitor cells. High-resolution MRI demonstrated concordant lesions indicative of deep gray matter injury. We also observed spinal, ocular, and neuromuscular pathology. Our data show that vascular compromise and neuroprogenitor cell dysfunction are hallmarks of CZS pathogenesis, suggesting novel strategies to prevent and to treat this disease.

Differential Rates of Perinatal Maturation of Human Primary and Nonprimary Auditory Cortex

Primary and nonprimary cerebral cortex mature along different timescales; however, the differences between the rates of maturation of primary and nonprimary cortex are unclear. Cortical maturation can be measured through changes in tissue microstructure detectable by diffusion magnetic resonance imaging (MRI). In this study, diffusion tensor imaging (DTI) was used to characterize the maturation of Heschl's gyrus (HG), which contains both primary auditory cortex (pAC) and nonprimary auditory cortex (nAC), in 90 preterm infants between 26 and 42 weeks postmenstrual age (PMA). The preterm infants were in different acoustical environments during their hospitalization: 46 in open ward beds and 44 in single rooms. A control group consisted of 15 term-born infants. Diffusion parameters revealed that (1) changes in cortical microstructure that accompany cortical maturation had largely already occurred in pAC by 28 weeks PMA, and (2) rapid changes were taking place in nAC between 26 and 42 weeks PMA. At term equivalent PMA, diffusion parameters for auditory cortex were different between preterm infants and term control infants, reflecting either delayed maturation or injury. No effect of room type was observed. For the preterm group, disturbed maturation of nonprimary (but not primary) auditory cortex was associated with poorer language performance at age two years.

Investigating the maturation of microstructure and radial orientation in the preterm human cortex with diffusion MRI

Preterm birth disrupts and alters the complex developmental processes in the cerebral cortex. This disruption may be a contributing factor to widespread delay and cognitive difficulties in the preterm population. Diffusion-weighted magnetic resonance imaging (DW MRI) is a noninvasive imaging technique that makes inferences about cellular structures, at scales smaller than the imaging resolution. One established finding is that DW MRI shows a transient radial alignment in the preterm cortex. In this study, we quantify this maturational process with the "radiality index", a parameter that measures directional coherence, which we expect to change rapidly in the perinatal period. To measure this index, we used structural T²-weighted MRI to segment the cortex and generate cortical meshes. We obtained normal vectors for each face of the mesh and compared them to the principal diffusion direction, calculated by both the DTI and DIAMOND models, to generate the radiality index. The subjects included in this study were 89 infants born at fewer than 34 weeks completed gestation, each imaged at up to four timepoints between 27 and 42 weeks gestational age. In this manuscript, we quantify the longitudinal trajectory of radiality, fractional anisotropy and mean diffusivity from the DTI and DIAMOND models. For the radiality index and fractional anisotropy, the DIAMOND model offers improved sensitivity over the DTI model. The radiality index has a consistent progression across time, with the rate of change depending on the cortical lobe. The occipital lobe changes most rapidly, and the frontal and temporal least: this is commensurate with known developmental anatomy. Analysing the radiality index offers information complementary to other diffusion parameters.

Reproducibility of Brain MRI Segmentation Algorithms: Empirical Comparison of Local MAP PSTAPLE, FreeSurfer, and FSL-FIRST

BACKGROUND AND PURPOSE

Segmentation of human brain structures is crucial for the volumetric quantification of brain disease. Advances in algorithmic approaches have led to automated techniques that save time compared to interactive methods. Recently, the utility and accuracy of template library fusion algorithms, such as Local MAP PSTAPLE (PSTAPLE), have been demonstrated but there is little guidance regarding its reproducibility compared to single template-based algorithms such as FreeSurfer and FSL-FIRST.

METHODS

Eight repeated magnetic resonance imagings of 20 subjects were segmented using FreeSurfer, FSL-FIRST, and PSTAPLE. We reported the reproducibility of segmentation-derived volume measurements for brain structures and calculated sample size estimates for detecting hypothetical rates of tissue atrophy given the observed variances.

RESULTS

PSTAPLE had the most reproducible volume measurements for hippocampus, putamen, thalamus, caudate, pallidum, amygdala, Accumbens area, and cortical regions. FreeSurfer was most reproducible for brainstem. PSTAPLE was the most accurate algorithm in terms of several metrics include Dice's coefficient. The sample size estimates showed that a study utilizing PSTAPLE would require tens to hundreds less subjects than the other algorithms for detecting atrophy rates typically observed in brain disease.

CONCLUSIONS

PSTAPLE is a useful tool for automatic human brain segmentation due to its precision and accuracy, which enable the detection of the size of the effect typically reported for neurological disorders with a substantially reduced sample size, in comparison to the other tools we assessed. This enables randomized controlled trials to be executed with reduced cost and duration, in turn, facilitating the assessment of new therapeutic interventions.

Motion-Robust Spatially Constrained Parameter Estimation in Renal Diffusion-Weighted MRI by 3D Motion Tracking and Correction of Sequential Slices

In this work, we introduce a novel motion-robust spatially constrained parameter estimation (MOSCOPE) technique for kidney diffusion-weighted MRI. The proposed motion compensation technique does not require a navigator, trigger, or breath-hold but only uses the intrinsic features of the acquired data to track and compensate for motion to reconstruct precise models of the renal diffusion signal. We have developed a technique for physiological motion tracking based on robust state estimation and sequential registration of diffusion sensitized slices acquired within 200ms. This allows a sampling rate of 5Hz for state estimation in motion tracking that is sufficiently faster than both respiratory and cardiac motion rates in children and adults, which range between 0.8 to 0.2Hz, and 2.5 to 1Hz, respectively. We then apply the estimated motion parameters to data from each slice and use motion-compensated data for 1) robust intra-voxel incoherent motion (IVIM) model estimation in the kidney using a spatially constrained model fitting approach, and 2) robust weighted least squares estimation of the diffusion tensor model. Experimental results, including precision of IVIM model parameters using bootstrap-sampling and in-vivo whole kidney tractography, showed significant improvement in precision and accuracy of these models using the proposed method compared to models based on the original data and volumetric registration.

Diffusion MRI microstructure models with in vivo human brain Connectome data: results from a multi-group comparison

A large number of mathematical models have been proposed to describe the measured signal in diffusion-weighted (DW) magnetic resonance imaging (MRI). However, model comparison to date focuses only on specific subclasses, e.g. compartment models or signal models, and little or no information is available in the literature on how performance varies among the different types of models. To address this deficiency, we organized the 'White Matter Modeling Challenge' during the International Symposium on Biomedical Imaging (ISBI) 2015 conference. This competition aimed to compare a range of different kinds of models in their ability to explain a large range of measurable in vivo DW human brain data. Specifically, we assessed the ability of models to predict the DW signal accurately for new diffusion gradients and b values. We did not evaluate the accuracy of estimated model parameters, as a ground truth is hard to obtain. We used the Connectome scanner at the Massachusetts General Hospital, using gradient strengths of up to 300 mT/m and a broad set of diffusion times. We focused on assessing the DW signal prediction in two regions: the genu in the corpus callosum, where the fibres are relatively straight and parallel, and the fornix, where the configuration of fibres is more complex. The challenge participants had access to three-quarters of the dataset and their models were ranked on their ability to predict the remaining unseen quarter of the data. The challenge provided a unique opportunity for a quantitative comparison of diverse methods from multiple groups worldwide. The comparison of the challenge entries reveals interesting trends that could potentially influence the next generation of diffusion-based quantitative MRI techniques. The first is that signal models do not necessarily outperform tissue models; in fact, of those tested, tissue models rank highest on average. The second is that assuming a non-Gaussian (rather than purely Gaussian) noise model provides little improvement in prediction of unseen data, although it is possible that this may still have a beneficial effect on estimated parameter values. The third is that preprocessing the training data, here by omitting signal outliers, and using signal-predicting strategies, such as bootstrapping or cross-validation, could benefit the model fitting. The analysis in this study provides a benchmark for other models and the data remain available to build up a more complete comparison in the future.

Assessing the validity of the approximation of diffusion-weighted-MRI signals from crossing fascicles by sums of signals from single fascicles

PURPOSE

To assess the validity of the superposition approximation for crossing fascicles, i.e., the assumption that the total diffusion-weighted MRI signal is the sum of the signals arising from each fascicle independently, even when the fascicles intermingle in a voxel.

METHODS

Monte Carlo simulations were used to study the impact of the approximation on the diffusion-weighted MRI signal and to assess whether this approximate model allows microstructural features of interwoven fascicles to be accurately estimated, despite signal differences.

RESULTS

Small normalized signal differences were observed, typically 10-3-10-2. The use of the approximation had little impact on the estimation of the crossing angle, the axonal density index, and the radius index in clinically realistic scenarios wherein the acquisition noise was the predominant source of errors. In the absence of noise, large systematic errors due to the superposition approximation only persisted for the radius index, mainly driven by a low sensitivity of diffusion-weighted MRI signals to small radii in general.

CONCLUSION

The use of the superposition approximation rather than a model of interwoven fascicles does not adversely impact the estimation of microstructural features of interwoven fascicles in most current clinical settings. Magn Reson Med 79:2332-2345, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Block-Matching Distortion Correction of Echo-Planar Images With Opposite Phase Encoding Directions

By shortening the acquisition time of MRI, Echo Planar Imaging (EPI) enables the acquisition of a large number of images in a short time, compatible with clinical constraints as required for diffusion or functional MRI. However such images are subject to large, local distortions disrupting their correspondence with the underlying anatomy. The correction of those distortions is an open problem, especially in regions where large deformations occur. We propose a new block-matching registration method to perform EPI distortion correction based on the acquisition of two EPI with opposite phase encoding directions (PED). It relies on new transformations between blocks adapted to the EPI distortion model, and on an adapted optimization scheme to ensure an opposite symmetric transformation. We present qualitative and quantitative results of the block-matching correction using different metrics on a phantom dataset and on in-vivo data. We show the ability of the block-matching to robustly correct EPI distortion even in strongly affected areas.

A New Sparse Representation Framework for Reconstruction of an Isotropic High Spatial Resolution MR Volume From Orthogonal Anisotropic Resolution Scans

In magnetic resonance (MR), hardware limitations, scan time constraints, and patient movement often result in the acquisition of anisotropic 3-D MR images with limited spatial resolution in the out-of-plane views. Our goal is to construct an isotropic high-resolution (HR) 3-D MR image through upsampling and fusion of orthogonal anisotropic input scans. We propose a multiframe super-resolution (SR) reconstruction technique based on sparse representation of MR images. Our proposed algorithm exploits the correspondence between the HR slices and the low-resolution (LR) sections of the orthogonal input scans as well as the self-similarity of each input scan to train pairs of overcomplete dictionaries that are used in a sparse-land local model to upsample the input scans. The upsampled images are then combined using wavelet fusion and error backprojection to reconstruct an image. Features are learned from the data and no extra training set is needed. Qualitative and quantitative analyses were conducted to evaluate the proposed algorithm using simulated and clinical MR scans. Experimental results show that the proposed algorithm achieves promising results in terms of peak signal-to-noise ratio, structural similarity image index, intensity profiles, and visualization of small structures obscured in the LR imaging process due to partial volume effects. Our novel SR algorithm outperforms the nonlocal means (NLM) method using self-similarity, NLM method using self-similarity and image prior, self-training dictionary learning-based SR method, averaging of upsampled scans, and the wavelet fusion method. Our SR algorithm can reduce through-plane partial volume artifact by combining multiple orthogonal MR scans, and thus can potentially improve medical image analysis, research, and clinical diagnosis.

Temporal slice registration and robust diffusion-tensor reconstruction for improved fetal brain structural connectivity analysis

Diffusion weighted magnetic resonance imaging, or DWI, is one of the most promising tools for the analysis of neural microstructure and the structural connectome of the human brain. The application of DWI to map early development of the human connectome in-utero, however, is challenged by intermittent fetal and maternal motion that disrupts the spatial correspondence of data acquired in the relatively long DWI acquisitions. Fetuses move continuously during DWI scans. Reliable and accurate analysis of the fetal brain structural connectome requires careful compensation of motion effects and robust reconstruction to avoid introducing bias based on the degree of fetal motion. In this paper we introduce a novel robust algorithm to reconstruct in-vivo diffusion-tensor MRI (DTI) of the moving fetal brain and show its effect on structural connectivity analysis. The proposed algorithm involves multiple steps of image registration incorporating a dynamic registration-based motion tracking algorithm to restore the spatial correspondence of DWI data at the slice level and reconstruct DTI of the fetal brain in the standard (atlas) coordinate space. A weighted linear least squares approach is adapted to remove the effect of intra-slice motion and reconstruct DTI from motion-corrected data. The proposed algorithm was tested on data obtained from 21 healthy fetuses scanned in-utero at 22-38 weeks gestation. Significantly higher fractional anisotropy values in fiber-rich regions, and the analysis of whole-brain tractography and group structural connectivity, showed the efficacy of the proposed method compared to the analyses based on original data and previously proposed methods. The results of this study show that slice-level motion correction and robust reconstruction is necessary for reliable in-vivo structural connectivity analysis of the fetal brain. Connectivity analysis based on graph theoretic measures show high degree of modularity and clustering, and short average characteristic path lengths indicative of small-worldness property of the fetal brain network. These findings comply with previous findings in newborns and a recent study on fetuses. The proposed algorithm can provide valuable information from DWI of the fetal brain not available in the assessment of the original 2D slices and may be used to more reliably study the developing fetal brain connectome.

Automated template-based brain localization and extraction for fetal brain MRI reconstruction

Most fetal brain MRI reconstruction algorithms rely only on brain tissue-relevant voxels of low-resolution (LR) images to enhance the quality of inter-slice motion correction and image reconstruction. Consequently the fetal brain needs to be localized and extracted as a first step, which is usually a laborious and time consuming manual or semi-automatic task. We have proposed in this work to use age-matched template images as prior knowledge to automatize brain localization and extraction. This has been achieved through a novel automatic brain localization and extraction method based on robust template-to-slice block matching and deformable slice-to-template registration. Our template-based approach has also enabled the reconstruction of fetal brain images in standard radiological anatomical planes in a common coordinate space. We have integrated this approach into our new reconstruction pipeline that involves intensity normalization, inter-slice motion correction, and super-resolution (SR) reconstruction. To this end we have adopted a novel approach based on projection of every slice of the LR brain masks into the template space using a fusion strategy. This has enabled the refinement of brain masks in the LR images at each motion correction iteration. The overall brain localization and extraction algorithm has shown to produce brain masks that are very close to manually drawn brain masks, showing an average Dice overlap measure of 94.5%. We have also demonstrated that adopting a slice-to-template registration and propagation of the brain mask slice-by-slice leads to a significant improvement in brain extraction performance compared to global rigid brain extraction and consequently in the quality of the final reconstructed images. Ratings performed by two expert observers show that the proposed pipeline can achieve similar reconstruction quality to reference reconstruction based on manual slice-by-slice brain extraction. The proposed brain mask refinement and reconstruction method has shown to provide promising results in automatic fetal brain MRI segmentation and volumetry in 26 fetuses with gestational age range of 23 to 38 weeks.

Heavy Prenatal Alcohol Exposure is Related to Smaller Corpus Callosum in Newborn MRI Scans

BACKGROUND

Magnetic resonance imaging (MRI) studies have consistently demonstrated disproportionately smaller corpus callosa in individuals with a history of prenatal alcohol exposure (PAE) but have not previously examined the feasibility of detecting this effect in infants. Tissue segmentation of the newborn brain is challenging because analysis techniques developed for the adult brain are not directly transferable, and segmentation for cerebral morphometry is difficult in neonates, due to the latter's incomplete myelination. This study is the first to use volumetric structural MRI to investigate PAE effects in newborns using manual tracing and to examine the cross-sectional area of the corpus callosum (CC).

METHODS

Forty-three nonsedated infants born to 32 Cape Coloured heavy drinkers and 11 controls recruited prospectively during pregnancy were scanned using a custom-designed birdcage coil for infants, which increases signal-to-noise ratio almost 2-fold compared to the standard head coil. Alcohol use was ascertained prospectively during pregnancy, and fetal alcohol spectrum disorders diagnosis was conducted by expert dysmorphologists. Data were acquired using a multi-echo FLASH protocol adapted for newborns, and a knowledge-based procedure was used to hand-segment the neonatal brains.

RESULTS

CC was disproportionately smaller in alcohol-exposed neonates than controls after controlling for intracranial volume. By contrast, CC area was unrelated to infant sex, gestational age, age at scan, or maternal smoking, marijuana, or methamphetamine use during pregnancy.

CONCLUSIONS

Given that midline craniofacial anomalies have been recognized as a hallmark of fetal alcohol syndrome in humans and animal models since this syndrome was first identified, the CC deficit identified here in newborns may support early identification of a range of midline structural impairments. Smaller CC during the newborn period may provide an early indicator of fetal alcohol-related cognitive deficits that have been linked to this critically important brain structure in childhood and adolescence.

Longitudinal multiple sclerosis lesion segmentation: Resource and challenge

In conjunction with the ISBI 2015 conference, we organized a longitudinal lesion segmentation challenge providing training and test data to registered participants. The training data consisted of five subjects with a mean of 4.4 time-points, and test data of fourteen subjects with a mean of 4.4 time-points. All 82 data sets had the white matter lesions associated with multiple sclerosis delineated by two human expert raters. Eleven teams submitted results using state-of-the-art lesion segmentation algorithms to the challenge, with ten teams presenting their results at the conference. We present a quantitative evaluation comparing the consistency of the two raters as well as exploring the performance of the eleven submitted results in addition to three other lesion segmentation algorithms. The challenge presented three unique opportunities: (1) the sharing of a rich data set; (2) collaboration and comparison of the various avenues of research being pursued in the community; and (3) a review and refinement of the evaluation metrics currently in use. We report on the performance of the challenge participants, as well as the construction and evaluation of a consensus delineation. The image data and manual delineations will continue to be available for download, through an evaluation website2 as a resource for future researchers in the area. This data resource provides a platform to compare existing methods in a fair and consistent manner to each other and multiple manual raters.

Active delineation of Meyer's loop using oriented priors through MAGNEtic tractography (MAGNET)

Streamline tractography algorithms infer connectivity from diffusion MRI (dMRI) by following diffusion directions which are similarly aligned between neighboring voxels. However, not all white matter (WM) fascicles are organized in this manner. For example, Meyer's loop is a highly curved portion of the optic radiation (OR) that exhibits a narrow turn, kissing and crossing pathways, and changes in fascicle dispersion. From a neurosurgical perspective, damage to Meyer's loop carries a potential risk of inducing vision deficits to the patient, especially during temporal lobe resection surgery. To prevent such impairment, achieving an accurate delineation of Meyer's loop with tractography is thus of utmost importance. However, current algorithms tend to under-estimate the full extent of Meyer's loop, mainly attributed to the aforementioned rule for connectivity which requires a direction to be chosen across a field of orientations. In this article, it was demonstrated that MAGNEtic Tractography (MAGNET) can benefit Meyer's loop delineation by incorporating anatomical knowledge of the expected fiber orientation to overcome local ambiguities. A new ROI-mechanism was proposed which supplies additional information to streamline reconstruction algorithms by the means of oriented priors. Their results showed that MAGNET can accurately generate Meyer's loop in all of our 15 child subjects (8 males; mean age 10.2 years ± 3.1). It effectively improved streamline coverage when compared with deterministic tractography, and significantly reduced the distance between the anterior-most portion of Meyer's loop and the temporal pole by 16.7 mm on average, a crucial landmark used for preoperative planning of temporal lobe surgery. Hum Brain Mapp 38:509-527, 2017. © 2016 Wiley Periodicals, Inc.

Dynamic Electrical Source Imaging (DESI) of Seizures and Interictal Epileptic Discharges Without Ensemble Averaging

We propose an algorithm for electrical source imaging of epileptic discharges that takes a data-driven approach to regularizing the dynamics of solutions. The method is based on linear system identification on short time segments, combined with a classical inverse solution approach. Whereas ensemble averaging of segments or epochs discards inter-segment variations by averaging across them, our approach explicitly models them. Indeed, it may even be possible to avoid the need for the time-consuming process of marking epochs containing discharges altogether. We demonstrate that this approach can produce both stable and accurate inverse solutions in experiments using simulated data and real data from epilepsy patients. In an illustrative example, we show that we are able to image propagation using this approach. We show that when applied to imaging seizure data, our approach reproducibly localized frequent seizure activity to within the margins of surgeries that led to patients' seizure freedom. The same approach could be used in the planning of epilepsy surgeries, as a way to localize potentially epileptogenic tissue that should be resected.

Evaluation of motion and its effect on brain magnetic resonance image quality in children

BACKGROUND

Motion artifacts pose significant problems for the acquisition of MR images in pediatric populations.

OBJECTIVE

To evaluate temporal motion metrics in MRI scanners and their effect on image quality in pediatric populations in neuroimaging studies.

MATERIALS AND METHODS

We report results from a large pediatric brain imaging study that shows the effect of motion on MRI quality. We measured motion metrics in 82 pediatric patients, mean age 13.4 years, in a T1-weighted brain MRI scan. As a result of technical difficulties, 5 scans were not included in the subsequent analyses. A radiologist graded the images using a 4-point scale ranging from clinically non-diagnostic because of motion artifacts to no motion artifacts. We used these grades to correlate motion parameters such as maximum motion, mean displacement from a reference point, and motion-free time with image quality.

RESULTS

Our results show that both motion-free time (as a ratio of total scan time) and average displacement from a position at a fixed time (when the center of k-space was acquired) were highly correlated with image quality, whereas maximum displacement was not as good a predictor. Among the 77 patients whose motion was measured successfully, 17 had average displacements of greater than 0.5 mm, and 11 of those (14.3%) resulted in non-diagnostic images. Similarly, 14 patients (18.2%) had less than 90% motion-free time, which also resulted in non-diagnostic images.

CONCLUSION

We report results from a large pediatric study to show how children and young adults move in the MRI scanner and the effect that this motion has on image quality. The results will help the motion-correction community in better understanding motion patterns in pediatric populations and how these patterns affect MR image quality.

A new neonatal cortical and subcortical brain atlas: the Melbourne Children's Regional Infant Brain (M-CRIB) atlas

Investigating neonatal brain structure and function can offer valuable insights into behaviour and cognition in healthy and clinical populations; both at term age, and longitudinally in comparison with later time points. Parcellated brain atlases for adult populations are readily available, however warping infant data to adult template space is not ideal due to morphological and tissue differences between these groups. Several parcellated neonatal atlases have been developed, although there remains strong demand for manually parcellated ground truth data with detailed cortical definition. Additionally, compatibility with existing adult atlases is favourable for use in longitudinal investigations. We aimed to address these needs by replicating the widely-used Desikan-Killiany (2006) adult cortical atlas in neonates. We also aimed to extend brain coverage by complementing this cortical scheme with basal ganglia, thalamus, cerebellum and other subcortical segmentations. Thus, we have manually parcellated these areas volumetrically using high-resolution neonatal T₂-weighted MRI scans, and initial automated and manually edited tissue classification, providing 100 regions in all. Linear and nonlinear T₂-weighted structural templates were also generated. In this paper we provide manual parcellation protocols, and present the parcellated probability maps and structural templates together as the Melbourne Children's Regional Infant Brain (M-CRIB) atlas.

Motion-Robust Diffusion-Weighted Brain MRI Reconstruction Through Slice-Level Registration-Based Motion Tracking

This work proposes a novel approach for motion-robust diffusion-weighted (DW) brain MRI reconstruction through tracking temporal head motion using slice-to-volume registration. The slice-level motion is estimated through a filtering approach that allows tracking the head motion during the scan and correcting for out-of-plane inconsistency in the acquired images. Diffusion-sensitized image slices are registered to a base volume sequentially over time in the acquisition order where an outlier-robust Kalman filter, coupled with slice-to-volume registration, estimates head motion parameters. Diffusion gradient directions are corrected for the aligned DWI slices based on the computed rotation parameters and the diffusion tensors are directly estimated from the corrected data at each voxel using weighted linear least squares. The method was evaluated in DWI scans of adult volunteers who deliberately moved during scans as well as clinical DWI of 28 neonates and children with different types of motion. Experimental results showed marked improvements in DWI reconstruction using the proposed method compared to the state-of-the-art DWI analysis based on volume-to-volume registration. This approach can be readily used to retrieve information from motion-corrupted DW imaging data.