Lesion size and shape in central vein sign assessment for multiple sclerosis diagnosis: An in vivo and postmortem MRI study.

BACKGROUND: The "central vein sign" (CVS), a linear hypointensity on T2*-weighted imaging corresponding to a central vein/venule, is associated with multiple sclerosis (MS) lesions. The effect of lesion-size exclusion criteria on MS diagnostic accuracy has not been extensively studied. OBJECTIVE: Investigate the optimal lesion-size exclusion criteria for CVS use in MS diagnosis. METHODS: Cross-sectional study of 163 MS and 51 non-MS, and radiological/histopathological correlation of 5 MS and 1 control autopsy cases. The effects of lesion-size exclusion on MS diagnosis using the CVS, and intralesional vein detection on histopathology were evaluated. RESULTS: CVS+ lesions were larger compared to CVS- lesions, with effect modification by MS diagnosis (mean difference +7.7 mm3, p = 0.004). CVS percentage-based criteria with no lesion-size exclusion showed the highest diagnostic accuracy in differentiating MS cases. However, a simple count of three or more CVS+ lesions greater than 3.5 mm is highly accurate and can be rapidly implemented (sensitivity 93%; specificity 88%). On magnetic resonance imaging (MRI)-histopathological correlation, the CVS had high specificity for identifying intralesional veins (0/7 false positives). CONCLUSION: Lesion-size measures add important information when using CVS+ lesion counts for MS diagnosis. The CVS is a specific biomarker corresponding to intralesional veins on histopathology.

Gradient nonlinearity effects on upper cervical spinal cord area measurement from 3D T1 -weighted brain MRI acquisitions.

PURPOSE: To explore (i) the variability of upper cervical cord area (UCCA) measurements from volumetric brain 3D T1 -weighted scans related to gradient nonlinearity (GNL) and subject positioning; (ii) the effect of vendor-implemented GNL corrections; and (iii) easily applicable methods that can be used to retrospectively correct data. METHODS: A multiple sclerosis patient was scanned at seven sites using 3T MRI scanners with the same 3D T1 -weighted protocol without GNL-distortion correction. Two healthy subjects and a phantom were additionally scanned at a single site with varying table positions. The 2D and 3D vendor-implemented GNL-correction algorithms and retrospective methods based on (i) phantom data fit, (ii) normalization with C2 vertebral body diameters, and (iii) the Jacobian determinant of nonlinear registrations to a template were tested. RESULTS: Depending on the positioning of the subject, GNL introduced up to 15% variability in UCCA measurements from volumetric brain T1 -weighted scans when no distortion corrections were used. The 3D vendor-implemented correction methods and the three proposed methods reduced this variability to less than 3%. CONCLUSIONS: Our results raise awareness of the significant impact that GNL can have on quantitative UCCA studies, and point the way to prospectively and retrospectively managing GNL distortions in a variety of settings, including clinical environments. Magn Reson Med 79:1595-1601, 2018. © 2017 International Society for Magnetic Resonance in Medicine.

Robust skull stripping using multiple MR image contrasts insensitive to pathology.

Automatic skull-stripping or brain extraction of magnetic resonance (MR) images is often a fundamental step in many neuroimage processing pipelines. The accuracy of subsequent image processing relies on the accuracy of the skull-stripping. Although many automated stripping methods have been proposed in the past, it is still an active area of research particularly in the context of brain pathology. Most stripping methods are validated on T1-w MR images of normal brains, especially because high resolution T1-w sequences are widely acquired and ground truth manual brain mask segmentations are publicly available for normal brains. However, different MR acquisition protocols can provide complementary information about the brain tissues, which can be exploited for better distinction between brain, cerebrospinal fluid, and unwanted tissues such as skull, dura, marrow, or fat. This is especially true in the presence of pathology, where hemorrhages or other types of lesions can have similar intensities as skull in a T1-w image. In this paper, we propose a sparse patch based Multi-cONtrast brain STRipping method (MONSTR),2 where non-local patch information from one or more atlases, which contain multiple MR sequences and reference delineations of brain masks, are combined to generate a target brain mask. We compared MONSTR with four state-of-the-art, publicly available methods: BEaST, SPECTRE, ROBEX, and OptiBET. We evaluated the performance of these methods on 6 datasets consisting of both healthy subjects and patients with various pathologies. Three datasets (ADNI, MRBrainS, NAMIC) are publicly available, consisting of 44 healthy volunteers and 10 patients with schizophrenia. Other three in-house datasets, comprising 87 subjects in total, consisted of patients with mild to severe traumatic brain injury, brain tumors, and various movement disorders. A combination of T1-w, T2-w were used to skull-strip these datasets. We show significant improvement in stripping over the competing methods on both healthy and pathological brains. We also show that our multi-contrast framework is robust and maintains accurate performance across different types of acquisitions and scanners, even when using normal brains as atlases to strip pathological brains, demonstrating that our algorithm is applicable even when reference segmentations of pathological brains are not available to be used as atlases.

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.

Optical coherence tomography reflects brain atrophy in multiple sclerosis: A four-year study.

OBJECTIVE: The aim of this work was to determine whether atrophy of specific retinal layers and brain substructures are associated over time, in order to further validate the utility of optical coherence tomography (OCT) as an indicator of neuronal tissue damage in patients with multiple sclerosis (MS). METHODS: Cirrus high-definition OCT (including automated macular segmentation) was performed in 107 MS patients biannually (median follow-up: 46 months). Three-Tesla magnetic resonance imaging brain scans (including brain-substructure volumetrics) were performed annually. Individual-specific rates of change in retinal and brain measures (estimated with linear regression) were correlated, adjusting for age, sex, disease duration, and optic neuritis (ON) history. RESULTS: Rates of ganglion cell + inner plexiform layer (GCIP) and whole-brain (r = 0.45; p < 0.001), gray matter (GM; r = 0.37; p < 0.001), white matter (WM; r = 0.28; p = 0.007), and thalamic (r = 0.38; p < 0.001) atrophy were associated. GCIP and whole-brain (as well as GM and WM) atrophy rates were more strongly associated in progressive MS (r = 0.67; p < 0.001) than relapsing-remitting MS (RRMS; r = 0.33; p = 0.007). However, correlation between rates of GCIP and whole-brain (and additionally GM and WM) atrophy in RRMS increased incrementally with step-wise refinement to exclude ON effects; excluding eyes and then patients (to account for a phenotype effect), the correlation increased to 0.45 and 0.60, respectively, consistent with effect modification. In RRMS, lesion accumulation rate was associated with GCIP (r = -0.30; p = 0.02) and inner nuclear layer (r = -0.25; p = 0.04) atrophy rates. INTERPRETATION: Over time GCIP atrophy appears to mirror whole-brain, and particularly GM, atrophy, especially in progressive MS, thereby reflecting underlying disease progression. Our findings support OCT for clinical monitoring and as an outcome in investigative trials.

Thalamic lesions in multiple sclerosis by 7T MRI: Clinical implications and relationship to cortical pathology.

OBJECTIVE: Pathology in both cortex and deep gray matter contribute to disability in multiple sclerosis (MS). We used the increased signal-to-noise ratio of 7-tesla (7T) MRI to visualize small lesions within the thalamus and to relate this to clinical information and cortical lesions. METHODS: We obtained 7T MRI scans on 34 MS cases and 15 healthy volunteers. Thalamic lesion number and volume were related to demographic data, clinical disability measures, and lesions in cortical gray matter. RESULTS: Thalamic lesions were found in 24/34 of MS cases. Two lesion subtypes were noted: discrete, ovoid lesions, and more diffuse lesional areas lining the periventricular surface. The number of thalamic lesions was greater in progressive MS compared to relapsing-remitting (mean ±SD, 10.7 ±0.7 vs. 3.0 ±0.7, respectively, p < 0.001). Thalamic lesion burden (count and volume) correlated with EDSS score and measures of cortical lesion burden, but not with white matter lesion burden or white matter volume. CONCLUSIONS: Using 7T MRI allows identification of thalamic lesions in MS, which are associated with disability, progressive disease, and cortical lesions. Thalamic lesion analysis may be a simpler, more rapid estimate of overall gray matter lesion burden in MS.

Constitutive and conditional deletion reveals distinct phenotypes driven by developmental versus neurotransmitter actions of the neuropeptide PACAP.

Neuropeptides may exert trophic effects during development, and then neurotransmitter roles in the developed nervous system. One way to associate peptide-deficiency phenotypes with either role is first to assess potential phenotypes in so-called constitutive knockout mice, and then proceed to specify, regionally and temporally, where and when neuropeptide expression is required to prevent these phenotypes. We have previously demonstrated that the well-known constellation of behavioral and metabolic phenotypes associated with constitutive pituitary adenylate cyclase-activating peptide (PACAP) knockout mice are accompanied by transcriptomic alterations of two types: those that distinguish the PACAP-null phenotype from wild-type (WT) in otherwise quiescent mice (cPRGs), and gene induction that occurs in response to acute environmental perturbation in WT mice that do not occur in knockout mice (aPRGs). Comparing constitutive PACAP knockout mice to a variety of temporally and regionally specific PACAP knockouts, we show that the prominent hyperlocomotor phenotype is a consequence of early loss of PACAP expression, is associated with Fos overexpression in hippocampus and basal ganglia, and that a thermoregulatory effect previously shown to be mediated by PACAP-expressing neurons of medial preoptic hypothalamus is independent of PACAP expression in those neurons in adult mice. In contrast, PACAP dependence of weight loss/hypophagia triggered by restraint stress, seen in constitutive PACAP knockout mice, is phenocopied in mice in which PACAP is deleted after neuronal differentiation. Our results imply that PACAP has a prominent role as a trophic factor early in development determining global central nervous system characteristics, and in addition a second, discrete set of functions as a neurotransmitter in the fully developed nervous system that support physiological and psychological responses to stress.

Dissociable control of motivation and reinforcement by distinct ventral striatal dopamine receptors.

Dopamine release in striatal circuits, including the nucleus accumbens (NAc), tracks separable features of reward such as motivation and reinforcement. However, the cellular and circuit mechanisms by which dopamine receptors transform dopamine release into distinct constructs of reward remain unclear. Here, we show that dopamine D3 receptor (D3R) signaling in the NAc drives motivated behavior by regulating local NAc microcircuits. Furthermore, D3Rs co-express with dopamine D1 receptors (D1Rs), which regulate reinforcement, but not motivation. Paralleling dissociable roles in reward function, we report non-overlapping physiological actions of D3R and D1R signaling in NAc neurons. Our results establish a novel cellular framework wherein dopamine signaling within the same NAc cell type is physiologically compartmentalized via actions on distinct dopamine receptors. This structural and functional organization provides neurons in a limbic circuit with the unique ability to orchestrate dissociable aspects of reward-related behaviors that are relevant to the etiology of neuropsychiatric disorders.

In vivo MRI is sensitive to remyelination in a nonhuman primate model of multiple sclerosis.

Remyelination is crucial to recover from inflammatory demyelination in multiple sclerosis (MS). Investigating remyelination in vivo using magnetic resonance imaging (MRI) is difficult in MS, where collecting serial short-interval scans is challenging. Using experimental autoimmune encephalomyelitis (EAE) in common marmosets, a model of MS that recapitulates focal cerebral inflammatory demyelinating lesions, we investigated whether MRI is sensitive to, and can characterize, remyelination. In six animals followed with multisequence 7 T MRI, 31 focal lesions, predicted to be demyelinated or remyelinated based on signal intensity on proton density-weighted images, were subsequently assessed with histopathology. Remyelination occurred in four of six marmosets and 45% of lesions. Radiological-pathological comparison showed that MRI had high statistical sensitivity (100%) and specificity (90%) for detecting remyelination. This study demonstrates the prevalence of spontaneous remyelination in marmoset EAE and the ability of in vivo MRI to detect it, with implications for preclinical testing of pro-remyelinating agents.

Central Vein Sign Profile of Newly Developing Lesions in Multiple Sclerosis: A 3-Year Longitudinal Study.

BACKGROUND AND OBJECTIVES: The central vein sign (CVS), a central linear hypointensity within lesions on T2*-weighted imaging, has been established as a sensitive and specific biomarker for the diagnosis of multiple sclerosis (MS). However, the CVS has not yet been comprehensively studied in newly developing MS lesions. We aimed to identify the CVS profiles of new white matter lesions in patients with MS followed over time and investigate demographic and clinical risk factors associated with new CVS+ or CVS- lesion development. METHODS: In this retrospective longitudinal cohort study, adults from the NIH MS Natural History Study were considered for inclusion. Participants with new T2 or enhancing lesions were identified through review of the radiology report and/or longitudinal subtraction imaging. Each new lesion was evaluated for the CVS. Clinical characteristics were identified through chart review. RESULTS: A total of 153 adults (95 relapsing-remitting MS, 27 secondary progressive MS, 16 primary progressive MS, 5 clinically isolated syndrome, and 10 healthy; 67% female) were included. Of this cohort, 96 had at least 1 new T2 or contrast-enhancing lesion during median 3.1 years (Q1-Q3: 0.7-6.3) of follow-up; lesions eligible for CVS evaluation were found in 62 (65%). Of 233 new CVS-eligible lesions, 159 (68%) were CVS+, with 30 (48%) individuals having only CVS+, 12 (19%) only CVS-, and 20 (32%) both CVS+ and CVS- lesions. In gadolinium-enhancing (Gd+) lesions, the CVS+ percentage increased from 102/152 (67%) at the first time point where the lesion was observed, to 92/114 (82%) after a median follow-up of 2.8 years. Younger age (OR = 0.5 per 10-year increase, 95% CI = 0.3-0.8) and higher CVS+ percentage at baseline (OR = 1.4 per 10% increase, 95% CI = 1.1-1.9) were associated with increased likelihood of new CVS+ lesion development. DISCUSSION: In a cohort of adults with MS followed over a median duration of 3 years, most newly developing T2 or enhancing lesions were CVS+ (68%), and nearly half (48%) developed new CVS+ lesions only. Importantly, the effects of edema and T2 signal changes can obscure small veins in Gd+ lesions; therefore, caution and follow-up is necessary when determining their CVS status. TRIAL REGISTRATION INFORMATION: Clinical trial registration number NCT00001248. CLASSIFICATION OF EVIDENCE: This study provides Class III evidence that younger age and higher CVS+ percentage at baseline are associated with new CVS+ lesion development.

Pupal behavior emerges from unstructured muscle activity in response to neuromodulation in Drosophila.

Identifying neural substrates of behavior requires defining actions in terms that map onto brain activity. Brain and muscle activity naturally correlate via the output of motor neurons, but apart from simple movements it has been difficult to define behavior in terms of muscle contractions. By mapping the musculature of the pupal fruit fly and comprehensively imaging muscle activation at single-cell resolution, we here describe a multiphasic behavioral sequence in Drosophila. Our characterization identifies a previously undescribed behavioral phase and permits extraction of major movements by a convolutional neural network. We deconstruct movements into a syllabary of co-active muscles and identify specific syllables that are sensitive to neuromodulatory manipulations. We find that muscle activity shows considerable variability, with sequential increases in stereotypy dependent upon neuromodulation. Our work provides a platform for studying whole-animal behavior, quantifying its variability across multiple spatiotemporal scales, and analyzing its neuromodulatory regulation at cellular resolution.

How do we find out how the brain works? One way is to use imaging techniques to visualise an animal's brain in action as it performs simple behaviours: as the animal moves, parts of its brain light up under the microscope. For laboratory animals like fruit flies, which have relatively small brains, this lets us observe their brain activity right down to the level of individual brain cells. The brain directs movements via collective activity of the body's muscles. Our ability to track the activity of individual muscles is, however, more limited than our ability to observe single brain cells: even modern imaging technology still cannot monitor the activity of all the muscle cells in an animal's body as it moves about. Yet this is precisely the information that scientists need to fully understand how the brain generates behaviour. Fruit flies perform specific behaviours at certain stages of their life cycle. When the fly pupa begins to metamorphose into an adult insect, it performs a fixed sequence of movements involving a set number of muscles, which is called the pupal ecdysis sequence. This initial movement sequence and the rest of metamorphosis both occur within the confines of the pupal case, which is a small, hardened shell surrounding the whole animal. Elliott et al. set out to determine if the fruit fly pupa's ecdysis sequence could be used as a kind of model, to describe a simple behaviour at the level of individual muscles. Imaging experiments used fly pupae that were genetically engineered to produce an activity-dependent fluorescent protein in their muscle cells. Pupal cases were treated with a chemical to make them transparent, allowing easy observation of their visually 'labelled' muscles. This yielded a near-complete record of muscle activity during metamorphosis. Initially, individual muscles became active in small groups. The groups then synchronised with each other over the different regions of the pupa's body to form distinct movements, much as syllables join to form words. This synchronisation was key to progression through metamorphosis and was co-ordinated at each step by specialised nerve cells that produce or respond to specific hormones. These results reveal how the brain might direct muscle activity to produce movement patterns. In the future, Elliott et al. hope to compare data on muscle activity with comprehensive records of brain cell activity, to shed new light on how the brain, muscles, and other factors work together to control behaviour.

Distributed deep learning across multisite datasets for generalized CT hemorrhage segmentation.

PURPOSE: As deep neural networks achieve more success in the wide field of computer vision, greater emphasis is being placed on the generalizations of these models for production deployment. With sufficiently large training datasets, models can typically avoid overfitting their data; however, for medical imaging it is often difficult to obtain enough data from a single site. Sharing data between institutions is also frequently nonviable or prohibited due to security measures and research compliance constraints, enforced to guard protected health information (PHI) and patient anonymity. METHODS: In this paper, we implement cyclic weight transfer with independent datasets from multiple geographically disparate sites without compromising PHI. We compare results between single-site learning (SSL) and multisite learning (MSL) models on testing data drawn from each of the training sites as well as two other institutions. RESULTS: The MSL model attains an average dice similarity coefficient (DSC) of 0.690 on the holdout institution datasets with a volume correlation of 0.914, respectively corresponding to a 7% and 5% statistically significant improvement over the average of both SSL models, which attained an average DSC of 0.646 and average correlation of 0.871. CONCLUSIONS: We show that a neural network can be efficiently trained on data from two physically remote sites without consolidating patient data to a single location. The resulting network improves model generalization and achieves higher average DSCs on external datasets than neural networks trained on data from a single source.

Extracting 2D weak labels from volume labels using multiple instance learning in CT hemorrhage detection.

Multiple instance learning (MIL) is a supervised learning methodology that aims to allow models to learn instance class labels from bag class labels, where a bag is defined to contain multiple instances. MIL is gaining traction for learning from weak labels but has not been widely applied to 3D medical imaging. MIL is well-suited to clinical CT acquisitions since (1) the highly anisotropic voxels hinder application of traditional 3D networks and (2) patch-based networks have limited ability to learn whole volume labels. In this work, we apply MIL with a deep convolutional neural network to identify whether clinical CT head image volumes possess one or more large hemorrhages (> 20cm3), resulting in a learned 2D model without the need for 2D slice annotations. Individual image volumes are considered separate bags, and the slices in each volume are instances. Such a framework sets the stage for incorporating information obtained in clinical reports to help train a 2D segmentation approach. Within this context, we evaluate the data requirements to enable generalization of MIL by varying the amount of training data. Our results show that a training size of at least 400 patient image volumes was needed to achieve accurate per-slice hemorrhage detection. Over a five-fold cross-validation, the leading model, which made use of the maximum number of training volumes, had an average true positive rate of 98.10%, an average true negative rate of 99.36%, and an average precision of 0.9698. The models have been made available along with source code1 to enabled continued exploration and adaption of MIL in CT neuroimaging.

Progressive multifocal leukoencephalopathy lesion and brain parenchymal segmentation from MRI using serial deep convolutional neural networks.

Progressive multifocal leukoencephalopathy (PML) is a rare opportunistic brain infection caused by the JC virus and associated with substantial morbidity and mortality. Accurate MRI assessment of PML lesion burden and brain parenchymal atrophy is of decisive value in monitoring the disease course and response to therapy. However, there are currently no validated automatic methods for quantification of PML lesion burden or associated parenchymal volume loss. Furthermore, manual brain or lesion delineations can be tedious, require the use of valuable time resources by radiologists or trained experts, and are often subjective. In this work, we introduce JCnet (named after the causative viral agent), an end-to-end, fully automated method for brain parenchymal and lesion segmentation in PML using consecutive 3D patch-based convolutional neural networks. The network architecture consists of multi-view feature pyramid networks with hierarchical residual learning blocks containing embedded batch normalization and nonlinear activation functions. The feature maps across the bottom-up and top-down pathways of the feature pyramids are merged, and an output probability membership generated through convolutional pathways, thus rendering the method fully convolutional. Our results show that this approach outperforms and improves longitudinal consistency compared to conventional, state-of-the-art methods of healthy brain and multiple sclerosis lesion segmentation, utilized here as comparators given the lack of available methods validated for use in PML. The ability to produce robust and accurate automated measures of brain atrophy and lesion segmentation in PML is not only valuable clinically but holds promise toward including standardized quantitative MRI measures in clinical trials of targeted therapies. Code is available at: https://github.com/omarallouz/JCnet.

Evaluating White Matter Lesion Segmentations with Refined Sørensen-Dice Analysis.

The Sørensen-Dice index (SDI) is a widely used measure for evaluating medical image segmentation algorithms. It offers a standardized measure of segmentation accuracy which has proven useful. However, it offers diminishing insight when the number of objects is unknown, such as in white matter lesion segmentation of multiple sclerosis (MS) patients. We present a refinement for finer grained parsing of SDI results in situations where the number of objects is unknown. We explore these ideas with two case studies showing what can be learned from our two presented studies. Our first study explores an inter-rater comparison, showing that smaller lesions cannot be reliably identified. In our second case study, we demonstrate fusing multiple MS lesion segmentation algorithms based on the insights into the algorithms provided by our analysis to generate a segmentation that exhibits improved performance. This work demonstrates the wealth of information that can be learned from refined analysis of medical image segmentations.

Distributed deep learning for robust multi-site segmentation of CT imaging after traumatic brain injury.

Machine learning models are becoming commonplace in the domain of medical imaging, and with these methods comes an ever-increasing need for more data. However, to preserve patient anonymity it is frequently impractical or prohibited to transfer protected health information (PHI) between institutions. Additionally, due to the nature of some studies, there may not be a large public dataset available on which to train models. To address this conundrum, we analyze the efficacy of transferring the model itself in lieu of data between different sites. By doing so we accomplish two goals: 1) the model gains access to training on a larger dataset that it could not normally obtain and 2) the model better generalizes, having trained on data from separate locations. In this paper, we implement multi-site learning with disparate datasets from the National Institutes of Health (NIH) and Vanderbilt University Medical Center (VUMC) without compromising PHI. Three neural networks are trained to convergence on a computed tomography (CT) brain hematoma segmentation task: one only with NIH data, one only with VUMC data, and one multi-site model alternating between NIH and VUMC data. Resultant lesion masks with the multi-site model attain an average Dice similarity coefficient of 0.64 and the automatically segmented hematoma volumes correlate to those done manually with a Pearson correlation coefficient of 0.87, corresponding to an 8% and 5% improvement, respectively, over the single-site model counterparts.

Deformable model reconstruction of the subarachnoid space.

The subarachnoid space is a layer in the meninges that surrounds the brain and is filled with trabeculae and cerebrospinal fluid. Quantifying the volume and thickness of the subarachnoid space is of interest in order to study the pathogenesis of neurodegenerative diseases and compare with healthy subjects. We present an automatic method to reconstruct the subarachnoid space with subvoxel accuracy using a nested deformable model. The method initializes the deformable model using the convex hull of the union of the outer surfaces of the cerebrum, cerebellum and brainstem. A region force is derived from the subject's Tl-weighted and T2-weighted MRI to drive the deformable model to the outer surface of the subarachnoid space. The proposed method is compared to a semi-automatic delineation from the subject's T2-weighted MRI and an existing multi-atlas-based method. A small pilot study comparing the volume and thickness measurements in a set of age-matched subjects with normal pressure hydrocephalus and healthy controls is presented to show the efficacy of the proposed method.

Joint Intensity Fusion Image Synthesis Applied to Multiple Sclerosis Lesion Segmentation.

We propose a new approach to Multiple Sclerosis lesion segmentation that utilizes synthesized images. A new method of image synthesis is considered: joint intensity fusion (JIF). JIF synthesizes an image from a library of deformably registered and intensity normalized atlases. Each location in the synthesized image is a weighted average of the registered atlases; atlas weights vary spatially. The weights are determined using the joint label fusion (JLF) framework. The primary methodological contribution is the application of JLF to MRI signal directly rather than labels. Synthesized images are then used as additional features in a lesion segmentation task using the OASIS classifier, a logistic regression model on intensities from multiple modalities. The addition of JIF synthesized images improved the Dice-Sorensen coefficient (relative to manually drawn gold standards) of lesion segmentations over the standard model segmentations by 0.0462 ± 0.0050 (mean ± standard deviation) at optimal threshold over all subjects and 10 separate training/testing folds.

Dice Overlap Measures for Objects of Unknown Number: Application to Lesion Segmentation.

The Dice overlap ratio is commonly used to evaluate the performance of image segmentation algorithms. While Dice overlap is very useful as a standardized quantitative measure of segmentation accuracy in many applications, it offers a very limited picture of segmentation quality in complex segmentation tasks where the number of target objects is not known a priori, such as the segmentation of white matter lesions or lung nodules. While Dice overlap can still be used in these applications, segmentation algorithms may perform quite differently in ways not reflected by differences in their Dice score. Here we propose a new set of evaluation techniques that offer new insights into the behavior of segmentation algorithms. We illustrate these techniques with a case study comparing two popular multiple sclerosis (MS) lesion segmentation algorithms: OASIS and LesionTOADS.

Random forest regression for magnetic resonance image synthesis.

By choosing different pulse sequences and their parameters, magnetic resonance imaging (MRI) can generate a large variety of tissue contrasts. This very flexibility, however, can yield inconsistencies with MRI acquisitions across datasets or scanning sessions that can in turn cause inconsistent automated image analysis. Although image synthesis of MR images has been shown to be helpful in addressing this problem, an inability to synthesize both T2-weighted brain images that include the skull and FLuid Attenuated Inversion Recovery (FLAIR) images has been reported. The method described herein, called REPLICA, addresses these limitations. REPLICA is a supervised random forest image synthesis approach that learns a nonlinear regression to predict intensities of alternate tissue contrasts given specific input tissue contrasts. Experimental results include direct image comparisons between synthetic and real images, results from image analysis tasks on both synthetic and real images, and comparison against other state-of-the-art image synthesis methods. REPLICA is computationally fast, and is shown to be comparable to other methods on tasks they are able to perform. Additionally REPLICA has the capability to synthesize both T2-weighted images of the full head and FLAIR images, and perform intensity standardization between different imaging datasets.