User-initialized active contour segmentation and golden-angle real-time cardiovascular magnetic resonance enable accurate assessment of LV function in patients with sinus rhythm and arrhythmias.

BackgroundData obtained during arrhythmia is retained in real-time cardiovascular magnetic resonance (rt-CMR), but there is limited and inconsistent evidence to show that rt-CMR can accurately assess beat-to-beat variation in left ventricular (LV) function or during an arrhythmia.MethodsMulti-slice, short axis cine and real-time golden-angle radial CMR data was collected in 22 clinical patients (18 in sinus rhythm and 4 patients with arrhythmia). A user-initialized active contour segmentation (ACS) software was validated via comparison to manual segmentation on clinically accepted software. For each image in the 2D acquisitions, slice volume was calculated and global LV volumes were estimated via summation across the LV using multiple slices. Real-time imaging data was reconstructed using different image exposure times and frame rates to evaluate the effect of temporal resolution on measured function in each slice via ACS. Finally, global volumetric function of ectopic and non-ectopic beats was measured using ACS in patients with arrhythmias.ResultsACS provides global LV volume measurements that are not significantly different from manual quantification of retrospectively gated cine images in sinus rhythm patients. With an exposure time of 95.2 ms and a frame rate of > 89 frames per second, golden-angle real-time imaging accurately captures hemodynamic function over a range of patient heart rates. In four patients with frequent ectopic contractions, initial quantification of the impact of ectopic beats on hemodynamic function was demonstrated.ConclusionUser-initialized active contours and golden-angle real-time radial CMR can be used to determine time-varying LV function in patients. These methods will be very useful for the assessment of LV function in patients with frequent arrhythmias

Efficient Generation of Shape-Based Reference Frames for the Corpus Callosum for DTI-based Connectivity Analysis

Yushkevich et.al. [17, 18] established a PDE-based deformable modeling approach called continuous medial representation (cm-rep), in which the geometric relationship between the medial axis of a 3D object and its boundary is captured. Continuous medial description of an object not only provides useful shape features for object characterization and comparison; it also imposes a shape-based reference frame on the interior of that object. Such a reference frame provides a useful means of representing different instances of an anatomical structure using a common canonical parametrization domain. This paper presents an efficient method to construct continuous medial shape models for 2D objects. A closed form solution for the ordinary differential equation (ODE) is derived via Pythagorean hodograph (PH) curves. That closed form solution reduces the computation complexity from solving an ODE system to pure algebraic manipulation. Using this method, we generate shape-based reference frames, and demonstrate how they can be applied to the analysis of anatomical connectivity of corpora callosa, obtained by fiber tracking in diffusion tensor magnetic resonance imaging (DTI) in a chromosome 22q11.2 deletion syndrome study

NODEO: A Neural Ordinary Differential Equation Based Optimization Framework for Deformable Image Registration

Deformable image registration (DIR), aiming to find spatial correspondence between images, is one of the most critical problems in the domain of medical image analysis. In this paper, we present a novel, generic, and accurate diffeomorphic image registration framework that utilizes neural ordinary differential equations (NODEs). We model each voxel as a moving particle and consider the set of all voxels in a 3D image as a high-dimensional dynamical system whose trajectory determines the targeted deformation field. Our method leverages deep neural networks for their expressive power in modeling dynamical systems, and simultaneously optimizes for a dynamical system between the image pairs and the corresponding transformation. Our formulation allows various constraints to be imposed along the transformation to maintain desired regularities. Our experiment results show that our method outperforms the benchmarks under various metrics. Additionally, we demonstrate the feasibility to expand our framework to register multiple image sets using a unified form of transformation,which could possibly serve a wider range of applications

Structure specific analysis of the hippocampus in temporal lobe epilepsy

The hippocampus is a major structure of interest affected by temporal lobe epilepsy (TLE). Region of interest (ROI)-based analysis has traditionally been used to study hippocampal involvement in TLE, although spatial variation of structural and functional pathology have been known to exist within the ROI. In this article, structure-specific analysis (Yushkevich et al. (2007) Neuroimage 35:1516–1530) is applied to the study of both structure and function in TLE patients. This methodology takes into account information about the spatial correspondence of voxels within ROIs on left and right sides of the same subject as well as between subjects. Hippocampal thickness is studied as a measure of structural integrity, and functional activation in a functional magnetic resonance imaging (fMRI) experiment in which subjects performed a memory encoding task is studied as a measure of functional integrity. Pronounced disease-related decrease in thickness is found in posterior and anterior hippocampus. A region in the body also shows increased thickness in patients' healthy hippocampi compared with controls. Functional activation in diseased hippocampi is reduced in the body region compared to controls, whereas a region in the tail showing greater right-lateralized activation in controls also shows greater activation in healthy hippocampi compared with the diseased side in patients. Summary measurements generated by integrating quantities of interest over the entire hippocampus can also be used, as is done in conventional ROI analysis. © 2009 Wiley-Liss, Inc.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/63055/1/20620_ftp.pd

Regional Deep Atrophy: a Self-Supervised Learning Method to Automatically Identify Regions Associated With Alzheimer's Disease Progression From Longitudinal MRI

Longitudinal assessment of brain atrophy, particularly in the hippocampus, is a well-studied biomarker for neurodegenerative diseases, such as Alzheimer's disease (AD). In clinical trials, estimation of brain progressive rates can be applied to track therapeutic efficacy of disease modifying treatments. However, most state-of-the-art measurements calculate changes directly by segmentation and/or deformable registration of MRI images, and may misreport head motion or MRI artifacts as neurodegeneration, impacting their accuracy. In our previous study, we developed a deep learning method DeepAtrophy that uses a convolutional neural network to quantify differences between longitudinal MRI scan pairs that are associated with time. DeepAtrophy has high accuracy in inferring temporal information from longitudinal MRI scans, such as temporal order or relative inter-scan interval. DeepAtrophy also provides an overall atrophy score that was shown to perform well as a potential biomarker of disease progression and treatment efficacy. However, DeepAtrophy is not interpretable, and it is unclear what changes in the MRI contribute to progression measurements. In this paper, we propose Regional Deep Atrophy (RDA), which combines the temporal inference approach from DeepAtrophy with a deformable registration neural network and attention mechanism that highlights regions in the MRI image where longitudinal changes are contributing to temporal inference. RDA has similar prediction accuracy as DeepAtrophy, but its additional interpretability makes it more acceptable for use in clinical settings, and may lead to more sensitive biomarkers for disease monitoring in clinical trials of early AD.Comment: Submitted to NeuroImage for revie

In-vivo heterogeneous functional and residual strains in human aortic valve leaflets

Residual and physiological functional strains in soft tissues are known to play an important role in modulating organ stress distributions. Yet, no known comprehensive information on residual strains exist, or non-invasive techniques to quantify in-vivo deformations for the aortic valve (AV) leaflets. Herein we present a completely non-invasive approach for determining heterogeneous strains – both functional and residual – in semilunar valves and apply it to normal human AV leaflets. Transesophageal 3D echocardiographic (3DE) images of the AV were acquired from open-heart transplant patients, with each AV leaflet excised after heart explant and then imaged in a flattened configuration ex-vivo. Using an established spline parameterization of both 3DE segmentations and digitized ex-vivo images (Aggarwal et al., 2014), surface strains were calculated for deformation between the ex-vivo and three in-vivo configurations: fully open, just-coapted, and fully-loaded. Results indicated that leaflet area increased by an average of 20% from the ex-vivo to in-vivo open states, with a highly heterogeneous strain field. The increase in area from open to just-coapted state was the highest at an average of 25%, while that from just-coapted to fully-loaded remained almost unaltered. Going from the ex-vivo to in-vivo mid-systole configurations, the leaflet area near the basal attachment shrank slightly, whereas the free edge expanded by ~10%. This was accompanied by a 10° −20° shear along the circumferential-radial direction. Moreover, the principal stretches aligned approximately with the circumferential and radial directions for all cases, with the highest stretch being along the radial direction. Collectively, these results indicated that even though the AV did not support any measurable pressure gradient in the just-coapted state, the leaflets were significantly pre-strained with respect to the excised state. Furthermore, the collagen fibers of the leaflet were almost fully recruited in the just-coapted state, making the leaflet very stiff with marginal deformation under full pressure. Lastly, the deformation was always higher in the radial direction and lower along the circumferential one, the latter direction made stiffer by the preferential alignment of collagen fibers. These results provide significant insight into the distribution of residual strains and the in-vivo strains encountered during valve opening and closing in AV leaflets, and will form an important component of the tool that can evaluate valve׳s functional properties in a non-invasive manner

Multimodal image analysis and subvalvular dynamics in ischemic mitral regurgitation

Background: The exact geometric pathogenesis of leaflet tethering in ischemic mitral regurgitation (IMR) and the relative contribution of each component of the mitral valve complex (MVC) remain largely unknown. In this study, we sought to further elucidate mitral valve (MV) leaflet remodeling and papillary muscle dynamics in an ovine model of IMR with magnetic resonance imaging (MRI) and 3-dimensional echocardiography (3DE). Methods: Multimodal imaging combining 3DE and MRI was used to analyze the MVC at baseline, 30 minutes post–myocardial infarction (MI), and 12 weeks post-MI in ovine IMR models. Advanced 3D imaging software was used to trace the MVC from each modality, and the tracings were verified against resected specimens. Results: 3DE MV remodeling was regionally heterogenous and observed primarily in the anterior leaflet, with significant increases in surface area, especially in A2 and A3. The posterior leaflet was significantly shortened in P2 and P3. Mean posteromedial papillary muscle (PMPM) volume was decreased from 1.9 ± 0.2 cm3 at baseline to 0.9 ± 0.3 cm3 at 12 weeks post-MI (P <.05). At 12 weeks post-MI, the PMPM was predominately displaced horizontally and outward along the intercommissural axis with minor apical displacement. The subvalvular contribution to tethering is a combination of unilateral movement, outward displacement, and degeneration of the PMPM. These findings have led to a proposed new framework for characterizing PMPM dynamics in IMR. Conclusions: This study provides new insights into the complex interrelated and regionally heterogenous valvular and subvalvular mechanisms involved in the geometric pathogenesis of IMR tethering