Search

The last few years have seen a boom in using generative models to augment real datasets, as synthetic data can effectively model real data distributions and provide privacy-preserving, shareable datasets that can be used to train deep learning models. However, most of these methods are 2D and provide synthetic datasets that come, at most, with categorical annotations. The generation of paired images and segmentation samples that can be used in downstream, supervised segmentation tasks remains fairly uncharted territory. This work proposes a two-stage generative model capable of producing 2D and 3D semantic label maps and corresponding multi-modal images. We use a latent diffusion model for label synthesis and a VAE-GAN for semantic image synthesis. Synthetic datasets provided by this model are shown to work in a wide variety of segmentation tasks, supporting small, real datasets or fully replacing them while maintaining good performance. We also demonstrate its ability to improve downstream performance on out-of-distribution data., Competing Interests: Declaration of competing interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Virginia Fernandez reports financial support, administrative support, and travel were provided by EPSRC Centre for Doctoral Training in Smart Medical Imaging. Tom Vercauteren reports a relationship with Hypervision Surgical that includes: board membership and equity or stocks. If there are other authors, they declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Combining multi-site data can strengthen and uncover trends, but is a task that is marred by the influence of site-specific covariates that can bias the data and, therefore, any downstream analyses. Post-hoc multi-site correction methods exist but have strong assumptions that often do not hold in real-world scenarios. Algorithms should be designed in a way that can account for site-specific effects, such as those that arise from sequence parameter choices, and in instances where generalisation fails, should be able to identify such a failure by means of explicit uncertainty modelling. This body of work showcases such an algorithm that can become robust to the physics of acquisition in the context of segmentation tasks while simultaneously modelling uncertainty. We demonstrate that our method not only generalises to complete holdout datasets, preserving segmentation quality but does so while also accounting for site-specific sequence choices, which also allows it to perform as a harmonisation tool., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023. Published by Elsevier B.V.)

We describe CounterSynth, a conditional generative model of diffeomorphic deformations that induce label-driven, biologically plausible changes in volumetric brain images. The model is intended to synthesise counterfactual training data augmentations for downstream discriminative modelling tasks where fidelity is limited by data imbalance, distributional instability, confounding, or underspecification, and exhibits inequitable performance across distinct subpopulations. Focusing on demographic attributes, we evaluate the quality of synthesised counterfactuals with voxel-based morphometry, classification and regression of the conditioning attributes, and the Fréchet inception distance. Examining downstream discriminative performance in the context of engineered demographic imbalance and confounding, we use UK Biobank and OASIS magnetic resonance imaging data to benchmark CounterSynth augmentation against current solutions to these problems. We achieve state-of-the-art improvements, both in overall fidelity and equity. The source code for CounterSynth is available at https://github.com/guilherme-pombo/CounterSynth., Competing Interests: Declaration of Competing Interest The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: Guilherme Pombo reports financial support was provided by Wellcome Trust. Guilherme Pombo reports financial support was provided by NIHR UCLH Biomedical Research Centre., (Copyright © 2022 The Authors. Published by Elsevier B.V. All rights reserved.)

Accurate regional brain quantitative PET measurements, particularly when using partial volume correction, rely on robust image registration between PET and MR images. We argue here that the precision, and hence the uncertainty, of MR-PET image registration is mainly driven by the registration implementation and the quality of PET images due to their lower resolution and higher noise compared to the structural MR images. We propose a dedicated uncertainty analysis for quantifying the precision of MR-PET registration, centred around the bootstrap resampling of PET list-mode events to generate multiple PET image realisations with different noise (count) levels. The effects of PET image reconstruction parameters, such as the use of attenuation and scatter corrections and different number of iterations, on the precision and accuracy of MR-PET registration were investigated. In addition, the performance of four software packages with their default settings for rigid inter-modality image registration were considered: NiftyReg, Vinci, FSL and SPM. Four distinct PET image distributions made of two early time frames (similar to cortical FDG) and two late frames using two amyloid PET dynamic acquisitions of one amyloid positive and one amyloid negative participants were investigated. For the investigated four PET frames, the biggest impact on the uncertainty was observed between registration software packages (up to 10-fold difference in precision) followed by the reconstruction parameters. On average, the lowest uncertainty for different PET frames and brain regions was observed with SPM and two iterations of fully quantitative image reconstruction. The observed uncertainty for the varying PET count-level (from 5% to 60%) was slightly lower than for the reconstruction parameters. We also observed that the registration uncertainty in quantitative PET analysis depends on amyloid status of the considered PET frames, with increased uncertainty (up to three times) when using post-reconstruction partial volume correction. This analysis is applicable for PET data obtained from either PET/MR or PET/CT scanners., (Copyright © 2021 The Author(s). Published by Elsevier Inc. All rights reserved.)

Down Syndrome is a chromosomal disorder that affects the development of cerebellar cortical lobules. Impaired neurogenesis in the cerebellum varies among different types of neuronal cells and neuronal layers. In this study, we developed an imaging analysis framework that utilizes gadolinium-enhanced ex vivo mouse brain MRI. We extracted the middle Purkinje layer of the mouse cerebellar cortex, enabling the estimation of the volume, thickness, and surface area of the entire cerebellar cortex, the internal granular layer, and the molecular layer in the Tc1 mouse model of Down Syndrome. The morphometric analysis of our method revealed that a larger proportion of the cerebellar thinning in this model of Down Syndrome resided in the inner granule cell layer, while a larger proportion of the surface area shrinkage was in the molecular layer., (Copyright © 2020 The Author(s). Published by Elsevier Inc. All rights reserved.)

Spinal cord atrophy measurements obtained from structural magnetic resonance imaging (MRI) are associated with disability in many neurological diseases and serve as in vivo biomarkers of neurodegeneration. Longitudinal spinal cord atrophy rate is commonly determined from the numerical difference between two volumes (based on 3D surface fitting) or two cross-sectional areas (CSA, based on 2D edge detection) obtained at different time-points. Being an indirect measure, atrophy rates are susceptible to variable segmentation errors at the edge of the spinal cord. To overcome those limitations, we developed a new registration-based pipeline that measures atrophy rates directly. We based our approach on the generalised boundary shift integral (GBSI) method, which registers 2 scans and uses a probabilistic XOR mask over the edge of the spinal cord, thereby measuring atrophy more accurately than segmentation-based techniques. Using a large cohort of longitudinal spinal cord images (610 subjects with multiple sclerosis from a multi-centre trial and 52 healthy controls), we demonstrated that GBSI is a sensitive, quantitative and objective measure of longitudinal spinal cord volume change. The GBSI pipeline is repeatable, reproducible, and provides more precise measurements of longitudinal spinal cord atrophy than segmentation-based methods in longitudinal spinal cord atrophy studies., (Copyright © 2020 The Authors. Published by Elsevier Inc. All rights reserved.)

Objective: Hippocampal sclerosis (HS) is the most common cause of drug-resistant temporal lobe epilepsy, and its accurate detection is important to guide epilepsy surgery. Radiological features of HS include hippocampal volume loss and increased T2 signal, which can both be quantified to help improve detection. In this work, we extend these quantitative methods to generate cross-sectional area and T2 profiles along the hippocampal long axis to improve the localization of hippocampal abnormalities., Methods: T1-weighted and T2 relaxometry data from 69 HS patients (32 left, 32 right, 5 bilateral) and 111 healthy controls were acquired on a 3-T magnetic resonance imaging (MRI) scanner. Automated hippocampal segmentation and T2 relaxometry were performed and used to calculate whole-hippocampal volumes and to estimate quantitative T2 (qT2) values. By generating a group template from the controls, and aligning this so that the hippocampal long axes were along the anterior-posterior axis, we were able to calculate hippocampal cross-sectional area and qT2 by a slicewise method to localize any volume loss or T2 hyperintensity. Individual patient profiles were compared with normative data generated from the healthy controls., Results: Profiling of hippocampal volumetric and qT2 data could be performed automatically and reproducibly. HS patients commonly showed widespread decreases in volume and increases in T2 along the length of the affected hippocampus, and focal changes may also be identified. Patterns of atrophy and T2 increase in the left hippocampus were similar between left, right, and bilateral HS. These profiles have potential to distinguish between sclerosis affecting volume and qT2 in the whole or parts of the hippocampus, and may aid the radiological diagnosis in uncertain cases or cases with subtle or focal abnormalities where standard whole-hippocampal measurements yield normal values., Significance: Hippocampal profiling of volumetry and qT2 values can help spatially localize hippocampal MRI abnormalities and work toward improved sensitivity of subtle focal lesions., (© 2019 The Authors. Epilepsia published by Wiley Periodicals, Inc. on behalf of International League Against Epilepsy.)

Breast cancer image fusion consists of registering and visualizing different sets of a patient synchronized torso and radiological images into a 3D model. Breast spatial interpretation and visualization by the treating physician can be augmented with a patient-specific digital breast model that integrates radiological images. But the absence of a ground truth for a good correlation between surface and radiological information has impaired the development of potential clinical applications. A new image acquisition protocol was designed to acquire breast Magnetic Resonance Imaging (MRI) and 3D surface scan data with surface markers on the patient's breasts and torso. A patient-specific digital breast model integrating the real breast torso and the tumor location was created and validated with a MRI/3D surface scan fusion algorithm in 16 breast cancer patients. This protocol was used to quantify breast shape differences between different modalities, and to measure the target registration error of several variants of the MRI/3D scan fusion algorithm. The fusion of single breasts without the biomechanical model of pose transformation had acceptable registration errors and accurate tumor locations. The performance of the fusion algorithm was not affected by breast volume. Further research and virtual clinical interfaces could lead to fast integration of this fusion technology into clinical practice., Competing Interests: Declaration of competing interest None., (Copyright © 2020 The Authors. Published by Elsevier Ltd.. All rights reserved.)

Introduction: To date, the clinical relevance of comorbid amyloid-β (Aβ) pathology in patients with vascular cognitive disorders (VCD) is largely unknown., Methods: We included 218 VCD patients with available cerebrospinal fluid Aβ 42 levels. Patients were divided into Aβ+ mild-VCD (n = 84), Aβ- mild-VCD (n = 68), Aβ+ major-VCD (n = 31), and Aβ- major-VCD (n = 35). We measured depression with the Geriatric Depression Scale, cognition with a neuropsychological test battery and derived white matter hyperintensities (WMH) and gray matter atrophy from MRI., Results: Aβ- patients showed more depressive symptoms than Aβ+. In the major-VCD group, Aβ- patients performed worse on attention (P = .02) and executive functioning (P = .008) than Aβ+. We found no cognitive differences in patients with mild VCD. In the mild-VCD group, Aβ- patients had more WMH than Aβ+ patients, whereas conversely, in the major-VCD group, Aβ+ patients had more WMH. Atrophy patterns did not differ between Aβ+ and Aβ- VCD group., Discussion: Comorbid Aβ pathology affects the manifestation of VCD, but effects differ by severity of VCD., (© 2019 The Authors. Alzheimer's & Dementia published by Wiley Periodicals, Inc. on behalf of Alzheimer's Association.)

Positron emission tomography/magnetic resonance imaging (PET/MRI) potentially offers several advantages over positron emission tomography/computed tomography (PET/CT), for example, no CT radiation dose and soft tissue images from MR acquired at the same time as the PET. However, obtaining accurate linear attenuation correction (LAC) factors for the lung remains difficult in PET/MRI. LACs depend on electron density and in the lung, these vary significantly both within an individual and from person to person. Current commercial practice is to use a single-valued population-based lung LAC, and better estimation is needed to improve quantification. Given the under-appreciation of lung attenuation estimation as an issue, the inaccuracy of PET quantification due to the use of single-valued lung LACs, the unique challenges of lung estimation, and the emerging status of PET/MRI scanners in lung disease, a review is timely. This paper highlights past and present methods, categorizing them into segmentation, atlas/mapping, and emission-based schemes. Potential strategies for future developments are also presented., (© 2019 The Authors. Medical Physics published by Wiley Periodicals, Inc. on behalf of American Association of Physicists in Medicine.)

Objective: To find the covert patterns of abnormality in patients with unilateral temporal lobe epilepsy (TLE) and visually normal brain magnetic resonance images (MRI-negative), comparing them to those with visible abnormalities (MRI-positive)., Methods: We used multimodal brain MRI from patients with unilateral TLE and employed contemporary machine learning methods to predict the known laterality of seizure onset in 104 subjects (82 MRI-positive, 22 MRI-negative). A visualization approach entitled "Importance Maps" was developed to highlight image features predictive of seizure laterality in both the MRI-positive and MRI-negative cases., Results: Seizure laterality could be predicted with an area under the receiver operating characteristic curve of 0.981 (95% confidence interval [CI] =0.974-0.989) in MRI-positive and 0.842 (95% CI = 0.736-0.949) in MRI-negative cases. The known image features arising from the hippocampus were the leading predictors of seizure laterality in the MRI-positive cases, whereas widespread temporal lobe abnormalities were revealed in the MRI-negative cases., Significance: Covert abnormalities not discerned on visual reading were detected in MRI-negative TLE, with a spatial pattern involving the whole temporal lobe, rather than just the hippocampus. This suggests that MRI-negative TLE may be associated with subtle but widespread temporal lobe abnormalities. These abnormalities merit close inspection and postacquisition processing if there is no overt lesion., (© 2019 The Authors. Epilepsia published by Wiley Periodicals, Inc. on behalf of International League Against Epilepsy.)

Quantification of cerebral white matter hyperintensities (WMH) of presumed vascular origin is of key importance in many neurological research studies. Currently, measurements are often still obtained from manual segmentations on brain MR images, which is a laborious procedure. The automatic WMH segmentation methods exist, but a standardized comparison of the performance of such methods is lacking. We organized a scientific challenge, in which developers could evaluate their methods on a standardized multi-center/-scanner image dataset, giving an objective comparison: the WMH Segmentation Challenge. Sixty T1 + FLAIR images from three MR scanners were released with the manual WMH segmentations for training. A test set of 110 images from five MR scanners was used for evaluation. The segmentation methods had to be containerized and submitted to the challenge organizers. Five evaluation metrics were used to rank the methods: 1) Dice similarity coefficient; 2) modified Hausdorff distance (95th percentile); 3) absolute log-transformed volume difference; 4) sensitivity for detecting individual lesions; and 5) F1-score for individual lesions. In addition, the methods were ranked on their inter-scanner robustness; 20 participants submitted their methods for evaluation. This paper provides a detailed analysis of the results. In brief, there is a cluster of four methods that rank significantly better than the other methods, with one clear winner. The inter-scanner robustness ranking shows that not all the methods generalize to unseen scanners. The challenge remains open for future submissions and provides a public platform for method evaluation.

There are numerous challenges to identifying, developing and implementing quantitative techniques for use in clinical radiology, suggesting the need for a common translational pathway. We developed the quantitative neuroradiology initiative (QNI), as a model framework for the technical and clinical validation necessary to embed automated segmentation and other image quantification software into the clinical neuroradiology workflow. We hypothesize that quantification will support reporters with clinically relevant measures contextualized with normative data, increase the precision of longitudinal comparisons, and generate more consistent reporting across levels of radiologists' experience. The QNI framework comprises the following steps: (1) establishing an area of clinical need and identifying the appropriate proven imaging biomarker(s) for the disease in question; (2) developing a method for automated analysis of these biomarkers, by designing an algorithm and compiling reference data; (3) communicating the results via an intuitive and accessible quantitative report; (4) technically and clinically validating the proposed tool pre-use; (5) integrating the developed analysis pipeline into the clinical reporting workflow; and (6) performing in-use evaluation. We will use current radiology practice in dementia as an example, where radiologists have established visual rating scales to describe the degree and pattern of atrophy they detect. These can be helpful, but are somewhat subjective and coarse classifiers, suffering from floor and ceiling limitations. Meanwhile, several imaging biomarkers relevant to dementia diagnosis and management have been proposed in the literature; some clinically approved radiology software tools exist but in general, these have not undergone rigorous clinical validation in high volume or in tertiary dementia centres. The QNI framework aims to address this need. Quantitative image analysis is developing apace within the research domain. Translating quantitative techniques into the clinical setting presents significant challenges, which must be addressed to meet the increasing demand for accurate, timely and impactful clinical imaging services.

We assessed whether HIV status was associated with white matter hyperintensities (WMH), a neuroimaging correlate of cerebral small vessel disease (CSVD), in men aged ≥50 years. A cross-sectional substudy was nested within a larger cohort study. Virologically suppressed men living with HIV (MLWH) and demographically matched HIV-negative men aged ≥50 underwent magnetic resonance imaging (MRI) at 3 Tesla. Sequences included volumetric three-dimensional (3D) T1-weighted, fluid-attenuated inversion recovery and pseudocontinuous arterial spin labeling. Regional segmentation by automated image processing algorithms was used to extract WMH volume (WMHV) and resting cerebral blood flow (CBF). The association between HIV status and WMHV as a proportion of intracranial volume (ICV; log-transformed) was estimated using a multivariable linear regression model. Thirty-eight MLWH [median age 59 years (interquartile range, IQR 55-64)] and 37 HIV-negative [median 58 years (54-63)] men were analyzed. MLWH had median CD4 + count 570 (470-700) cells/μL and a median time since diagnosis of 20 (14-24) years. Framingham 10-year risk of cardiovascular disease was 6.5% in MLWH and 7.4% in controls. Two (5%) MLWH reported a history of stroke or transient ischemic attack and five (13%) reported coronary heart disease compared with none of the controls. The total WMHV in MLWH was 1,696 μL (IQR 1,229-3,268 μL) or 0.10% of ICV compared with 1,627 μL (IQR 1,032-3,077 μL), also 0.10% of ICV in the HIV-negative group ( p  = .43). In the multivariable model, WMHV/ICV was not associated with HIV status ( p  = .86). There was an age-dependent decline in cortical CBF [-3.9 mL/100 mL/min per decade of life (95% confidence interval 1.1-6.7 mL)] but no association between CBF and HIV status ( p  > .2 in all brain regions analyzed). In conclusion, we found no quantitative MRI evidence of an increased burden of CSVD in MLWH aged 50 years and older.

Disease progression modeling (DPM) using longitudinal data is a challenging machine learning task. Existing DPM algorithms neglect temporal dependencies among measurements, make parametric assumptions about biomarker trajectories, do not model multiple biomarkers jointly, and need an alignment of subjects' trajectories. In this paper, recurrent neural networks (RNNs) are utilized to address these issues. However, in many cases, longitudinal cohorts contain incomplete data, which hinders the application of standard RNNs and requires a pre-processing step such as imputation of the missing values. Instead, we propose a generalized training rule for the most widely used RNN architecture, long short-term memory (LSTM) networks, that can handle both missing predictor and target values. The proposed LSTM algorithm is applied to model the progression of Alzheimer's disease (AD) using six volumetric magnetic resonance imaging (MRI) biomarkers, i.e., volumes of ventricles, hippocampus, whole brain, fusiform, middle temporal gyrus, and entorhinal cortex, and it is compared to standard LSTM networks with data imputation and a parametric, regression-based DPM method. The results show that the proposed algorithm achieves a significantly lower mean absolute error (MAE) than the alternatives with p < 0.05 using Wilcoxon signed rank test in predicting values of almost all of the MRI biomarkers. Moreover, a linear discriminant analysis (LDA) classifier applied to the predicted biomarker values produces a significantly larger area under the receiver operating characteristic curve (AUC) of 0.90 vs. at most 0.84 with p < 0.001 using McNemar's test for clinical diagnosis of AD. Inspection of MAE curves as a function of the amount of missing data reveals that the proposed LSTM algorithm achieves the best performance up until more than 74% missing values. Finally, it is illustrated how the method can successfully be applied to data with varying time intervals. This paper shows that built-in handling of missing values in training an LSTM network benefits the application of RNNs in neurodegenerative disease progression modeling in longitudinal cohorts., (Copyright © 2019 Elsevier B.V. All rights reserved.)

Multi-Atlas based Segmentation (MAS) algorithms have been successfully applied to many medical image segmentation tasks, but their success relies on a large number of atlases and good image registration performance. Choosing well-registered atlases for label fusion is vital for an accurate segmentation. This choice becomes even more crucial when the segmentation involves organs characterized by a high anatomical and pathological variability. In this paper, we propose a new genetic atlas selection strategy (GAS) that automatically chooses the best subset of atlases to be used for segmenting the target image, on the basis of both image similarity and segmentation overlap. More precisely, the key idea of GAS is that if two images are similar, the performances of an atlas for segmenting each image are similar. Since the ground truth of each atlas is known, GAS first selects a predefined number of similar images to the target, then, for each one of them, finds a near-optimal subset of atlases by means of a genetic algorithm. All these near-optimal subsets are then combined and used to segment the target image. GAS was tested on single-label and multi-label segmentation problems. In the first case, we considered the segmentation of both the whole prostate and of the left ventricle of the heart from magnetic resonance images. Regarding multi-label problems, the zonal segmentation of the prostate into peripheral and transition zone was considered. The results showed that the performance of MAS algorithms statistically improved when GAS is used., (Copyright © 2018 Elsevier B.V. All rights reserved.)

Background Cardiovascular risk factors vary between ethnicities but little is known about their differential effects on white matter hyperintensities ( WMH ), an indicator of brain aging and burden of cerebrovascular disease. Methods and Results Brain magnetic resonance imaging scans from 213 people of South Asian and 256 of European ethnicity (total=469) were analyzed for global and regional WMH load. Associations with cardiovascular risk factors and a composite cardiovascular risk score (National Cholesterol Education Programme Adult Treatment Panel III) were compared by ethnicity, diabetes mellitus, smoking, and hypertension status. Distributional patterns of WMH were similar by ethnicity but the vulnerability to specific risk factors differed. Associations between WMH and age or National Cholesterol Education Programme Adult Treatment Panel III scores were stronger in South Asians compared with Europeans. For instance, a year of age led to an excess of 3.8% (confidence interval=[0.2, 7.6]; P=0.04) of WMH load in frontal regions in South Asians compared with Europeans. In the diabetic subgroup, South Asians had more WMH than Europeans (+63.3%, confidence interval=[14.1, 133.9]; P=0.007), particularly in the deeper regions (+102% confidence interval=[24, 329]; P=0.004). In the population as a whole, diabetes mellitus was not, or only weakly, related to an increase in WMH volume (12.4%, confidence interval=[-10.7, 41.3]; P=0.32), and diabetes mellitus duration was a positive predictor of frontal periventricular WMH load in Europeans but not in South Asians. In turn, diastolic blood pressure was positively associated with WMH volumes in South Asians but not in Europeans. Hypertension was not associated with WMH load ( P=0.9). Conclusions Distribution patterns of WMH are similar in South Asians and Europeans but older age and higher cardiovascular risk are associated with more WMH in South Asians.

Owing to its excellent soft-tissue contrast, magnetic resonance (MR) imaging has found an increased application in radiation therapy (RT). By harnessing these properties for treatment planning, automated segmentation methods can alleviate the manual workload burden to the clinical workflow. We investigated atlas-based segmentation methods of organs at risk (OARs) in the head and neck (H&N) region using one approach that selected the most similar atlas from a library of segmented images and two multi-atlas approaches. The latter were based on weighted majority voting and an iterative atlas-fusion approach called STEPS. We built the atlas library from pre-treatment T1-weighted MR images of 12 patients with manual contours of the parotids, spinal cord and mandible, delineated by a clinician. Following a leave-one-out cross-validation strategy, we measured the geometric accuracy by calculating Dice similarity coefficients (DSC), standard and 95% Hausdorff distances (HD and HD95), and the mean surface distance (MSD), whereby the manual contours served as the gold standard. To benchmark the algorithm, we determined the inter-observer variability (IOV) between three observers. To investigate the dosimetric effect of segmentation inaccuracies, we implemented an auto-planning strategy within the treatment planning system Monaco (Elekta AB, Stockholm, Sweden). For each set of auto-segmented OARs, we generated a plan for a 9-beam step and shoot intensity modulated RT treatment, designed according to our institution's clinical H&N protocol. Superimposing the dose distributions on the gold standard OARs, we calculated dose differences to OARs caused by delineation differences between auto-segmented and gold standard OARs. We investigated the correlations between geometric and dosimetric differences. The mean DSC was larger than 0.8 and the mean MSD smaller than 2 mm for the multi-atlas approaches, resulting in a geometric accuracy comparable to previously published results and within the range of the IOV. While dosimetric differences could be as large as 23% of the clinical goal, treatment plans fulfilled all imposed clinical goals for the gold standard OARs. Correlations between geometric and dosimetric measures were low with R 2   <  0.5. The geometric accuracy and the ability to achieve clinically acceptable treatment plans indicate the suitability of using atlas-based contours for RT treatment planning purposes. The low correlations between geometric and dosimetric measures suggest that geometric measures alone are not sufficient to predict the dosimetric impact of segmentation inaccuracies on treatment planning for the data utilised in this study.

Background and Purpose: Visual rating scales have limited capacities to depict the regional distribution of cerebral white matter hyperintensities (WMH). We present a regional-zonal volumetric analysis alongside a visualization tool to compare and deconstruct visual rating scales., Materials and Methods: 3D T1-weighted, T2-weighted spin-echo and FLAIR images were acquired on a 3T system, from 82 elderly participants in a population-based study. Images were automatically segmented for WMH. Lobar boundaries and distance to ventricular surface were used to define white matter regions. Regional-zonal WMH loads were displayed using bullseye plots. Four raters assessed all images applying three scales. Correlations between visual scales and regional WMH as well as inter and intra-rater variability were assessed. A multinomial ordinal regression model was used to predict scores based on regional volumes and global WMH burdens., Results: On average, the bullseye plot depicted a right-left symmetry in the distribution and concentration of damage in the periventricular zone, especially in frontal regions. WMH loads correlated well with the average visual rating scores (e.g. Kendall's tau [Volume, Scheltens]=0.59 CI=[0.53 0.62]). Local correlations allowed comparison of loading patterns between scales and between raters. Regional measurements had more predictive power than global WMH burden (e.g. frontal caps prediction with local features: ICC=0.67 CI=[0.53 0.77], global volume=0.50 CI=[0.32 0.65], intra-rater=0.44 CI=[0.23 0.60])., Conclusion: Regional-zonal representation of WMH burden highlights similarities and differences between visual rating scales and raters. The bullseye infographic tool provides a simple visual representation of regional lesion load that can be used for rater calibration and training., (Copyright © 2017 The Authors. Published by Elsevier Masson SAS.. All rights reserved.)

In this paper we present a hierarchical generative model of medical image data, which can capture simultaneously the variability of both signal intensity and anatomical shapes across large populations. Such a model has a direct application for learning average-shaped probabilistic tissue templates in a fully automated manner. While in principle the generality of the proposed Bayesian approach makes it suitable to address a wide range of medical image computing problems, our work focuses primarily on neuroimaging applications. In particular we validate the proposed method on both real and synthetic brain MR scans including the cervical cord and demonstrate that it yields accurate alignment of brain and spinal cord structures, as compared to state-of-the-art tools for medical image registration. At the same time we illustrate how the resulting tissue probability maps can readily be used to segment, bias correct and spatially normalise unseen data, which are all crucial pre-processing steps for MR imaging studies., (Copyright © 2017 The Authors. Published by Elsevier Inc. All rights reserved.)

The aim of this study was to assess whether the use of accelerated MRI scans in place of non-accelerated scans influenced brain volume and atrophy rate measures in controls and subjects with mild cognitive impairment and Alzheimer's disease. We used data from 861 subjects at baseline, 573 subjects at 6 months and 384 subjects at 12 months from the Alzheimer's Disease Neuroimaging Initiative (ADNI). We calculated whole-brain, ventricular and hippocampal atrophy rates using the k-means boundary shift integral (BSI). Scan quality was visually assessed and the proportion of good quality accelerated and non-accelerated scans compared. We also compared MMSE scores, vascular burden and age between subjects with poor quality scans with those with good quality scans. Finally, we estimated sample size requirements for a hypothetical clinical trial when using atrophy rates from accelerated scans and non-accelerated scans. No significant differences in whole-brain, ventricular and hippocampal volumes and atrophy rates were found between accelerated and non-accelerated scans. Twice as many non-accelerated scan pairs suffered from at least some motion artefacts compared with accelerated scan pairs (p ≤ 0.001), which may influence the BSI. Subjects whose accelerated scans had significant motion had a higher mean vascular burden and age (p ≤ 0.05) whilst subjects whose non-accelerated scans had significant motion had poorer MMSE scores (p ≤ 0.05). No difference in estimated sample size requirements was found when using accelerated vs. non-accelerated scans. Accelerated scans reduce scan time and are better tolerated. Therefore it may be advantageous to use accelerated over non-accelerated scans in clinical trials that use ADNI-type protocols, especially in more cognitively impaired subjects.

Background and Purpose: Computed tomography (CT) imaging is the current gold standard for radiotherapy treatment planning (RTP). The establishment of a magnetic resonance imaging (MRI) only RTP workflow requires the generation of a synthetic CT (sCT) for dose calculation. This study evaluates the feasibility of using a multi-atlas sCT synthesis approach (sCT a ) for head and neck and prostate patients., Material and Methods: The multi-atlas method was based on pairs of non-rigidly aligned MR and CT images. The sCT a was obtained by registering the MRI atlases to the patient's MRI and by fusing the mapped atlases according to morphological similarity to the patient. For comparison, a bulk density assignment approach (sCT bda ) was also evaluated. The sCT bda was obtained by assigning density values to MRI tissue classes (air, bone and soft-tissue). After evaluating the synthesis accuracy of the sCTs (mean absolute error), sCT-based delineations were geometrically compared to the CT-based delineations. Clinical plans were re-calculated on both sCTs and a dose-volume histogram and a gamma analysis was performed using the CT dose as ground truth., Results: Results showed that both sCTs were suitable to perform clinical dose calculations with mean dose differences less than 1% for both the planning target volume and the organs at risk. However, only the sCT a provided an accurate and automatic delineation of bone., Conclusions: Combining MR delineations with our multi-atlas CT synthesis method could enable MRI-only treatment planning and thus improve the dosimetric and geometric accuracy of the treatment, and reduce the number of imaging procedures., (Copyright © 2017 Associazione Italiana di Fisica Medica. All rights reserved.)

Aim: To accurately quantify the radioactivity concentration measured by PET, emission data need to be corrected for photon attenuation; however, the MRI signal cannot easily be converted into attenuation values, making attenuation correction (AC) in PET/MRI challenging. In order to further improve the current vendor-implemented MR-AC methods for absolute quantification, a number of prototype methods have been proposed in the literature. These can be categorized into three types: template/atlas-based, segmentation-based, and reconstruction-based. These proposed methods in general demonstrated improvements compared to vendor-implemented AC, and many studies report deviations in PET uptake after AC of only a few percent from a gold standard CT-AC. Using a unified quantitative evaluation with identical metrics, subject cohort, and common CT-based reference, the aims of this study were to evaluate a selection of novel methods proposed in the literature, and identify the ones suitable for clinical use., Methods: In total, 11 AC methods were evaluated: two vendor-implemented (MR-AC DIXON and MR-AC UTE ), five based on template/atlas information (MR-AC SEGBONE (Koesters et al., 2016), MR-AC ONTARIO (Anazodo et al., 2014), MR-AC BOSTON (Izquierdo-Garcia et al., 2014), MR-AC UCL (Burgos et al., 2014), and MR-AC MAXPROB (Merida et al., 2015)), one based on simultaneous reconstruction of attenuation and emission (MR-AC MLAA (Benoit et al., 2015)), and three based on image-segmentation (MR-AC MUNICH (Cabello et al., 2015), MR-AC CAR-RiDR (Juttukonda et al., 2015), and MR-AC RESOLUTE (Ladefoged et al., 2015)). We selected 359 subjects who were scanned using one of the following radiotracers: [ 18 F]FDG (210), [ 11 C]PiB (51), and [ 18 F]florbetapir (98). The comparison to AC with a gold standard CT was performed both globally and regionally, with a special focus on robustness and outlier analysis., Results: The average performance in PET tracer uptake was within ±5% of CT for all of the proposed methods, with the average±SD global percentage bias in PET FDG uptake for each method being: MR-AC DIXON (-11.3±3.5)%, MR-AC UTE (-5.7±2.0)%, MR-AC ONTARIO (-4.3±3.6)%, MR-AC MUNICH (3.7±2.1)%, MR-AC MLAA (-1.9±2.6)%, MR-AC SEGBONE (-1.7±3.6)%, MR-AC UCL (0.8±1.2)%, MR-AC CAR-RiDR (-0.4±1.9)%, MR-AC MAXPROB (-0.4±1.6)%, MR-AC BOSTON (-0.3±1.8)%, and MR-AC RESOLUTE (0.3±1.7)%, ordered by average bias. The overall best performing methods (MR-AC BOSTON , MR-AC MAXPROB , MR-AC RESOLUTE and MR-AC UCL , ordered alphabetically) showed regional average errors within ±3% of PET with CT-AC in all regions of the brain with FDG, and the same four methods, as well as MR-AC CAR-RiDR , showed that for 95% of the patients, 95% of brain voxels had an uptake that deviated by less than 15% from the reference. Comparable performance was obtained with PiB and florbetapir., Conclusions: All of the proposed novel methods have an average global performance within likely acceptable limits (±5% of CT-based reference), and the main difference among the methods was found in the robustness, outlier analysis, and clinical feasibility. Overall, the best performing methods were MR-ACBOSTON, MR-ACMAXPROB, MR-ACRESOLUTE and MR-ACUCL, ordered alphabetically. These methods all minimized the number of outliers, standard deviation, and average global and local error. The methods MR-ACMUNICH and MR-ACCAR-RiDR were both within acceptable quantitative limits, so these methods should be considered if processing time is a factor. The method MR-ACSEGBONE also demonstrates promising results, and performs well within the likely acceptable quantitative limits. For clinical routine scans where processing time can be a key factor, this vendor-provided solution currently outperforms most methods. With the performance of the methods presented here, it may be concluded that the challenge of improving the accuracy of MR-AC in adult brains with normal anatomy has been solved to a quantitatively acceptable degree, which is smaller than the quantification reproducibility in PET imaging., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Objective: To assess the association between proximity to the inner (ventricular and aqueductal) and outer (pial) surfaces of the brain and the distribution of normal appearing white matter (NAWM) and grey matter (GM) abnormalities, and white matter (WM) lesions, in multiple sclerosis (MS)., Methods: 67 people with relapse-onset MS and 30 healthy controls were included in the study. Volumetric T1 images and high-resolution (1 mm 3 ) magnetisation transfer ratio (MTR) images were acquired and segmented into 12 bands between the inner and outer surfaces of the brain. The first and last bands were discarded to limit partial volume effects with cerebrospinal fluid. MTR values were computed for all bands in supratentorial NAWM, cerebellar NAWM and brainstem NA tissue, and deep and cortical GM. Band WM lesion volumes were also measured., Results: Proximity to the ventricular surfaces was associated with progressively lower MTR values in the MS group but not in controls in supratentorial and cerebellar NAWM, brainstem NA and in deep and cortical GM. The density of WM lesions was associated with proximity to the ventricles only in the supratentorial compartment, and no link was found with distance from the pial surfaces., Conclusions: In MS, MTR abnormalities in NAWM and GM are related to distance from the inner and outer surfaces of the brain, and this suggests that there is a common factor underlying their spatial distribution. A similar pattern was not found for WM lesions, raising the possibility that different factors promote their formation., (Published by the BMJ Publishing Group Limited. For permission to use (where not already granted under a licence) please go to http://www.bmj.com/company/products-services/rights-and-licensing/.)

Axonal loss in the spinal cord is one of the main contributing factors to irreversible clinical disability in multiple sclerosis (MS). In vivo axonal loss can be assessed indirectly by estimating a reduction in the cervical cross-sectional area (CSA) of the spinal cord over time, which is indicative of spinal cord atrophy, and such a measure may be obtained by means of image segmentation using magnetic resonance imaging (MRI). In this work, we propose a new fully automated spinal cord segmentation technique that incorporates two different multi-atlas segmentation propagation and fusion techniques: The Optimized PatchMatch Label fusion (OPAL) algorithm for localising and approximately segmenting the spinal cord, and the Similarity and Truth Estimation for Propagated Segmentations (STEPS) algorithm for segmenting white and grey matter simultaneously. In a retrospective analysis of MRI data, the proposed method facilitated CSA measurements with accuracy equivalent to the inter-rater variability, with a Dice score (DSC) of 0.967 at C2/C3 level. The segmentation performance for grey matter at C2/C3 level was close to inter-rater variability, reaching an accuracy (DSC) of 0.826 for healthy subjects and 0.835 people with clinically isolated syndrome MS., Competing Interests: FP, MJC, MY, LH, ET, HK, ML, CWK and SO declare no competing financial interests. DM has received honoraria through payments to his employer, UCL Institute of Neurology, for Advisory Committee and/or Consultancy advice in multiple sclerosis studies from Biogen Idec, GlaxoSmithKline, Novartis, Merck, Chugai, Mitsubishi Pharma Europe and Bayer Schering Pharma. He has also received compensation through payments to my employer for performing central MRI analysis of multiple sclerosis trials from GlaxoSmithKline, Biogen Idec, Novartis and Merck. OC is a consultant for Novartis, Biogen-Idec, Genzyme and General Electric, and all the payments are made to the UCL Institute of Neurology. She is an associate editor of Neurology.

Unlabelled: In this work, we assessed the feasibility of attenuation correction (AC) based on a multi-atlas-based method (m-Atlas) by comparing it with a clinical AC method (single-atlas-based method [s-Atlas]), on a time-of-flight (TOF) PET/MRI scanner., Methods: We enrolled 15 patients. The median patient age was 59 y (age range, 31-80). All patients underwent clinically indicated whole-body (18)F-FDG PET/CT for staging, restaging, or follow-up of malignant disease. All patients volunteered for an additional PET/MRI scan of the head (no additional tracer being injected). For each patient, 3 AC maps were generated. Both s-Atlas and m-Atlas AC maps were generated from the same patient-specific LAVA-Flex T1-weighted images being acquired by default on the PET/MRI scanner during the first 18 s of the PET scan. An s-Atlas AC map was extracted by the PET/MRI scanner, and an m-Atlas AC map was created using a Web service tool that automatically generates m-Atlas pseudo-CT images. For comparison, the AC map generated by PET/CT was registered and used as a gold standard. PET images were reconstructed from raw data on the TOF PET/MRI scanner using each AC map. All PET images were normalized to the SPM5 PET template, and (18)F-FDG accumulation was quantified in 67 volumes of interest (VOIs; automated anatomic labeling atlas). Relative (%diff) and absolute differences (|%diff|) between images based on each atlas AC and CT-AC were calculated. (18)F-FDG uptake in all VOIs and generalized merged VOIs were compared using the paired t test and Bland-Altman test., Results: The range of error on m-Atlas in all 1,005 VOIs was -4.99% to 4.09%. The |%diff| on the m-Atlas was improved by about 20% compared with s-Atlas (s-Atlas vs. m-Atlas: 1.49% ± 1.06% vs. 1.21% ± 0.89%, P < 0.01). In generalized VOIs, %diff on m-Atlas in the temporal lobe and cerebellum was significantly smaller (s-Atlas vs. m-Atlas: temporal lobe, 1.49% ± 1.37% vs. -0.37% ± 1.41%, P < 0.01; cerebellum, 1.55% ± 1.97% vs. -1.15% ± 1.72%, P < 0.01)., Conclusion: The errors introduced using either s-Atlas or m-Atlas did not exceed 5% in any brain region investigated. When compared with the clinical s-Atlas, m-Atlas is more accurate, especially in regions close to the skull base., (© 2016 by the Society of Nuclear Medicine and Molecular Imaging, Inc.)

Introduction: Infants born extremely preterm (<28 weeks of gestation) are at risk of significant neurodevelopmental sequelae. In these infants birth coincides with a period of rapid brain growth and development, when the brain is also vulnerable to a range of insults. Mapping these changes is crucial for identifying potential biomarkers to predict early impairment., Methods: In this study we use surface-based spectral matching techniques to find an intrasubject longitudinal surface correspondence between the white-grey matter boundary at 30 and 40 weeks equivalent gestational age in nine extremely preterm born infants., Results: Using the resulting surface correspondence, we identified regions that undergo more cortical folding of the white-grey matter boundary during the preterm period by looking at changes in well-known curvature measures. We performed Hotelling T(2) statistics to evaluate the significance of our findings., Discussion: The prefrontal and temporal lobes exhibit most development during the preterm period, especially in the left hemisphere. Such correspondences are a promising result as longitudinal measurements of change in cortical folding could provide insightful information about the mechanical properties of the underlying tissue and may be useful in inferring changes during growth and development in this vulnerable period.

Mouse models of Alzheimer's disease have served as valuable tools for investigating pathogenic mechanisms relating to neurodegeneration, including tau-mediated and neurofibrillary tangle pathology-a major hallmark of the disease. In this work, we have used multiparametric magnetic resonance imaging (MRI) in a longitudinal study of neurodegeneration in the rTg4510 mouse model of tauopathy, a subset of which were treated with doxycycline at different time points to suppress the tau transgene. Using this paradigm, we investigated the sensitivity of multiparametric MRI to both the accumulation and suppression of pathologic tau. Tau-related atrophy was discernible from 5.5 months within the cortex and hippocampus. We observed markedly less atrophy in the treated rTg4510 mice, which was enhanced after doxycycline intervention from 3.5 months. We also observed differences in amide proton transfer, cerebral blood flow, and diffusion tensor imaging parameters in the rTg4510 mice, which were significantly less altered after doxycycline treatment. We propose that these non-invasive MRI techniques offer insight into pathologic mechanisms underpinning Alzheimer's disease that may be important when evaluating emerging therapeutics targeting one of more of these processes., (Copyright © 2016 The Authors. Published by Elsevier Inc. All rights reserved.)

Structural MRI is widely used for investigating brain atrophy in many neurodegenerative disorders, with several research groups developing and publishing techniques to provide quantitative assessments of this longitudinal change. Often techniques are compared through computation of required sample size estimates for future clinical trials. However interpretation of such comparisons is rendered complex because, despite using the same publicly available cohorts, the various techniques have been assessed with different data exclusions and different statistical analysis models. We created the MIRIAD atrophy challenge in order to test various capabilities of atrophy measurement techniques. The data consisted of 69 subjects (46 Alzheimer's disease, 23 control) who were scanned multiple (up to twelve) times at nine visits over a follow-up period of one to two years, resulting in 708 total image sets. Nine participating groups from 6 countries completed the challenge by providing volumetric measurements of key structures (whole brain, lateral ventricle, left and right hippocampi) for each dataset and atrophy measurements of these structures for each time point pair (both forward and backward) of a given subject. From these results, we formally compared techniques using exactly the same dataset. First, we assessed the repeatability of each technique using rates obtained from short intervals where no measurable atrophy is expected. For those measures that provided direct measures of atrophy between pairs of images, we also assessed symmetry and transitivity. Then, we performed a statistical analysis in a consistent manner using linear mixed effect models. The models, one for repeated measures of volume made at multiple time-points and a second for repeated "direct" measures of change in brain volume, appropriately allowed for the correlation between measures made on the same subject and were shown to fit the data well. From these models, we obtained estimates of the distribution of atrophy rates in the Alzheimer's disease (AD) and control groups and of required sample sizes to detect a 25% treatment effect, in relation to healthy ageing, with 95% significance and 80% power over follow-up periods of 6, 12, and 24months. Uncertainty in these estimates, and head-to-head comparisons between techniques, were carried out using the bootstrap. The lateral ventricles provided the most stable measurements, followed by the brain. The hippocampi had much more variability across participants, likely because of differences in segmentation protocol and less distinct boundaries. Most methods showed no indication of bias based on the short-term interval results, and direct measures provided good consistency in terms of symmetry and transitivity. The resulting annualized rates of change derived from the model ranged from, for whole brain: -1.4% to -2.2% (AD) and -0.35% to -0.67% (control), for ventricles: 4.6% to 10.2% (AD) and 1.2% to 3.4% (control), and for hippocampi: -1.5% to -7.0% (AD) and -0.4% to -1.4% (control). There were large and statistically significant differences in the sample size requirements between many of the techniques. The lowest sample sizes for each of these structures, for a trial with a 12month follow-up period, were 242 (95% CI: 154 to 422) for whole brain, 168 (95% CI: 112 to 282) for ventricles, 190 (95% CI: 146 to 268) for left hippocampi, and 158 (95% CI: 116 to 228) for right hippocampi. This analysis represents one of the most extensive statistical comparisons of a large number of different atrophy measurement techniques from around the globe. The challenge data will remain online and publicly available so that other groups can assess their methods., (Copyright © 2015. Published by Elsevier Inc.)

This paper introduces a novel method for inferring spatially varying regularisation in non-linear registration. This is achieved through full Bayesian inference on a probabilistic registration model, where the prior on the transformation parameters is parameterised as a weighted mixture of spatially localised components. Such an approach has the advantage of allowing the registration to be more flexibly driven by the data than a traditional globally defined regularisation penalty, such as bending energy. The proposed method adaptively determines the influence of the prior in a local region. The strength of the prior may be reduced in areas where the data better support deformations, or can enforce a stronger constraint in less informative areas. Consequently, the use of such a spatially adaptive prior may reduce unwanted impacts of regularisation on the inferred transformation. This is especially important for applications where the deformation field itself is of interest, such as tensor based morphometry. The proposed approach is demonstrated using synthetic images, and with application to tensor based morphometry analysis of subjects with Alzheimer's disease and healthy controls. The results indicate that using the proposed spatially adaptive prior leads to sparser deformations, which provide better localisation of regional volume change. Additionally, the proposed regularisation model leads to more data driven and localised maps of registration uncertainty. This paper also demonstrates for the first time the use of Bayesian model comparison for selecting different types of regularisation., (Copyright © 2015 The Authors. Published by Elsevier B.V. All rights reserved.)

Positron Emission Tomography/Magnetic Resonance Imaging (PET/MR) scanners are expected to offer a new range of clinical applications. Attenuation correction is an essential requirement for quantification of PET data but MRI images do not directly provide a patient-specific attenuation map. Methods We further validate and extend a Computed Tomography (CT) and attenuation map (μ-map) synthesis method based on pre-acquired MRI-CT image pairs. The validation consists of comparing the CT images synthesised with the proposed method to the original CT images. PET images were acquired using two different tracers ((18)F-FDG and (18)F-florbetapir). They were then reconstructed and corrected for attenuation using the synthetic μ-maps and compared to the reference PET images corrected with the CT-based μ-maps. During the validation, we observed that the CT synthesis was inaccurate in areas such as the neck and the cerebellum, and propose a refinement to mitigate these problems, as well as an extension of the method to multi-contrast MRI data. Results With the improvements proposed, a significant enhancement in CT synthesis, which results in a reduced absolute error and a decrease in the bias when reconstructing PET images, was observed. For both tracers, on average, the absolute difference between the reference PET images and the PET images corrected with the proposed method was less than 2%, with a bias inferior to 1%. Conclusion With the proposed method, attenuation information can be accurately derived from MRI images by synthesising CT using routine anatomical sequences. MRI sequences, or combination of sequences, can be used to synthesise CT images, as long as they provide sufficient anatomical information.

Preterm birth is a major public health concern, with the severity and occurrence of adverse outcome increasing with earlier delivery. Being born preterm disrupts a time of rapid brain development: in addition to volumetric growth, the cortex folds, myelination is occurring and there are changes on the cellular level. These neurological events have been imaged non-invasively using diffusion-weighted (DW) MRI. In this population, there has been a focus on examining diffusion in the white matter, but the grey matter is also critically important for neurological health. We acquired multi-shell high-resolution diffusion data on 12 infants born at ≤ 28 weeks of gestational age at two time-points: once when stable after birth, and again at term-equivalent age. We used the Neurite Orientation Dispersion and Density Imaging model (NODDI) (Zhang et al., 2012) to analyse the changes in the cerebral cortex and the thalamus, both grey matter regions. We showed region-dependent changes in NODDI parameters over the preterm period, highlighting underlying changes specific to the microstructure. This work is the first time that NODDI parameters have been evaluated in both the cortical and the thalamic grey matter as a function of age in preterm infants, offering a unique insight into neuro-development in this at-risk population., (Copyright © 2015 The Authors. Published by Elsevier Inc. All rights reserved.)

As the number of people diagnosed with Alzheimer's disease (AD) reaches epidemic proportions, there is an urgent need to develop effective treatment strategies to tackle the social and economic costs of this fatal condition. Dozens of candidate therapeutics are currently being tested in clinical trials, and compounds targeting the aberrant accumulation of tau proteins into neurofibrillary tangles (NFTs) are the focus of substantial current interest. Reliable, translatable biomarkers sensitive to both tau pathology and its modulation by treatment along with animal models that faithfully reflect aspects of the human disease are urgently required. Magnetic resonance imaging (MRI) is well established as a valuable tool for monitoring the structural brain changes that accompany AD progression. However the descent into dementia is not defined by macroscopic brain matter loss alone: non-invasive imaging measurements sensitive to protein accumulation, white matter integrity and cerebral haemodynamics probe distinct aspects of AD pathophysiology and may serve as superior biomarkers for assessing drug efficacy. Here we employ a multi-parametric array of five translatable MRI techniques to characterise the in vivo pathophysiological phenotype of the rTg4510 mouse model of tauopathy (structural imaging, diffusion tensor imaging (DTI), arterial spin labelling (ASL), chemical exchange saturation transfer (CEST) and glucose CEST). Tau-induced pathological changes included grey matter atrophy, increased radial diffusivity in the white matter, decreased amide proton transfer and hyperperfusion. We demonstrate that the above markers unambiguously discriminate between the transgenic group and age-matched controls and provide a comprehensive profile of the multifaceted neuropathological processes underlying the rTg4510 model. Furthermore, we show that ASL and DTI techniques offer heightened sensitivity to processes believed to precede detectable structural changes and, as such, provides a platform for the study of disease mechanisms and therapeutic intervention., (Copyright © 2015 Elsevier Inc. All rights reserved.)

A number of algorithms for brain segmentation in preterm born infants have been published, but a reliable comparison of their performance is lacking. The NeoBrainS12 study (http://neobrains12.isi.uu.nl), providing three different image sets of preterm born infants, was set up to provide such a comparison. These sets are (i) axial scans acquired at 40 weeks corrected age, (ii) coronal scans acquired at 30 weeks corrected age and (iii) coronal scans acquired at 40 weeks corrected age. Each of these three sets consists of three T1- and T2-weighted MR images of the brain acquired with a 3T MRI scanner. The task was to segment cortical grey matter, non-myelinated and myelinated white matter, brainstem, basal ganglia and thalami, cerebellum, and cerebrospinal fluid in the ventricles and in the extracerebral space separately. Any team could upload the results and all segmentations were evaluated in the same way. This paper presents the results of eight participating teams. The results demonstrate that the participating methods were able to segment all tissue classes well, except myelinated white matter., (Copyright © 2014 Elsevier B.V. All rights reserved.)

The advent of large of multi-modal imaging databases opens up the opportunity to learn how local intensity patterns covariate between multiple modalities. These models can then be used to describe expected intensities in an unseen image modalities given one or multiple observations, or to detect deviations (e.g. pathology) from the expected intensity patterns. In this work, we propose a template-based multi-modal generative mixture-model of imaging data and apply it to the problems of inlier/outlier pattern classification and image synthesis. Results on synthetic and patient data demonstrate that the proposed method is able to synthesise unseen data and accurately localise pathological regions, even in the presence of large abnormalities. It also demonstrates that the proposed model can provide accurate and uncertainty-aware intensity estimates of expected imaging patterns.

We present a framework for intrinsic comparison of surface metric structures and curvatures. This work parallels the work of Kurtek et al. on parameterization-invariant comparison of genus zero shapes. Here, instead of comparing the embedding of spherically parameterized surfaces in space, we focus on the first fundamental form. To ensure that the distance on spherical metric tensor fields is invariant to parameterization, we apply the conjugation-invariant metric arising from the L2 norm on symmetric positive definite matrices. As a reparameterization changes the metric tensor by a congruent Jacobian transform, this metric perfectly suits our purpose. The result is an intrinsic comparison of shape metric structure that does not depend on the specifics of a spherical mapping. Further, when restricted to tensors of fixed volume form, the manifold of metric tensor fields and its quotient of the group of unitary diffeomorphisms becomes a proper metric manifold that is geodesically complete. Exploiting this fact, and augmenting the metric with analogous metrics on curvatures, we derive a complete Riemannian framework for shape comparison and reconstruction. A by-product of our framework is a near-isometric and curvature-preserving mapping between surfaces. The correspondence is optimized using the fast spherical fluid algorithm. We validate our framework using several subcortical boundary surface models from the ADNI dataset.

Attenuation correction is an essential requirement for quantification of positron emission tomography (PET) data. In PET/CT acquisition systems, attenuation maps are derived from computed tomography (CT) images. However, in hybrid PET/MR scanners, magnetic resonance imaging (MRI) images do not directly provide a patient-specific attenuation map. The aim of the proposed work is to improve attenuation correction for PET/MR scanners by generating synthetic CTs and attenuation maps. The synthetic images are generated through a multi-atlas information propagation scheme, locally matching the MRI-derived patient's morphology to a database of MRI/CT pairs, using a local image similarity measure. Results show significant improvements in CT synthesis and PET reconstruction accuracy when compared to a segmentation method using an ultrashort-echo-time MRI sequence and to a simplified atlas-based method.

Objective: Minimally invasive autopsy using post-mortem magnetic resonance imaging (MRI) is a valid alternative to conventional autopsy in fetuses and infants. Estimation of brain weight is an integral part of autopsy, but manual segmentation of organ volumes on MRI is labor intensive and prone to errors, therefore unsuitable for routine clinical practice. In this paper we aim to show that volumetric measurements of the post-mortem fetal and neonatal brain can be accurately estimated using semi-automatic techniques and a high correlation can be found with the weights measured from conventional autopsy results., Methods: The brains of 17 newborn subjects, part of Magnetic Resonance Imaging Autopsy Study (MaRIAS), were segmented from post-mortem MR images into cerebrum, cerebellum and brainstem using a publicly available neonate brain atlas and semi-automatic segmentation algorithm. The results of the segmentation were averaged to create a new atlas, which was then used for the automated atlas-based segmentation of 17 MaRIAS fetus subjects. As validation, we manually segmented the MR images from 8 subjects of each cohort and compared them with the automatic ones. The semi-automatic estimation of cerebrum weight was compared with the results of the conventional autopsy., Results: The Dice overlaps between the manual and automatic segmentations are 0.991 and 0.992 for cerebrum, 0.873 and 0.888 for cerebellum and 0.819 and 0.815 for brainstem, for newborns and fetuses, respectively. Excellent agreement was obtained between the estimated MR weights and autopsy gold standard ones: mean absolute difference of 5 g and 2% maximum error for the fetus cohort and mean absolute difference of 20 g and 11% maximum error for the newborn one., Conclusions: The high correlation between the obtained segmentation and autopsy weights strengthens the idea of using post-mortem MRI as an alternative for conventional autopsy of the brain.

Background: Pathological abnormalities including demyelination and neuronal loss are reported in the outer cortex in multiple sclerosis (MS)., Objective: We investigated for in vivo evidence of outer cortical abnormalities by measuring the magnetisation transfer ratio (MTR) in MS patients of different subgroups., Methods: Forty-four relapsing-remitting (RR) (mean age 41.9 years, median Expanded Disability Status Scale (EDSS) 2.0), 25 secondary progressive (SP) (54.1 years, EDSS 6.5) and 19 primary progressive (PP) (53.1 years, EDSS 6.0) MS patients and 35 healthy control subjects (mean age 39.2 years) were studied. Three-dimensional (3D) 1×1×1mm(3) T1-weighted images and MTR data were acquired. The cortex was segmented, then subdivided into outer and inner bands, and MTR values were calculated for each band., Results: In a pairwise analysis, mean outer cortical MTR was lower than mean inner cortical MTR in all MS groups and controls (p<0.001). Compared with controls, outer cortical MTR was decreased in SPMS (p<0.001) and RRMS (p<0.01), but not PPMS. Outer cortical MTR was lower in SPMS than PPMS (p<0.01) and RRMS (p<0.01)., Conclusions: Lower outer than inner cortical MTR in healthy controls may reflect differences in myelin content. The lowest outer cortical MTR was seen in SPMS and is consistent with more extensive outer cortical (including subpial) pathology, such as demyelination and neuronal loss, as observed in post-mortem studies of SPMS patients., (© The Author(s) 2014.)

Multi-atlas segmentation propagation has evolved quickly in recent years, becoming a state-of-the-art methodology for automatic parcellation of structural images. However, few studies have applied these methods to preclinical research. In this study, we present a fully automatic framework for mouse brain MRI structural parcellation using multi-atlas segmentation propagation. The framework adopts the similarity and truth estimation for propagated segmentations (STEPS) algorithm, which utilises a locally normalised cross correlation similarity metric for atlas selection and an extended simultaneous truth and performance level estimation (STAPLE) framework for multi-label fusion. The segmentation accuracy of the multi-atlas framework was evaluated using publicly available mouse brain atlas databases with pre-segmented manually labelled anatomical structures as the gold standard, and optimised parameters were obtained for the STEPS algorithm in the label fusion to achieve the best segmentation accuracy. We showed that our multi-atlas framework resulted in significantly higher segmentation accuracy compared to single-atlas based segmentation, as well as to the original STAPLE framework.

The detection of abnormal intensities in brain images caused by the presence of pathologies is currently under great scrutiny. Selecting appropriate models for pathological data is of critical importance for an unbiased and biologically plausible model fit, which in itself enables a better understanding of the underlying data and biological processes. Besides, it impacts on one's ability to extract pathologically meaningful imaging biomarkers. With this aim in mind, this work proposes a fully unsupervised hierarchical model selection framework for neuroimaging data which permits the stratification of different types of abnormal image atterns without prior knowledge about the subject's pathological status.

Multiple Sclerosis lesions influence the process of image analysis, leading to tissue segmentation problems and biased morphometric estimates. With the aim of reducing this bias, existing techniques fill segmented lesions as normal appearing white matter. However, due to lesion segmentation errors or the presence of neighbouring structures, such as the ventricles and deep grey matter structures, filling all lesions as white matter like intensities is prone to introduce errors and artefacts. In this paper, we present a novel lesion filling strategy based on in-painting techniques for image completion. This technique makes use of a patch-based Non-Local Means algorithm that fills the lesions with the most plausible texture, rather than normal appearing white matter. We demonstrate that this strategy introduces less bias and fewer artefacts and spurious edges than previous techniques. The advantages of the proposed methodology are that it preserves both anatomical structure and signal-to-noise characteristics even when the lesions are neighbouring grey matter and cerebrospinal fluid, and avoids excess blurring or rasterisation due to the choice of segmentation plane, and lesion shape, size and/or position.

There is considerable interest in designing therapeutic studies of individuals at risk of Alzheimer disease (AD) to prevent the onset of symptoms. Cortical β-amyloid plaques, the first stage of AD pathology, can be detected in vivo using positron emission tomography (PET), and several studies have shown that ~1/3 of healthy elderly have significant β-amyloid deposition. Here we assessed whether asymptomatic amyloid-PET-positive controls have increased rates of brain atrophy, which could be harnessed as an outcome measure for AD prevention trials. We assessed 66 control subjects (age = 73.5±7.3 yrs; MMSE = 29±1.3) from the Australian Imaging Biomarkers & Lifestyle study who had a baseline Pittsburgh Compound B (PiB) PET scan and two 3T MRI scans ~18-months apart. We calculated PET standard uptake value ratios (SUVR), and classified individuals as amyloid-positive/negative. Baseline and 18-month MRI scans were registered, and brain, hippocampal, and ventricular volumes and annualized volume changes calculated. Increasing baseline PiB-PET measures of β-amyloid load correlated with hippocampal atrophy rate independent of age (p = 0.014). Twenty-two (1/3) were PiB-positive (SUVR>1.40), the remaining 44 PiB-negative (SUVR≤1.31). Compared to PiB-negatives, PiB-positive individuals were older (76.8±7.5 vs. 71.7±7.5, p<0.05) and more were APOE4 positive (63.6% vs. 19.2%, p<0.01) but there were no differences in baseline brain, ventricle or hippocampal volumes, either with or without correction for total intracranial volume, once age and gender were accounted for. The PiB-positive group had greater total hippocampal loss (0.06±0.08 vs. 0.02±0.05 ml/yr, p = 0.02), independent of age and gender, with non-significantly higher rates of whole brain (7.1±9.4 vs. 4.7±5.5 ml/yr) and ventricular (2.0±3.0 vs. 1.1±1.0 ml/yr) change. Based on the observed effect size, recruiting 384 (95%CI 195-1080) amyloid-positive subjects/arm will provide 80% power to detect 25% absolute slowing of hippocampal atrophy rate in an 18-month treatment trial. We conclude that hippocampal atrophy may be a feasible outcome measure for secondary prevention studies in asymptomatic amyloidosis.

The combination of functional and anatomical imaging technologies such as Positron Emission Tomography (PET) and Computed Tomography (CT) has shown its value in the preclinical and clinical fields. In PET/CT hybrid acquisition systems, CT-derived attenuation maps enable a more accurate PET reconstruction. However, CT provides only very limited soft-tissue contrast and exposes the patient to an additional radiation dose. In comparison, Magnetic Resonance Imaging (MRI) provides good soft-tissue contrast and the ability to study functional activation and tissue microstructures, but does not directly provide patient-specific electron density maps for PET reconstruction. The aim of the proposed work is to improve PET/MR reconstruction by generating synthetic CTs and attenuation-maps. The synthetic images are generated through a multi-atlas information propagation scheme, locally matching the MRI-derived patient's morphology to a database of pre-acquired MRI/CT pairs. Results show improvements in CT synthesis and PET reconstruction accuracy when compared to a segmentation method using an Ultrashort-Echo-Time MRI sequence.

X-ray mammography is routinely used in national screening programmes and as a clinical diagnostic tool. Magnetic Resonance Imaging (MRI) is commonly used as a complementary modality, providing functional information about the breast and a 3D image that can overcome ambiguities caused by the superimposition of fibro-glandular structures associated with X-ray imaging. Relating findings between these modalities is a challenging task however, due to the different imaging processes involved and the large deformation that the breast undergoes. In this work we present a registration method to determine spatial correspondence between pairs of MR and X-ray images of the breast, that is targeted for clinical use. We propose a generic registration framework which incorporates a volume-preserving affine transformation model and validate its performance using routinely acquired clinical data. Experiments on simulated mammograms from 8 volunteers produced a mean registration error of 3.8±1.6mm for a mean of 12 manually identified landmarks per volume. When validated using 57 lesions identified on routine clinical CC and MLO mammograms (n=113 registration tasks) from 49 subjects the median registration error was 13.1mm. When applied to the registration of an MR image to CC and MLO mammograms of a patient with a localisation clip, the mean error was 8.9mm. The results indicate that an intensity based registration algorithm, using a relatively simple transformation model, can provide radiologists with a clinically useful tool for breast cancer diagnosis., (Copyright © 2012 Elsevier B.V. All rights reserved.)

This work explores a fully-automated algorithm for estimation of the uptake of radio-pharmaceutical in brain MR-PET imaging. The algorithm is based on a model of the pharmaceutical uptake coupled with probabilistic models of the PET and MR acquisition systems. In contrast to algorithms that attempt to correct for the partial volume effect (PVE), the problem is tackled here in the reconstruction by means of a probabilistic model of the pharmaceutical uptake. We make use of hybrid Bayesian networks to describe the joint probabilistic model and to obtain an efficient optimisation algorithm. We describe solutions adopted in order to mitigate the effect of local maxima and to reduce the sensitivity to the initialisation of the parameters, rendering the algorithm fully automatic. The algorithm is evaluated on simulated MR-PET data and on the reconstruction of clinical PET FDG acquisitions.

Cortical thinning is a widely used and powerful biomarker for measuring disease progression in Alzheimer's disease (AD). However, there has been little work on the effect of atrophy on the cortical folding patterns. In this study, we examined whether the cortical folding could be used as a biomarker of AD. Cortical folding metrics were computed on 678 patients from the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort. For each subject, the boundary between grey matter and white matter was extracted using a level set technique. At each point on the boundary two metrics characterising folding, curvedness and shape index, were generated. Joint histograms using these metrics were calculated for five regions of interest (ROIs): frontal, temporal, occipital, and parietal lobes as well as the cingulum. Pixelwise statistical maps were generated from the joint histograms using permutations tests. In each ROI, a significant reduction was observed between controls and AD in areas associated with the sulcal folds, suggesting a sulcal opening associated with neurodegeneration. When comparing to MCI patients, the regions of significance were smaller but overlapping with those regions found comparing controls to AD. It indicates that the differences in cortical folding are progressive and can be detected before formal diagnosis of AD. Our preliminary analysis showed a viable signal in the cortical folding patterns for Alzheimer's disease that should be explored further.

Preterm birth increases the risk of perinatal brain injury and is believed to initiate a cascade of processes causing white matter damage resulting in subsequent neurological deficit; neonatal magnetic resonance imaging provides a number of potential biomarkers of this deficit. In this work we unify measures of the cortical folding pattern and of white matter integrity to establish correlation between grey and white matter derived properties. Diffusion weighted MRI has revealed that the cortical grey matter in the extremely preterm period exhibits a strong transient radial organisation suggesting neuronal axons are orientated towards the underlying white matter. This effect is lost during cortical maturation and is considered no longer visible on MRI at term equivalent age. Here we show that, in a group of 19 infants, radial organisation in the cortical grey matter remains detectable at term-equivalent age and that there is a strong anterior-posterior asymmetry. A group of three infants with moderate or severe abnormal white matter abnormality have significantly higher cortical grey matter radial organisation (p < 0.02), higher grey matter FA (p < 0.01) and a lower measure of cortical complexity (p < 0.03) than infants with normal or mild abnormal white matter abnormality; all measures associated with the preterm phenotype before term equivalent age. The novel combination of state-of-the-art imaging techniques, analysing grey-matter based spatial characteristics, may provide insight into the mechanism of neurodevelopmental deficits seen in infants with abnormal MR imaging at term equivalent age.

Homogenising the availability of manually generated information in large databases has been a key challenge of medical imaging for many years. Due to the time consuming nature of manually segmenting, parcellating and localising landmarks in medical images, these sources of information tend to be scarce and limited to small, and sometimes morphologically similar, subsets of data. In this work we explore a new framework where these sources of information can be propagated to morphologically dissimilar images by diffusing and mapping the information through intermediate steps. The spatially variant data embedding uses the local morphology and intensity similarity between images to diffuse the information only between locally similar images. This framework can thus be used to propagate any information from any group of subject to every other subject in a database with great accuracy. Comparison to state-of-the-art propagation methods showed highly statistically significant (p < 10(-4)) improvements in accuracy when propagating both structural parcelations and brain segmentations geodesically.